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FOREWORD 

This report summarizes the work conducted by McDonnell Douglas Astronautics 
Company-East (MDAC-E) in St, Louis, Missouri for the Structures and Mechanics 
Division of the NASA Johnson Space Center (NASA-JSC) under Contract NAS9-14012, 

"Data Correlation and Analysis of Arc Tunnel and Wind Timnel Test of RSI Joints 
and Gaps, Phase II." This final report consists of two volumes: Volume I - 
Technical Report and Volume II - Data Base, Part 1 and Part 2. The period of perfor- 
mance was from 20 May 1974 thru 19 May 1975. 

Mr. Donald J. Tillian was the NASA Technical Monitor for this study; Messrs. 

H. E. Christensen and H. W. Kipp were the MDAC Principal Investigator and Study 
Manager, respectively. Significant contributions to this study were made by 
A. E. Bruns, M. B. Donovan, L. H. Ebbesmeyer, E. A. Eiswirth and T. W. Parkinson. 

The cooperation of numerous NASA Personnel at Ames Research Center, Johnson Space 
Center and Langley Research Center in providing experimental data, supplemental 
calculations and valuable counsel was instrumental to the successful completion 
of this study. Data used in Section 4.2 are based on an analysis originally 
performed for Rockwell International Corp. under the direction of P. C. Merhoff. 

We are appreciative of the cocporation from the following for supplying test data; 

C. D. Scott and L. P. Murray of JSC, W. K. Lockman and F. J. Centolanzi of Ames, 

D. A. Throckmoi ton and I. Weinstein of LaRC, and G. W. Mauss and C. B. Blumer 
of Rockwell International. Special acknowledgement is made of assistance and 
support provided by the NaSA Ames Research Center under Contract NAS2— 7897 (Rev 4) 
for the RSI sizing calculations using current Shuttle design philosophy which 

are presented in Section 7 of this report. 

The International System of Units is used as the primary system for all 
resulus reported herein. The results are also reported in the British Engineering 
System of Units which was used for calculations made during the course of this 
study. 
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ABSTRACT 

Heat transfer data measured In gaps representative of those being employed 
for joints In the Space Shuttle RSI thermal protection systemr have been assimilated, 
analyzed and correlated. The study reported herein Is the second phase of an earlier 
study of gap heating phenomena reported In Reference 1. Portions of these results 
are Included herein for completenes>^ . 

The body of data under study was obtained In six NASA facilities, the Ames 
3.5 Foot Hypersonic Wind Tunnel, the Ames 20 Megawatt Turbulent Duct, the JSC 10 
Megawatt Arc Tunnel, the LaRC Mach 10 Continuous Flow Hypersonic Tunnel, the LaRC 
Mach 8 Variable Density Tunnel, and the LaRC 8 Foot High Temperature Structures 
Tunnel. Several types of gap were Investigated with emphasis on simple butt joints. 
Gap widths ranged from 0.0 to 0.76 cm and depths ranged from 1 to 6 cm. Laminar, 
transitional and turbulent boundary layer flows over the gap opening were investi- 
gated. The angle between gap axis and external flow was varied between 0 and tt/ 2 
radians. The "contoured" cross section gap performed significantly better than all 
other wide gaps and slightly better than all other narrow gap geometries. Three- 
dimensional heating variations were observed within gaps in the absence of external 
flow pressure gradients. Interactions between heating within gaps and heating of 
adjacent top tile surfaces were observed in several tests. Also, gaps aligned with 
the flow were observed to promote boundary layer transition in tests conducted in 
the Ames 3.5 Foot Hypersonic Wind Tunnel. Heat transfer correlation equations were 
obtained for many of the tests. TPS thickness requirements with and without gaps 
were computed for a current Shuttle entry trajectory. Experimental data employed in 
the study are summarized in Volume II. 
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The reusable surface insulation (RSI)* thencal protection system (TPS) for 
Space Shuttle requires gaps at RSI Joints to accommodate structural deflection 
resulting from loads and thermal expansion. In addition, allowance must be made 
for manufacturing and assembly tolerances. At room temperature, gap widths under 
current consideration range from 0.127 cm to 0.254 cm (0.050" to 0.10") ^ 0.038 cm to 
0.076 cm (0.015" to 0.030”). In orbital operation, these may shrink to near zero 
during cold soak or grow by as much as 25%. Candidate tile edge radii range from 
0.076 to 0.254 cm (0.030" to 0.1"). 

The successful application of RSI material for Shuttle thermal protection is 
significantly affected by entry heating within the RSI gaps. Gap width, depth, 
cross section geometry, gap orientation, boundary layer state and surface mismatch 
are all known to affect convective heating within the gap and heat leakage to the 
protected substructure. For instance, present study results indicate a 0.254 cm 
wide flush transverse butt gap increases TPS thickness requirements by approximately 
33% above the thickness required for a TPS without gaps. 

During 1972 and 1973, extensive tests of various gap configurations were run 
by NASA to provide a data base for accurate assessment of gap heating were 

taken in both wind tunnels and in arc tunnels. A large segment of these data were 
analyzed and correlated to obtain methods for predicting heating in RSI gaps on 
Shuttle. Based on the correlations, the effect of gap heating on Shuttle lower surface 
TPS requirements was determined for a typical reentry trajectory. The results of 
these Phase I studies are documented in Reference 1 . 

Subsequently, additional tests were run by NASA to clarify prior data and to 
address unresolved questions. The analysis of these data on a consistent basis with 
the earlier analysis was the objective of the Phase II study documented herein. 

The test data sources utilized in this study are identified in Figure 1. Data sources 
used in the earlier Phase I study are also shown, identified with an asterisk. 

NASA is continuously Improving and adding test facilities. Tlroughout this report 
we refer to testing in the JSC 10 Ifcj Arc Tunnel facility using a laminar duct con- 
figuration. This duct was operated at TP2 (test position number 2). Similarly, the 
NASA personnel at Ames have several duct test facilities, and the test programs 
evaluated in this study were performed using the 2"x9" turbulent duct facility. 


* Also referred to as HRSI (High Temperature Reusable Surface Insulation) 
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The following major tasks were performed during the study: 
o Assimilation of gap heating data from NASA facilities 

o Analysis of the data to determine heating rates and sensitivities; comparison 
of the various candidate joints 

o Correlation of the assimilated data e^nd the development of a gap heating 
procedure which was applied to a current Shuttle trajectory. 

This volume describes the assimilation, analyses and correlations resulting from 
the study as well as the conclusions derived therefrom. Volume II of this report 
presents the basic gap heating data Including information relating to each test 
facility, run schedule, test conditions and models. 
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2.0 SUMMARY 

This study of heat transfer within RSI gaps builds upon the Phase I study re- 
sults documented in Reference 1 . As in Phase I, the Phase II work described herein 
was performed in three tasks, namely, data assimilation, data analysis and correla- 
tion of results. Highlights of these tasks are summarized in this section. 

The data assimilation task entailed compilation and evaluation of the gap 
heating data from the new sources listed in Figure 1. Also, errors in some of the 
thin skin heat transfer data assimilated during Phase I were eliminated by a new 
data reduction employing corrected specific heat values. In addition, skin conduc- 
tion corrections were established for portions of the data. The experimental data 
were read from magnetic tapes provided by the various test facilities and trans- 
C'ibed into a uniform format in a gap heating data bank. To facilitate data re- 
trieval and analysis, 24 attribute words were assigned to each heating data point. 
These attribute words consisted of information such as test and geometry identifiers, 
instrumentation locations, flow orientation, invlscld flow conditions and boundary 
layer parameters. The SELECT computer program, prepared during the Phase I study 
was used to assess the data bank and to prepare specific data for the subsequent 
multiple regression analysis. The assimilated data are compiled in a test data 
document which is the second volume of this report. 

The data analysis task entailed performing a mumber of diverse subtasks. In- 
cluded were reduction of temperature histories measured on RSI tiles in the JSC 
10 MW Tunnel and in the LaRC 8 Foot HTSf to heat flux by means of an inverse solu- 
tion. Subsequently these data were incorporated in the data bank. Additional sub- 
tasks included graphic data presentation, data-theory comparison, sensitivity 
analysis and boundary layer calculations. End results of the analysis task in- 
cluded the identification of significant phenomena observed in the test program 
and the preparation of data for correlation. 

Reinforcement of many of the Phase I conclusions resulted from the additional 
data obtained during Phase II. 

Among the more significant observations made during Phase II is that gaps 
aligned with the flow promote transition. Examination of the data taken in the 
Ames 3.5 Foot H.W.T. indicates that the transition onset occurs within the gap. 

The parametric variation of edge radius tested in the JSC 10 MW indicates that 
increased convective heating results from Increasing edge radius. Coupling these 
data with a complete RSI tile heat transfer analysis Indicates, however, that the 
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Incrasscd convective heating Is more than offset by Increased radiation r<- 
obtalned from the larger radius. 

Correlation equations for several new classes of gap heating were obtained 
during Phase II and Improvements were made on several developed during Phase I. 

The development of correlation equations benefited significantly by the use oi mul- 
tiple regression analysis. The correlation of transverse gap data obtained during 
Phase I has been updated by Including a larger body of data. An Improved correla- 
tion of In-line gaps to Include flow Incidence up to 51 degrees was developed using 
data from the Ames 3.5 Foot HWT. A correlation for the effect of edge radius on 
In-line and transverse gap heating was developed using data obtained In the JSC 10 
MW Tunnel. 

An automated gap heating subroutine Incorporating the above correlations has 
been developed. This subroutine Is described In Section 6 of this report. 

A re-assessment of the Influence of gap heating on TPS requirement was made 
during Phase II. The assessment Included a parametric evaluation of gap width and 
radius effect. The Increase In TPS thickness caused by the presence of gaps was 
found to vary between 12% and 42%. The calculations are described In Section 7. 
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3. TEST PROGRAMS AND DATA ASSIMILATION 

Ten test programs constitute the principal data '''\ae for this study as 
summarized in Figure 1 . Four of these ten programs were previously reported in 
detail in Reference 1 , but are summarized herein for reader convenience. In 
addition, supplenental runs were made for the previous tests conducted in the LaRC 
Mach 10 Continuois Flow Hypersonic Tunnel (CFHT) . The primary data assimilated 
during Phase II resulted from test programs conducted in the JSC Laminar Duct, the 
LaRC 8 foot High Temperature Structures Tunnel (HTST) , the Ames 20 MW Turbulent 
Duct, the JSC 10 MW Laminar Duct and the Ames 3.5 foot Hypervelocity Wind Tunnel 
(HWT) facility. 

A general description of each test program and the data assimilated follows. 
Additional infomution related to the test facilities, models, test conditions and 
data appears in Volume II of this report and Volume II of Reference 1 (Phase I 
data) . The terminology used in this report to describe gap configurations and tile 
arrangements is depicted in Figure 2 . 

3.1 Data Previously Assimilated (Phase I) - Gap heating tests were peT-*' ed 
in the channel nozzle of the JSC 10 MW Arc Tunnel to provide heating dat e 

presence of a high enthalpy laminar boundary layer. The tests employed .. . of 

Mullite RSI tile which were heavily instrumented on gap surfaces and in depth at 
the center of one tile in Cach array. Thirteen models employing a variety of gap 
and tile configurations were tested. The gaps between the tiles were adjustable 
to study the effects of gap width using consistent sets of instrumentation. Four 
gap settings were employed (0.127, 0.254, 0.381 and 0.762 cm) with tile thickness 
of 3.175, 5.08 and 6.35 cm. butt joint step heights of j<).381 cm wore also in- 
vestigated. Figure 3 summarizes the matrix of configurations tested in the 
channel nozzle. The figure also depicts the arc tunnel and shows a 5.08 cm butt 
joint model. Temperature response data from 36 gap locations were analyzed to 
obtain comparisons of heating rate distributions using the inverse solution tech- 
nique . 

Gap heating tests were conducted in the LaRC Mach 10 CFHT to provide data in 
the presence of a relatively thick turbulent boundary layer. The tests employed a 
wall-mounted stainless steel, thin skin "tile" model in which the instrumented thin 
skin tile was surrounded by an array of un instrumented RSI tiles, "he model was 
mounted on a turntable to permit variation of flow orientation relative to the tile 
array. A total of 157 runs was made tt a unit Reynolds number per meter of 
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3.2x10^. The test matrix included nine flow orientations (0 to tt/2 radians), four 
gap widths (0.127, 0.229, 0.457 and 0.711 cm) and three step heights (0, +0.254 and 
-0.168 cm). Tile thickness was 6.35 cm. Figure 4 shows the instrumented tile, 
the surrounding array of uninst rumen ted RSI tiles and the installation in the CFHT 
tunnel wall. A photograph of the installation is shown in Figure 5 which also 
summarizes the test matrix. Temperature response data from 81 thermocouples were 
analyzed to obtain heat fluxes and to determine the effects of inline versus stag- 
gered tiles, gap width, flow orientation and steps. 

Gap heating tests were conducted by C. B. Johnson (NASA, Langley Research Center) 
in the LaRC Ilach 8 Variable Density Tunnel (VDT) to provide data in the presence of 
laminar and turbulent boundary layers. The tests employed models which simulated 
thin skin tiles which were mounted in a curved plate. Models were tested both in 
the free stream and mounted flush with the tunnel wall. In each model position, the 
test section unit Reynolds number was varied over the range of 1.1x10^ to 40x10^ 
per meter. Both in-line and staggered tile configurations were evaluated at gap 
widths of 0.159, 0.317 and 0.476 cm. The tile geometry is shown in Figure 6 . 
Temperature data from 22 test runs was assimilated and analyzed to obtain heating 
patterns on in-line and staggered tiles. 

Gap heating tests were conducted by W. K. Lockman (NASA, Ames Research Center) 
and C. B. Blumer (Rockwell International) in the Ames 3.5 Foot H.W.T. 
facility to provide data in the presence of laminar, transitional, and turbulent 
boundary layers. The model consisted of a 68.6x152.4 cm carrier plate into which 
61.0x106.7 cm instrumented thin skin test articles were inserted. The five insert 
configurations used in the test program are shown in Figure 7 . The configurations 
included a flat calibration plate, a single transverse gap with and without surface 
steps, multiple transverse gaps, staggered tiles, and skewed intersecting gaps. 

Five gap settings (0, 0.064, 0.127, 0.254 and 0.508) were investigated with tile 
thickness of 1, 2 and 4 cm. Butt joint step heights of ^.159 and jHD.318 cm were 
also investigated. Figure 7 summarizes the variables tested in the HWT facility. 
Dajta from 81 test runs were assimilated and analyzed for gap heating distributions. 

Analysis of the above data is documented in Sf-'.tion 4.1 and correlation in 
Section 5. The facility, model description, test conditions and data are documented 
in Volume II of Reference 1 . 

3.2 Silica RSI Tile Tests in the JSC Laminar Duct - A test program was conducted 
by G. Mauss (Rockwell International) in the JSC 10 MW Arc Tunnel on test panels 
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using Silica RSI klles fabricated and Instrumented by NASA Ames. The purpose of 
the program was to obtain comparative thermal performance data on overlap and butt 
joint designs and the effect of tile thickness and gap width under simulated laminar 
flow conditions. Test panel overall dimensions were 33x33 cm. 

Six runs from the test program were selected for analysis. All six runs were 
made with the staggered gap design (shown in Figure 8 ) consisting of 3 full and 
2 half tiles with a perimeter cf narrow guard tiles. Gap widths of 0.127 and .254 
cm were evaluated with RSI thicknesses of 2.54 and 5.08 cm. Four of the six runs 
were made on butt joint models and two runs on overlap joint models, both of which 
are shown in cross-section in Figure 9 which also includes gap thermocouple loca- 
tions. Thermocouple location and distances were measured from x-ray photographs 
of the instrumented tiles. 

Test conditions for the six runs analyzed are summarized in Figure 10 . Arcjet 
mass flow rate was 45.4 grams/sec in all cases and enthalpy varied from 8.58 x 10° to 
21.95 X 10^ J/Kg. All runs were conducted in the channel nozzle of the JSC 10 MW 
Arc Jet Facility. The test panels were moxmted in one wall of the channel nozzle 
and a calibration plate was mounted on the other wall. The calibration plate served 
as a basis for referencing data and as a basis for investigating anomalies. Test 
results used in the analysis consisted of gap and plug temperature histories, test 
section pressures and channel wall temperatures. Figure 11 is a typical set of 
gap temperature histories obtained from the JSC automated plotting output. Data 
were also recorded on magnetic tapes for direct reading into the CDC 6500 computer 
to obtain heating rate distributions via the inverse solution technique. 

Analysis of the above data is reported in Section 4.2. The facility, model 
description, test conditions and data are documented in Volur.e II of this report. 

3.3 Supplemental LaRC Mach 10 CFHT Tests - Additional gap heating tests were 
conducted in the LaRC Mach 10 CFHT to provide supplemental data to the previous tests 
reported in Section 3.1. The purpose of these additional tests was to substantiate 
the calibration data previously measured in the CFHT facility, to evaluate the effect 
of zero gap width, and to determine the effects on heating in the gaps between the 
tiles at aigher Reynolds number. 

The test program employed the same model as used in the previous tests (see 
Figure 4 ), namely, a wall-mounted thin skin "tile" model In which the instrumented 
thin skin tile was surrounded by an array of uninstrumented RSI tiles. Again, the 
model was mounted on a turntable to permit variation of flow orientation relative 
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TEST MATRIX 

SILICA RSI TILE TESTS IN JSC lAMINAR DUCT 


RUN 

NO. 

MODEL 

GAP 

WIDTH, cm 

m 

grams/sec 

"■•'■'.'•NT, 

ENTHALPY, 

J/kg 

528 

2.54 cm OVERLAP 

0.127 

45.4 


8.582X10^ 

533 

5.08 cm OVERLAP 

0.127 

45.4 

480 

21.950X10® 

533 

2.54 cm BUn 

0.127 

45.4 

375 

11.578X10® 

539 

2.54 cm BUTT 

0.254 

45.4 

375 

11.578X10® 

542 

5.08 cm BUTT 

0.127 

45.4 

480 

21.950X^0® 

544 

5.08 cm BUTT 

0.254 

45.4 

480 

21.950X106 
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to the tile array. Results from ten runs were assimilated into the Data Bank for 
analysis. The test runs were conducted at unit Reynolds number per meter of both 
3.28 X 10^ and 7.38 x 10^. The test matrix included flow orientations of 0 and 
v/2 radians, gap widths of 0, 0.127, and 0.229 cm, and gap depths of 0.0 and 6.35 
cm. Both in-line and staggered butt joint designs were investigated. Figure 12 
summarizes the test matrix. 

Temperature response data from these runs was analyzed to obtain heating rate 
distributions using the inverse solution technique as discussed in Section 4.3. 

The .acility, model description, test conditions and data are documented in Volume 
II of this report. 

3.4 Edge Radius Tests at NASA JSC - Gap heating tests were conducted by C. D. 
Scott (NASA, Johnson Space Center) in the JSC 10 MW Arc Tunnel to provide heating 
data in the presence of a high enthalpy laminar boundary layer. The primary purpose 
of the tests was to investigate the effect of tile edge radius on gap heating. 

The tests employed sets of thin skin metallic tiles mounted in a wedge test 
fixture. Four edge radii (0.157, 0.3175, 0.635, and 1.27 cm) were parametrically 
tested at gap widths of 0.127, 0.254, and 0.381 cm. It should be pointed out that 
the tile wi*-h the small edge radius was originally specified as a sharp edged tile 
and many of the correlations were initially developed using this specification. 

The joint configuration was an in-line butt, and the tile height was 4.1275 cm 
for all tests. The test panel configuration, edge radii and instrumentation are 
shown in Figure 13 . 

Test data were received in tabular form. The analysis of the heating distri- 
butions is reported in Section 4.4 and correlation of the data in Section 5.3. 
Correlations for edge radius effects have been developed for both the vertical 
wall, and the edge-and-wall geometries. Test related information and data are 
documented in Volume II of this report. 


3.5 Single In-Line Gap Tests at Ames 3.5 Foot H.W.T. - Gap heating tests were 
conducted by W. K. Lockman (NASA, Ames Research Center) and C. B. Blumer (Rockwell 
International) in the Ames 3.5 foot H.W.T. facility to obtain data in the presence 
of laminar, transitional and turbulent boundary layers. These tests were conducted 
on panels furnished by Rockwell International and provided additional data to the 
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CFHT TEST 92 - SUPPLEMENTAL RUNS - FEBRUARY 1974 


RUN 

GAP (cm) 

ORIENTATION 

Re /m 
00 

161 


.27 
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3.28x10® 

162 



Staggered 

7.38x10® 

163 



In-line 

3.28x10® 

164 

.127 

In-line 

7.38x10® 

165 

.2 

29 

Staggered 

3.28x10® 

166 



Staggered 

7.38x10® 

167 



In-line 

3.28x10® 

168 

.2 

29 

In-line 

7.38x10® 

169 

Closed 


7.38x10® 

170 

Closed 


3.28x10® 
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previous tjests at this facility described in Section 3.1. The primary tes: pMrpose 
was to investigate the effects of tile orientation for several gap width bettii.,.^s . 

The t^sts employed a thin skin model Inserted into a 68.6 x 152.4 t > c^rr^er 
plate. The joint configuration was a single, in-line, gap oi 30.48 and 101.6 cm 
length. Thermocouples were installed along the top of the panel and along the 
faces of the gap. The model is shown in Figure 14 . The test matrix included 
four orientations (0, 5, 10 and 15 degrees), three gap widths (0.063,. 0.127 and 
0.254 cm), four gap depths (0, 1.016, 2.037 and 3.81 cm), and two gap lengths 
(30.48 and 101.6 cm). Test nms were conducted at Reynolds numbers per meter of 
1.64 X lO'’, 3.28 X 10 and 6.56 x 10 . The complete matrix of 93 runs is listed 
in Figure 15 . 

Information from these runs was assimilated into the Data Bank. Analysis of 
the heating rate distributions are described in Section 4.5. Correlations for long 
in-line gaps are contained in Section 5.3 for both zero flow orientation and 
orientation up to 15 degrees. Information relating to the test and data are 
documented in Volume II of this report. 

3.6 Test of Large Gap Panel in LaRC 8 Foot HTST - Gap heating tests were per- 
formed in the LaRC 8 Foot HTST to obtain heating data on a large gap panel in the 
presence of a turbulent boundary layer. The same test panel was to be used in both 
the HTST and the AFFDL 50 MW Arc Tunnel Tests. The panel was originally scheduled 
for testing at AFFDL but was switched to the LaRC HTST due to fabrication difficul- 
ties with the test fixture by the vendor. The 50 MW tests were Intended to deter- 
mine the effects of higher enthalpy on the heating in a field of RSI gaps, and would 
have also provided a comparison of arc and wind tunnel gap heating data. The AFFDL 
test program was su'.isequently cancell* d. The test program at LaRC was under the 
direction of I. Weinstien. 

The test panel as shown in Figure 16 consisted of eleven LI 900 silica tiles 
with an interchangeable thin skin metallic ceiiter tile. Testing was planned for 
both the center LI 900 tile and the center metallic tile, but time permitted testing 
only of the metallic tile. The panel size was 46 x 46 x 6.5 cm. 

Both the LI 900 and the metallic tiles were heavily instrumented with the thermo- 
couple locations for the thin wall metallic center tile shown in Figure 17 . The 
RSI panel was mounted in a large test sled for free stream testing and permitted 
variation of angle of attack relative to the tile array. The panel was tested in 
both the in-line and staggered tile configurations. Gaps between the tiles were 
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MODEL OH-43 

TEST SUPERVISED BY C. B. BLUMER 


NO. (X10°/m) 


GAP 

WIDTH 

(cm) 

GAP 

DEPTH 

(cm) 

GAP 

LENGTH 

(cm) 

GAP 

ORIENT. 

ANGLE 

(DEG) 

T/C 

SCHEDULE 

MODULE 

POSITION 

SCHEDULE 

0.254 

1 .016 

101.6 

0 

3A 
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adjustable to study the effects of gap width, and shims were used under the tiles 
to study the effects of step heights. A companion plate that fitted into the same 
opening as the test panel was used for calibration purposes. 

The test r>«atrlx for the 8 Foot HIST program is summarized in Figure 18 . 

Five gap settings were employed (0, 0.10, 0.18, 0.30, and 0.41 cm) with the tile 

thickness of 6.35 cm. Step heights Investigated were 0 and -fO.254 cm. Tests were 

run at Reynolds number per meter of 1.9 x 10^ >‘nd 4.8 x 10^, while the test sled 

angle of attack was varied from 0 to 15 degrees. - 

Problems specifically related to the model and facility occurred during the 
program. Figure 19 is a view of the panel prior to the first test and Figure 
20 shows the panel following the same test. The time in the stream was 11 seconds 
and coating was removed in three locations from the LI 900 tiles. The erosion 
formed in the regions where the coating had been repaired prior to testing. Primary 
cause of the problem was that good bonding of the repair coating was prevented due 
to the silicone waterproofing on the tile. Also, the repaired coating could not be 
thermally fused in-place prior to testing. 

The leading edge of the transverse gaps in the panel also experienced progres- 
sive erosion during the tests. The erosion was caused by a fine alumina dust 
originating at the combustor liner. Exposure duration per run was reduced to four 
seconds to minimize this erosion problem. 

Finally, heating of the center, thin skin, metallic tile resulted in surface 
warping. Surface contour maps of the metallic tile are shown in Figures 21 and 
22 following test runs 4 and 11 respectively. 

Analysis of the abovf: data is documented in Section 4.6 and correlation in 
Section 5. Information relating to the test and data are documented in Volume II 
of this report. 

3.7 RSI Tile Test? in the Ames 20 MW Turbulent Duct - Gap heating tests were 
conducted by F. J. Centolanzi (NASA-Ames Research Center) in the Ames 20 MW Turb- 
ulent Flow Duct Facility to provide heating data in the presence of a high enthalpy 
turbulent boundary layer. The tests employed panels of Silica RSI tiles fabricated 
and instrumented by NASA Ames. The panels were placed in the wall of the Ames 2x9 
inch turbulent duct which employs an arc heated air stream to produce turbulent 
flow. A schematic of the 2x9 inch duct facility is shown in Figure 23 . Panel 
sizes of 20.4 x 25.4 cm and 20.4 x 50.8 cm can be tested at a Mach number of 3.5. 
Figure 24 shows e photograph of the facility. 
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The test program Investigated heating effects on a transverse butt Joint de- 
sign. Gap widths of 0, 0.127, 0.180, 0.254, 0.381 and 0.508 cm were evaluated for 
a RSI thickness of 5.08 cm and an edge radius of 0.635 cm. The test section 
Reynolds number was 0.3 x 10^ per meter. The test specimen Is shown In Figure A-25. 
As shown, the transverse gap was Instrumented on both the upstream and downstream 
sides with thermocouples. The panel was tested twice at each gap width such that 
thermocouple number 1 (Figure 25 ) was located upstream and downstream of the 
Instrumented gap. 

Temperature response data were assimilated and analyses were performed In- 
cluding temperature history comparison and Inverse solutions. The analyses are 
reported In Section 4.7. 
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A.O DATA ANALYSIS 


A.l Analysis of Data Assimilated During Phase I - During Phase I efforts, 
heat transfer data were assimilated from tests conducted at the NASA JSC 10 MW Arc 

Tunnel, the LaRC Mach 10 Continuous Flow Hypersonic Tunnel, the LaRC Mach 8 Vari- 
able Density Tunnel, and the Ames 3.5 Foot Hypersonic Wind Tunnel. A detailed dis- 
cussion of the analysis of these data is given in Section 4.0 of Reference 
1 . However, significant results and conclusions are included in this report for 
reader convenience and report completeness. 

4.1.1 Heat Protection Ability of Candidate Joints - The heat protection 

performance of candidate RSI joint configurations was compared based on 

maximum bondline heat-up (temperature rise) rates. The temperature responses were 

measured during tests in the NASA JSC 10 MW channel nozzle arc tunnel. Four joint 

configurations were tested (butt, contoured, inclined and overlap block). The butt 

joint was tested with forward- and aft-facing steps at the transverse joint and 

with gap wall emittances of 0.6 (white) and 0.9 (black). The c'her configurations 

were tested only with "white" walls. The term "white" walls refers to tiles having 

white gap walls (e “ 0,6) except for the first 0.635 centimeters down the gap, 

which is black (e - 0.9). Gap widths of 0.127, 0.234, 0.381 and 0.762 centimeters 

were tested for each combination of other test variables. Tile thicknesses of 3.18, 

5.08 and 6.35 centimeters were tested, but not for all joint configurations. The 

high cross range shuttle orbiter A2P entry heating rate-time history was simulated 

in the 10 MW channel nozzle for each test run. These test conditions resulted in 

a laminar boundary layer displacement thickness of approximately 1.02 cm, a Mach 

number of approximately 4.5, and a theoretical cold wall flat plate heating rate 

2 

of up to 27.23 watts/cm . The conclusions drawn from the joint configuration 
comparisons are summarized below. 

For the candidate joints, the rates of bondline heat-up in the transverse gap 
are shown in Figure 26 as functions of gap width. As seen, the contoured joint 
affords the best heat protection. At large gap widths the variation in heat protection 
ability among joint types is substantial with the forward-facing step model experiencing 
the largest temperature rise at the 0.762 cm gap width. This is in contrast to the 
aft-facing step which affords almost as much heat protec' ion as the contoured 
joint. Thermal response for transverse gaps indicates that for small gaps (less 
than 0.381 cm), the forward-facing side of the gap experiences higher bondline 
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temperatures than the shielded aft-facing side of the gap. However, for the widest 
gap (0.762 cm), bondline heat-up rate Is the same for both sides of the transverse 
gap- 

The data generated with butt joint models having tile thicknesses of 3.18 and 
6.35 centimeters and a gap wall emlttance of 0.6 (white) are presented In Figures 
27 and 28 . As was expected, these data show a sharp Increase in heat-up rate as 
the gap was opened. There was no clear differentiation of heating In the transverse 
gap as opposed to tne axial gap. The 6. 35 centimeter thick tiles (Figure 28) pro- 
vide more heat protection due not only to Increased Insulation, but also to a re- 
duced sensitivity to gap width for most locations. Figure 28 shows that the down- 
stream parallel gap location and both of the transverse locations were Insensitive 
to the presence of the gaps for widths of 0.381 centimeters or less. Bondline heat- 
up rate (6.35 cm tile) at the upstream parallel gap location shows approximately t’ae 
same sensitivity to Increased gap width as do the parallel gap measurements for the 
3.18 and 5.08 centimeter thick tiles. 

Similarly, a comparison was made of data generated with butt joint models 
having tile thicknesses of 3.18, 5.08 and 6.35 centimeters and a gap wall emlttance 
of 0.9. The joint bondline heat-up rate also shows a strong- sensitivity to gap 

m 

width. Neither the axial or transverse gap orientation appears to be consistently 
hotter. 

A direct comparison of the "black" and "white" coatings indicated the bondline 
heat-up rates for a given tile thickness are similar for each coating. In general, 
the data showed that bondline heat-up rates were equal or slightly lower for the 
high emlttance (e - 0.9) walls. Consequently, It was concluded that Increased gap 
wall emlttance has little effect on bondline heat-up rate. 

Either forward or aft-facing steps at the transverse gaps can be caused by 
manufacturing tolerances, structure deflection, etc. Both these configurations were 
tested by using two upstream butt-gap tiles of 5.08 centimeter thickness and two 
downstream tiles 5.46 centimeters thick to create forward-facing steps, or two down- 
stream tiles 4.70 centimeters thick for aft-facing steps. The results of these tests 
are illustrated In Figure 29 and compared with results >f the tests of 5.08 centi- 
meter thick butt-joint tiles with flush surfaces. For small gap widths, both forward 
and aft facing steps produced bondline heat-up rates lower than those of the flush 
tiles. The aft-facing steps provide better thermal protection due to the shielding 
effect while the forward-facing step is cooler because of the thicker tiles (5.46 cm 
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as opposed Co 5.08 cm). As the gap Is opened, the aft-facing step maintains lower 
response rates, but the forward- facing step experiences bondline haat-up rates which 
Increase to levels well above those of the flush tiles. 

Alternate gap configurations also offer an opportunity for reduction of gap 
heating. Three configurations, the contoured. Inclined and overlap block In addi- 
tion to the butt Joint were tested. In each case the hot external flow Is denied a 
direct path to the gap bottom. None of these configurations was an unqualified 
success, yet all achieved a reduction In bondline response at some Instrumented 
location. Test data for 5.08 centimeter thick tiles of each candidate Joint con- 
figuration are compared next and related to the 5.08 centimeter butt Joint. 

The measurements taken with the Inclined Joint tile set (Figure 30 ) exhibit 
a particularly wide range of sensitivity to gap width demanding on Instrument loca- 
tion. All locations are essentially equivalent to the butt Joint performance 
(Figure 28 ) at a gap width of 0.127 centimeter. As the gap Is opened, however, 
the heat-up rate of the upstream axial location quickly Increases to a level well 
In excess of the heat-up rates experienced during the butt Joint tests. The data 
taken at the bondline of the downstream-facing transverse gap wall remain at re- 
latively low levels for all gap widths and data taken at the other two locations 
are roughly equivalent to those of the butt Joint configurations. 

The contoured Joint configuration (which was more complicated) produced the 
least sensitivity to gap width of any of the configurations tested (Figure 31 ) . 

It is also the only configuration for which bondline heat-up rate varies signifi- 
cantly with location for the 0.127 centimeter gap width. This configuration did, 
in fact, experience higher bondline heat-up rates In the downstream axial ^ap at 
widths of 0.127 and 0.254 centimeter than does the butt Joint. With the largest 
gap v'idrh, though, the contoured Joint provided significantly improved heat pro- 
tection at all locations, compared to the butt Joint. 

The overlap block configuration creates a tortuous path for gas circulating 
from the surface to the bondline. Figure 32 , however, shows that the RSI filler 
block which Is used to create that devious path suffers quite a high bondline heat- 
up rate. By comparison one may see that the other configurations hold no advantage 
over the slmplf^ butt Joint. 

The results of these tests Indicate then, that If small gap widths can be 
achieved, little can be gained by use of Joint configurations more complex chan the 
butt Joint. If gap widths approach the local boundary layer displacement thickness 
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(about 1.016 centimeter for these tests), use of the contoured joint configuration 
may afford considerable relief and forward-facing steps at the joint may exact a 
considerable penalty. 


A. 1.2 Heating Rates in RSI Models of Gaps - Heating rates were calculated 
using RSI-gap thermal response histories measured during tests in the JSC 
10 MW Arc Tunnel channel nozzle. The test conditions, model description and data 
assimilated are dl& cussed in Section 3.1 of this report and in Volume II of Refer- 
ence 1 . The analyses concentrated on a description of radiant heat exchange 
within the gap, a graphical description of the various modes of heat transfer in 
the gap, evaluation of data uncertainties, and comparison of convective heating 
distributions for transverse and in-line gaps for both butt and inclined joint con- 
figurations . 

Radian c Heat Exchange - Radiation exchange between the faces of a gap is im- 
portant because relatively high temperatures (816°C or more) are experienced at 

gap depths of more than 1.27 centimeter. Preliminary data indicated that convec- 

2 

tive heating at such depths was on the order of 0.113 watts/cm . For example a 

difference between two infinite plates (each with a 0.6 emittance) at 816“C 

2 

produces about 0.113 watts/cm net heat transfer. The 8.3*C represents only a 1% 
change in gap wall temperature. Initial thermal modeling employed a minimal number (8) 
of nodes within the gap. Consequently the radiation exchange modeling was refined by 
increasing the number of nodes on the gap wall from 8 to 18. Smaller area nodes are 
particularly helpful because they allow view factors and the nodal tenqjeratures asso- 
ciated with them to approach the ideal condition of an infinitesimal element model. 

In modeling a 6.33 centimeter thick specimen, node lengths were decreased from 0.794 
to 0.353 centimeters, alleviating the situation where large nodes are used in an analysis 
of a small gap. In the ;;ase of large nodes, the only significant view factor may be 
with the opposite node, virtually eliminating the opportunity for emitted or reflected 
energy to be transferred down the gap. Node size in the thermal model was varied so 
that the smallest nodes occurred in the region of highest temperature. 

Graphical Description of Heat Transfer Modes - The heat transfer for the wall 
of a typical forward facing, transverse gap is segregated into its three components 
of convection, radiation and conduction in Figure 33 . The analyses used an 18 node 
model down the gap and the inverse solution method. Temperatures for the uninstru- 
mented wall of the gap were set equal to the values of the corresponding nodes on 
the instrumented wall. Heat fluxes less than zero indicate energy leaving the sur- 
face at that depth. The conductive flux consists of two parts, the conduction be- 
tween adjacent coating nodes and conduction into the RSI, normal to the surface of 
the gap wall. 
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The convective heating and radiative flux shown in Figure 33 are large at the 
top of the gap and their distributions approximately mirror one another. At appro- 
ximately 0.508 centimeters into the gap, the convective and radiative heat fluxes 
become minimal. The net conductive flux has the greatest magnitude near the top 
of the gap, as expected, but much smaller than either the radiative or convective 
components. The negative net conductive flux at the top of the gap is indicative 
of a temperature distribution which is approaching steady state and the remainder 
of the conduction curve is most likely due to the transient thermal response of 
the RSI. During peak heating, the gap wall is hotter at all locations than later- 
ally adjacent RSI material so the conduction in the normal direction is always 
negative. 

Impact of Uncertainties - Studies were conducted to determine the impact of 
uncertainties in gap heating rate distribution on Shuttle TPS performance. This is 
of importance since bondline temperature is a factor which must be controlled through 
proper TPS design, and the adhesives and strain isolation sponges available for this 
type of system have relatively low temperature capability compared with the RSI to 
which they are applied. The results of this analysis give an indication of the level 
at which gap heating ceases to be a significant factor in determining bondline tem- 
perature. Hypothetical cut-off points were assumed in the gap above which convective 
heating rates are known and below which they are uncertain. An extremely conserva- 
tive design approach would then be to assume the heating value remains constant from 
that point. 1'fhile a nonconservative approach would be to assume no heating below the 
cut-off point. The true heating rate distribution of course lies between these two 
assumptions, tor the test case studied and for gap depths greater than 2.54 cm, the 
two assumptions gave comparative results. Thus, if the heating rate distribution is 
known accurately down to a depth of 2.54 cm, the ill effects of uncertainties below 
that gap depth are minimal. 

Comparison of Heating Rate Distributions in RSI Models of Gaps - The convective 
heating analyses of gap tests contained in this section were performed on data ob- 
tained in the NASA-JSC 10 MW Arc Tunnel facility utilizing series of adjustable RSI 
models installed in one wall of a channel nozzle. Convective heating in the gaps 
between tiles was calculated using the MDC General Heat Transfer Program inverse 
solution technique. A description and the method of utilizing this technique are 
given in 3^ction 4.2.1 of Reference 1 . 

Convective heating results for the transverse and in-line gaps were obtained for 

the butt and inclined joint configuration. Four gap widths (0.127, 0.254, 0.381, and 
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0.762 centimeter) and three tile thicknesses (3.18, 5.C8, and 6.35 cm) were analyzed 
for the butt joint configuration. The inclined joint was analyzed for the 6.35 cm 
tile and all gap widths. 

Heating distributions for the downstream wall of the transverse gap for the 
butt joint configuration are shown in Figures 34 , 35 , and 36 . The data on 

Figures 34 and 35 are presented In rectilinear and semi-log coordinates to em- 
phasize the low magnitude of convective heating in the transverse gap at depths 
beyond 0.762 centimeter for the arc heater environment produced by the 10 MW facil- 
ity. In comparing these figures it is seen that heating drops off rapidly down the 
gap, heating Increases with gap width and penetrates deeper into wide gaps, and 
that for a wide gap (0.762 cm) the heating rate ratio near the top of the gap can 
be higher than 1.0. Figure 37 is another way of presenting the above data, as 
a function of gap width, and indicates that for most conditions increasing gap depth 
lowers the gap heating distribution. 

Similarly, in-line gap heating distributions were calculated for the butt joint 
configuration and gap depths of 3.18, 5.08, and 6.35 cm. The results are 
summarized in Figure 38 . A comparison of the data for the two type gaps (in-line 
versus transverse) indicates that heating in the in-line gaps is higher than for the 
transverse gaps at some combinations of gap width and depth. (See Figures 38 and 
37 .) A comparison of these figures also indicates the in-line gap is more sensi- 
tive to gap width than the transverse gap for gap depths greater than 0.2 cm. 

The effect of a forward facing step on heating distributions is shown in Figure 
39 for a 5.08 cm tile and the whole range of gap widths. For the narrow gaps 
(0.127, 0.254, 0.381 cm) the heating near the top of the gap stagnates, increasing 
the heating rates. Since the gap is small, flow within the gap is f^'^eded, 
therefore the heating falls off rapidly to 0.9 cm gap depth and then de- 
creases gradually below that point. For the wide gap (0.762 cm) the heating takes 
on a different distribution which is also significantly higher than for the narrow 
gaps. Since the gap is wide the heating near the top gets relief from the wide gap 
below, which causes heating to recirculate and penetrate deep into the pap Because 
the step distorts the gap heating distribution, a significant increase in heating 
occurs whan compared to a flush joint. Figure 40 shows a comparison of downstream 
transverse gap heating for a butt joint model with and without a forward facing step. 
As seen, heating along the gap wall of a forward facing step can be an order of magni- 
tude or more greater than for flush tiles. 
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Figure 36 
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Figure 37 
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Figure 38 
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Convective heating analyses for an Inclined joint model were also performed 
for both the downstream side of the transverse gap an * for an in-line gap upstream 
of the transverse gap. The tiles were 6.35 centimeters thick. The heating distri- 
butions were found to be very similar to the butt joint distributions. For both 
in-line and transverse gaps, the heating for the inclined joint was slightly greater 
and more sensitive to gap width from O.IC to 1.0 cm gap depths. 

4.1.3 Analysis of Mach 10 CFHT Test of Gap Model - Analy.-. es were perfomu-d on 
heat transfer dat obtained on a wall-mounted, thin skin tile model tested in the 
LaRC Mach 10 Continuous Flow Hypersonic Tunnel. Test conditions, model description 
and data assimilated are discussed in Section 3.1 of this report and in Volume II 
of Reference 1 . The test article consisted of a panel with six RSI tiles sur- 
roundin'^ a highly Instrumented thin skin tile. The test panel was located on a 
rotational plate on the tunnel sidewall such that the flow angle could be varied 
over the test panel. The tests performed during the initial program were condotted 
at a Mach number of 10 and a unit Reynolds number of 3.3 x 10^. Included in the 
data analyses are the following: 

a) fcivaluatlon of data reduction metnods including calibration plate heating 
and the effect of considering conduction on measured heating rates. 

b) In-line versus staggered tile heating patterns. 

c) Effect of gap width on tile heating patterns. 

d) Effect of flow angle on gap heating patterns. 

e) Effect of steps on tile heating patterns. 

Data Reduction Methods - LaRC performed calibration runs in support of a 
McDoni'.ell Douglas sponsored program to measure heat transfer data on a corrugated 
panel model mounted on the tunnel sidewall. As part of that effort, flat plate 
heat transfer data were tsken. Figure 41 presents the measured heat transfer 
distrlbutioT, in the vertical direction on the flat date mounted on the tunnel 
sidewall. The distribution shown is based upon three data runs. A significant 
variation in the heating across the flat plate is observed. This tunnel character- 
istic has been attributed to the square nozzle and test section which results in a 
slight flow convergence toward the center of the tunnel sidewall. Data taken on 
the corrugated panel exhibited a similar spanwise heating gradient to that observed 
on the flat plate. Normalizing the corrugated panel heat transfer coefficients by 
the flat plate coefficients 'esulted in successful collapsing of the data in the 
spanwise direction. Because of this previous experience in correlating the corru- 
gated panel data, the gap heating data, taken at the same test condition were 
normalized by the measured flat plate heat transfer coefficients. 
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Data reduction for the thin skin tile tests employed the slope of the 
temperature-tine curve at selected times to solve for the heat transfer coefficient 
using the following equation: 


pCp X (dTw/de) 

h * — (Taw ' -fw) where X ■ skin thickness 

Two times were selected for data reduction for each thermocouple during these tests . 
The first time selected was 0.50 seconds after the test article reached the tunnel 
wall since during the first 0.35 to 0.45 seconds the temperature readings were 
erratic. A second time was also selected which was 0.50 seconds after the first 
time or 1.0 seconds after tile insertion. The temperature-time derivative (dT/d0) 
was obtained by taking the slope of a least squares quadratic curve fit through ten 
seconds of data obtained for each thermocouple. The initial point of the curve fit 
interval is the time selected for data reduction; i.e., 0.50 seconds and 1.0 second. 
Ten seconds was selected as the curve fit interval to obtain data deep in the gaps 
where heating levels are low. The curve fit expressions are of the form: 

2 

Tw=a+b0+c6 and dT/d0 = b + 2c 0 


where a, b, and c are constants. Two heat transfer coefficients were computed fcr 
each value of dT/d0 based upon two values of adiabatic wall temperature, i.e., 

Taw/T^ = 0.895 and 1.0. However, all data in this section ^re based upon Taw/T^ = 
0.895 because the boundary layer was turbulent for all tests. Also all d^ta pre- 
sented in the following figures were evaluated at 0.50 seconds after test .irt ;le 
insertion was complete. 

The above method is graphically demonstrated in Figure 42 for two typically 
measured temperature histories on the downstream face of the thin skin tile . In 
the first curve 'TiC Channel 20), the thermocouple was located in the gap 0,414 cm 
from the surface, and the two slopes are nearly equal indicating a linear tempera- 
ture response. All thermocouples where a significant temperature rise occurred 
exhibited this type of temperature response. In the second curve (T/C Channel 55), 
the thermocouple is deeper within the gap (1.113 cm from the surface), and 
the temperature response was an order of magnitude lower than Channel 20. The 
temperature response was also low enough that "noise" from the data recording system 
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is observed in the data. The third uhrough the sixth rows of thermocouples had a 
temperature response so low that the "noise" from the data recording system was 
greater than the temperature rise over a ten secona time interval. Consequently, 
a significant number of thermocouples In the gaps had such a low signal to noise 
ratio that the data were not Included In the analysis. 

A study was made to examine the effect of including conduction In the calcu- 
lation of heating rates from wind tunnel tests that employ thin skin tiles. Data 
from Run 14 of the CFHT were selected for this purpose. An eight node thermal 
model was formulated to describe the heat storage and heat conduction characteris- 
tics of a section through the thin skin tile for the temperature distributions 
shown in Figure 43 . For each node in the thermal model there was a correspond- 
ing thermocouple on the thin skin tile which was used to define the nodes' tem- 
perature history. Hand fairings of the temperature histories were input into the 
General Heat Transfer Computer Program along with the thermal model descriptors and 
an in'^erse solution was performed to calculate a heating rate for each node. Fig- 
ure 43 depicts the spanwlse temperature distribution across the top and down the 
side of the tile at 0.5 seconds after test article insertion. Two distributions 
are shown corresponding to "as received data" and "revised temperature data" . The 
revised data resulted from applying a thermocouple calibration correction to 
the as-received-data. The correction was determined by NASA LaRC 
after initial efforts to examine the effects of conduction were unsuccessful (de- 
tailed discussion is given in Section 4.3 of Reference 1 ). The revised tempera- 
ture distribution is seen to be much smoother than ti tglnal distribution with 
the resultant elimination of the "knee" in the curve t.. existed at node 4. 

Figure 44 shows the effect on calculated heating rates of excluding and in- 
cluding conduction in the thermal model using the revised temperature distributions. 
The data ^hown in this figure is for 0.5 second after tile Insertion. The heating 
rates for the "no conduction" case consider only the heat storage term and are com- 
parable with the data reported by LaRC. The differences between these data are due 
to the techniques used in the fairing of the temperature histories. The technique 
used by LaRC consiited of least squares curve fit of the temperature histories while 
hand-faired histories were used with the eight node thermal model. When conduction 
is included in the thermal model analysis, the surface node near the top edge of the 
tile ihowed an increase in calculated heating rate of 8.5%. However, the calculated 
heating rates in the gap are generally lower when conduction is included. Conse- 
quently, it was concluded that the heating rates neglecting conduction (i.e., data 
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RSI GAP HEATING ANALYSIS - II 

VOLUME I 

obtained from LaRC) be used In the gap heating data correlation since these data 
give conservative results. 

Heating Patterns. In-line Versus Staggered Tiles - Data were taken on the thin 
skin tile with the two basic tile arrangements of staggered and In-line. Figure 
45 shows the orientation of the test article with respect to the flow for both 
arrangements. The In-line arrangement Is achieved by rotating the test article 90 
degrees from the staggered tile orientation. It should be noted that for the In- 
line arrangement, the tiles are In-line In the axial direction only. In the span- 
wise direction (normal to the flow) the tiles are staggered. Comparisons of the 
axial heating distributions for the staggered and In-line tile arrangements at 
y = 0.0, -3.8 and -7.3 cm are presented In Figures 45 , 46 , and 47 respectively. 

These figures present data for a tile gap width of 0.23 cm. Heating on both the 
upstream and downstream faces of the tile does not appear to be significantly af- 
fected by the surrounding tile arrangement for a gap width of 0.23 cm. Heating on 
the top surface of the thiu skin tile is higher (4% to 24%) for the staggered 
arrangement than for the In-line arrangement. This trend is most pronounced at 
the centerline (y = 0.0 cm) of the tile and decreases near the edge (y = -7.3 cm) 
of the tile. Also the magnitude of the heating on the top of both tile configura- 
tions decreases near the tile edges. 

Effect of Gap Width - Comparisons of heating distributions for four gap widths 
at y = 0.0, -3.8 and -7.3 cm are presented for the staggered tile arrangement in 
Figures 48 , 49 , and 50 , respectively. Data for gap widths of 0.13, 0.23, 

0.46, and 0.71 cm are shown in each figure. These figures show that the effect of 
gap width on tile heating changes with location on the tile. Figure 48 presents 
data along the centerline of the tile (y = 0.0). On the upstream face, the gap 
heating increases slightly with increasing gap width. This trend is reversed on 
the top surface and downstream face with the exception of the upstream edge of the 
tile top surface. Examination of the axial distribution at y = -3.8 cm shows a 
much greater increase in gap heating with increasing gap width on the upstream face 
than was shown at y = 0.0. On the top surface of the tile the heating at the up- 
stream edge of the tile increased dramatically with increasing gap width. Over the 
rest of the top surface, the heating appears essentially independent of gap width. 

On the downstream face the trend is mixed with the gap heating either increasing or 
decreasing with gap width depending on the depth into the gap. The data near the 
outer edge of the tile (y ■ -7.3 cm) show the heating on the upstream face, top 
surface, and downstream face all increasing significantly with increasing gap width. 
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In general for the three (y) stations examined, It appears that the gap heating 
is related to the heating on the top surface near the gaps. The same trends that 
apply to the gap heating also apply to the top surface heating near the tsdges of the 
tile. This can be seen in nlJ the figures presented for thp staggered tile arrange- 
ment . 


1 

I 

I 



Similar comparisons of heating <^i«t.rlbutlon8 for the in-line tile arrangement 
were made and arc repor*-cJ in Reference 1 (Section 4.3). 

Review of gap width effects on both staggered a.'.: in-line tile arrangements 
indicates chat the in-line tile arrangement results in lower and more uniform heat- 
ing on the top surface of the tile. The gap heating on the upstream faces of the 
tile do not appear to be significantly different for the two tile arrangements and 
a similar conclusion can be drawn for the downstream face. As expected, gap width 
significantly affects gap heating; however, the magnitude of the effect is dependent 
on the location in the gap. 

Ef f ect of Flow Orientation - The effect of flow angle on gap heating was 
examined on two faces ot 'In thin skin tile for 0.229 and 0.710 cm gaps. For ease 
of data handling, a consistent coordinate system was defined which is fixed in the 
thin skin tile. A fl.ov angle (y) was defined which varies from 0 radiant at the 
in-line cor'/igiTatlor. to r/? the staggered configuration. The analyses indicated 
that either an ir.-.,inf 'lie arrangement (y“ 0) or a staggered tile arrangement (y** 
tt/2) Is more desirable than other flow orie. cations . The spread in gap heating on 
the tile faces is minimized at y«0 and .'/2 radians, and the peak heatl g in the 
gaps are minima at these flow angles. These concluslv..'9 a' ^ i-asad on data for 
both gap widths (0.229 and 0.710 cm). Figure 51 shows the lci..:;i;icns and levels 
of maximum gap heating at various flow angles for the f . 710 cm jap. The data were 

measured approximately 0.3 cm from the tile surface, As seen, [ease 

gaps are minima at 0 and r/2 radians flow angles. 

The effect of gap width on heating rates at various fl w orlantatioiii u«.s al'o 
studied. The heating distributions along tile gaps at a flow angle of 't/4 radians 
is illustrated in Figure 52 . The data wes meas >re.. nr g." depth of 0.3 cm, and 
at gap widths of 0.23 cm and 0.71 cm. The heating rr'®- gr.^d Lents, <’ at occur at 
this flow angle become mere severe as the gap see-, in this figure. 

Effect of Steps - The effect of tile mismatch waa exao red Ly ■ ■ shims 

to raise and lower the thin skin tile. Tests were tun with th#' tiu ri > 3 h.~j 

cm above the surrounding ESI cilcs and lowered 0.168 cm below the surrounding 
tiles. Figure 53 presents in-line tile heating distributions (Y^ ■■ 0.0) 


77 

MCOOMMEt-L t>kHJCt MS MSTtrO^MUTtCS • MMST 



] 





1 


... 



RSI GAP HEATING ANALYSIS - II 



VOLUME I 


, 1 ... ..... 




1 

' 4 - 


REPORT MOC E1248 
JSC 09651 


LOCATIONS OF MAXIMUM GAP HEATING 



C CFHT DATA 
• M = 10 

oo 

ft GAP WIDTH = 0.71 CM 

ft T/C LOCATED 

APPROXIMATELY 0.3 CM 
FROM TILE SURFACE 

ft (X.XX) = HT/Hpp 


78 


Figure 51 


«9COOWWr*.L IMl^OLiftS AaT0K>MAUTICa COn^RANV • BMar 








c/> 



1 : _GAP (TOP) : : CORNER I ; 1 GAP (BO 

tn. :.il ; r YT 

4 8 12 16 20 24 

S -V DISTANCE FROM LEADING EDGE OF TILE (CM) 


.. GAP (BOTTOM) : 



Figure 52 


MCOO/V/VEl.f. OOUGLAS MtSTtfOM/^UTICS COMt*AM'V - EAST 






RSI GAP HEATING ANALYSIS - II 

VOLUME I 


REPORT MDC E1248 
JSC 09651 


IN-LINE TILE HEATING DISTRIBUTIONS 
IN CFHT (Yc - 0.0) 


FLOW 


IN-LINE 


B 

L_ 


a 


Jj£ 


a 



GAP WIDTH = 0.23 cm 

GAMMA » 0 RADIANS 

M = 10 
00 

Re /m = 3.3 X 10® 


nr hn ri in n — r~i 

FLUSH STEP-UP STEP-DOWN 


VIEW B-B 



X - DISTANCE FROM CENTER OF TILE - cm 


80 


Figure 53 


‘^COOMMCLl. OOUGLAO ASTHONAUTICS COM9>A^V • mA»T 





a 


RSI GAP HEATING ANALYSIS - II 

a 


VOLUME I 


REPORT MOC EI248 
JSC 09651 


for the step-up, flush and step-down tile configurations with a gap width of 0.23 
centimeters. The step-up configuration resulted in a 55S Ir crease in the peak heat- 
ing rate measured on the tile surface when compared to the flush tile. The peak 
heating location appears to move nearer the upstream edge of the tile as the tile 
is raised above the flush position. Raising the tile increased heating on the up- 
stream face of the gap and entire tile surface while having little effect on the 
downstream face. The step-down tile configuration resulted in lower heating on 
the upstream half of the thin skin tile and similar heating on the downstream half 
when compared with the flush tile. The peak heating rate was only 4.5% lower for 
the step-down tile. Figure 54 presents heating distributions for the three tile 
positions when the surrounding tiles are in the staggered configuration (y = it/2) . 
Similar conclusions can be drawn from the staggered tile data as were found for 
the in-line tile data. 
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4.1.4 Analysis of Mach 8 Variable Density Tunnel Tests - Heat transfer 
measurements were performed on a thin skin tile panel In the LaRc Mach 8 Variable 
Density Tunnel to provide data In the presence of laminar and turbulent boundary 
layers. Model test position included both on the centerline (free stream) and on 
the tunnel wall (flush) . Test conditions, model description and data are discussed 
in Section 3.1 of this report and in Volune II of Reference 1 . Test output was 
recorded on tape to facilitate analysis and for use in the Multiple Regression 
Analysis computer program. The analysis iicluded' evaluation of: 

a) a reference for data correlation 

b) heating patterns on in-line tiles 

c) he? ing patterns on staggered tiles 

d) effect of Reynolds number on heating rates 

e) effect of gap width on tile heating patterns 

Reference for Correlation - Heat transfer rates on the top of the test panel 
were examined to e tablish either a laminar or turbulent flow reference for data 
correlation. Initially, a laminar recovery factor was used to compute a heat 
transfer coefficient (H^^) as well as laminar boundary layer theory to compute a 
referenced transfer coefficient which was bcsed or. distance from the model 

leading edge. The ratio of H^ to indicates tie type of flow with a value of 

unity indicating laminar flow. Analysis indicated that for the free stream tests, 
conditions were laminar on the forward portion of the model at the low Reynolds 
number, and flow transition occurred on the aft portion of the panel. However, for 
the tunnel wall tests, Reynolds number and boundary layer thicknesses were character- 
istic of a turbulent boundary layer. For both the fr»e Ktream and the flush wall 
tests, the referenced coefficient was that measured 2.582 cm from the forward tile 
leading edge. 

Heating Patterns on In-Line Tiles - The measured heating distributions along 
the length of the panel for the in-line tile configuration are shown in Figure 55 
for the free stream tests. Data are presented for unit Reynolds number per meter 
ranging from 2.3 x 10^ to 21.8 x 10^. A variety of heating patterns is exh.oited 
on the cop of the tiles i.idicating laminar flow on tile #1 for the lower Reynolds 
number and transitional flow for the higher Reynolds numbers. For tile #2, the two 
lower Reynolds number tests exhibit a trend toward transitional flow, while for 
Re^/m = 21.8 x 10^ the data suggests that the flow is fully turbulent. The gap at 
the center of the panel did not affect the heating on the top of the panel except 
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in the Re^/m * 6.1 x 10 test where the gap produced transition onset. Transverse 
gap heating is relatively low compared to that measured on top of the panel and 
does not show a strong dependence on Reynolds number for the forward facing wall. 

The aft facing walls experienced a slight change in heating distribution with Reynolds 
number. The larger gaps (0.318 cm) experience higher heating than the smaller gap 
(0.159 cm). 

Heating distributions for the in-line gap model positioned flush with the tunnel 
; wall are shown in Figure 56 . Tests were performed at unit Reynolds number/m from 

1.16 X 10 to 41.4 X 10^. Heating on top of the tiles is relatively uniform with 
the higher Reynolds n'lmber data being almost constant. The lower Reynolds number 
data show a 10% heating increase on the tile top. Heating in the transverse gaps 
' does not decrease as rapidly with distance into the gap as for the free stream test. 

For both the free stream and tunnel wall position, the upstream side of the transverse 
gap experiences equal or higher heating than the downstream side of the gap. Again 
i, the heating is higher for the larger gap. 

Heating Patterns on Staggered Tiles - Heating distributions were also measured 
for a staggered tile ' onflguratlon with the model in the free stream position and 
flush with the tunnel w\.il. Figure 57 contains the heating data for the free 
stream tests. Heating distributions on the top surface of the tile.s are similar to 
those obtained for the in-line gap model (Figure 55 ) over the range of Reynolds 
numbers investigated. The staggered tile configuration had twice the gap width of 
that tested in the in-line model, and only the downstream face of the center trans- 
verse gap was instrumented. Gap heating at the center of the model (stagnation 
region) is also presented in Figure 57 and shows a strong dependence on Reynolds 
number. For Reynolds numbers greater than 11.2 x 10 per meter, heating in the gap 
is greater than on the top surface. 

Heating data for the staggered tile configuration flush with the tunnel wall 
(Figure 58 ) show that the turbulent boundary layer produced essentially uniform 
heating or the top of the panel similar to that experienced by the in-line 
tile model (Figure 56 ) . As with the free stream tests of the staggered tiles, 
the downstream face of the transverse gap experienced significant heating with a 
distribution strongly dependent on Reynolds number. 

Lateral heating distributions measured across the gap face of the downstream 
tile for the staggered tile configuration indicated that heating is almost constant 
across the half gap width and then decreases in the lateral direction similar to a 
"normal" distribution. 
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4.1.5 Analyses of Ames 1.5 Foot HWT Tests - Analyses were performed on data 
obtained from the Rockwell International gap heating tests conducted In NASA-Ames 
3,5 Foot Hypersonic Wind Tunnel. Data were obtained on thin skin test articles sub- 
jected to laminar, transitional, and turbulent flow conditions. Test conditions, 
model description and data are discussed In Section 3.1 of this report and In 
Volume II of Reference 1 . All tests used a flat plate model at zero angle-of- 
attack. Analyses Included evaluation of: 

a) calibration plate heating pattern 

b) heating d: strlbutlons In transverse gaps 

c) heating distributions for In-line gaps 

d) effect pf Reynolds number on heating 

e) eff' ct of gap width on heating patterns 

f) comparison of data trends to CFHT data 

Calibration runs were made to characterize the flow over the test configuration 
utilizing a srrcth heat transfer calibration plate. Three rows (located laterally 
at y = 0.0 and + 20.54 cm) of the rmc couples and three corresponding total temper- 
ature probes were used to measure the heating environment along and across the test 
article. Figure 59 presents the heat transfer coefficient along the calibration 
plate (y “ -20.54) for the four test unit Reynolds numbers. The heat transfer 
coefficient (H ) is based on a recovery factor of 0.874 The low Reynolds number 
heating data decreased approximately with the square root of distance along the panel 
which is characteristic of a laminar boundary layer. The higher i'.eynolds number 
data show a decrease and then a sharp rise in heating along the panel characteristic 
of transitional flow. 

The heat transfer data measured on the thin skin inserts with simulated RSI 
gaps were referenced to the flat plate calibration data. For each test condition, a 
two-dimensional interpolation in the x and y directions was performed on the cali- 
bration data to determine the reference heat transfer coefficient at the specific 
iocat-ons where gap heating data were measured. 

Gap Heating Distributions - Heating distributions in a single transverse ^;ip 
model are shown in Figure 60 for three Reynolds number conditions. Gap width was 
0.254 cm and gap depth was 2.03 cm. The heat transfer coefficients were normalized 
to the calibration plate coefficients at the same x and y locations. Heating dis- 
tributions are shown on the surface of the model both forward and aft of the gap 
well as both faces in the gap. The heating forward of the gap increases and th r 
decreases with distance for all Reynolds number conditions, and the heating parameter 
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(H^/H/pp) consistently decreases with Increasing Reynolds number. The heating In 
the gap shows consistent Reynolds number trends with the effect of Reynolds number 
much less than Is observed In the surface heating data. The surface heating aft 
of the gaps has Irregular distributions with an apparent Inconsistent Reynolds 
number trend. It can also be noted that the surface heating rates on the single 
t .tnsverse gap model were generally less than that measured on the calibration 
plate, l.e., < 1.0. 

A comparison of surface heating distributions on the downstream side of a 
transverse ga_. Is made In Figure 61 on data obtained from the Ames 3.5 foot HWT 
test and the LaRC CFHT test. The models used In the Ames tests had an ample 
number of thermocouples on the surface to define the heating rate distribution near 
the gaps while the model used In the CFHT tests did not. The Ames data indicate 
that the peak surface heating rate occurs approximately one edge radius dov^nstream 
of the tile leading edge. A dashed curve has been added to the CFHT portion of the 
figure to Illustrate a plausible heating distribution which would he consistent with 
the AMES data. In addition, the AMES data show a level of sensitivity to unit 
Reynolds number, i.e., as the Reynolds number is increased, the heating ratio (H/Hpp) 
Increases. 

Heating distributions for In-line gaps from the Ames 3.5 foot HWT and the CFHT 
were also compared (See Figure 62 ) . These data are for the top of the tile near 
the gap where effects of gap flow should be evident. Both sets of data show an 
increase in (H/H^p) with distance along the gap. Also, the in-line gap data increase 
with Reynolds number as do the other data measured at AMES . 

Heating distributions for the upstream and downstream sides of gaps oriented 
at 30 degrees and 60 degrees to the flow (Figure 63 and 64 ) show similar 
dependence on unit Reynolds number. As the Reynolds number increases, the heating 
distribution intensifies. At 30 degrees, heating on the panel is monotonic with 
distance with a sharp drop near the gap. Heating on the downstream side of the 
gap shows a definite enhancement due to the gap. Reynolds number effects are 
evident on the downstream side of the gap. For the 60 degree orientation, heating 
on both sides of the gap is dependent on Reynolds number. It should be noted that 
the 60 degree data are downstream of the disturbance caused by the 30 degree gap. 

Also shown on both figures are comparable data measured during the CFHT tests. 

At the 30 degree orientation, the heating distributions from the CFHT and the AMES 
3.5 foot tunnel are similar. At the 60 degree orientation the distribution measured 
at AMES Is higher than measured at the CFHT. This Is probably due to the transitional 
flow present in the AMES tests and the fully turbulent conditions in the CFHT. 
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4.2 Analysis of Silica RSI Testa In the JSC Laminar Duct - Analyses were per- 
formed on data from the Rockwell International gap heating tests In the NASA-JSC 
10 MW Arc Tunnel (Reference 2) . The testa eiqployed silica RSI tiles fabricated 
and Instrumented by NASA Ames. The test conditions, model description, and data 
assimilated are discussed In Section 3.2 and Volume II of this report. The purpose 
of the program was to obtain comparative thermal performance data on overlap and 
butt Joint designs, and to evaluate the effect of tile thickness and joint gap width 
under laminar flow conditions . 

Six runs from the test program were selected for analysis. All six runs were 
made with the staggered gap design consisting of three full and two half tiles with 
a perimeter of narrow guard tiles. Gap widths of 0.127 and 0.254 cm were evaluated 
with RSI thicknesses of 2.54 and 5.08 cm. Four of the six runs were made on butt 
joint models and two runs on overlap Joint models. Mass flow rate was 45.4 grams/ 
sec In all cases and enthalpy varied from 8.58 x 10^ to 21.95 x 10^ J/Kg. All runs 
were conducted In the channel nozzle of the JSC 10 MW Arc Jet Facility. The test 
panels were mounted In one wall of the channel nozzle and a calibration plate was 
mounted on the other wall. 

Included In the data analysis are the following: 

a) Selection of data for analysis 

b) Analysis method Including thermal models and thermal properties 

c) Butt joint heating patterns 

d) Overlap joint heating patterns 

e) Comparison of butt and overlap joint heating patterns, and comparison 
to other test data 

f) Sensitivity studies 

4.2.1 Selection of Data for Analysis - Six runs from the test program and one 
thermocouple stack from each of the six runs were selected for analysis. The initial 
set of selected runs and thermocouple stack (T/C) are listed in Figure 65 as well 
as joint type, tile thickness, and gap width. 

In reviewing the data It was found necessary to make certain substitutions as 
described below and Indicated In Figure 66 . In Run 528, T/C stack 271 was analyzed 
instead of 159. By symmetry, conditions should be identical at these locations but 
measurements at locations 272 and 156 were desired In the analysis and no analogous 
data were available for stack 159. 
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BUTT 
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IN-LINE 
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BUTT 
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In runs 538 and 539, the requested T/C stack 164 Is located in the upstream In- 
line gap but data for gap thermocouples 1 and 3 were missing. Since adequate tem- 
perature definition at and near the top of the gap is required for the inverse solu- 
tion no attempt was made to compute heating at location 164. All data appeared good 
at location 169, and this location was used Instead for the analysis of Runs 538 and 
539. Plug data for this tile were largely missing for both of these runs, with only 
T/C No. 7 appearing valid; hence, Indepth temperatures were computed and used in the 
place of the plug data. 

For runs 542 and 544, data at locations 262 were complete and these runs were 
analyzed as requested, with temperature boundaries defined by plug measurements at 
the center of the tile. 

4.2.2 Analysis Method - The MDAC HEATRAN inverse solution technique was used 
to obtain gap convective heating distributions from the test temperature data. Test 
results used in the analysis consisted of gap and plug temperature histories, teuc 
section pressures and channel wall temperature. The calibration plate served as a 
basis for referencing data and as a basis for investigating anomalies. Thermal 
models, material thermal properties and solution/test times are discussed below. 

Thermal Models - Two thermal models were developed for this analysis — one for 
the butt joint configuration and one for the overlap joint with filler bar. The 
butt joint model is shown in Figure 67 and the overlap joint model in Figure 68 . 
Both models are modifications of an existing model described in Reference 1 . 

In the butt joint model the RSI and waterproof coating are divided into 17 
layers vertically and each RSI layer was divided into two sections laterally. All 
temperature node locations are defined in Figure 67 . Node spacing is com i ant be- 
tween nodes 1 and 4, 4 and 15 and from 15 to 18; within these limits the spacing 
may be varied. The coating on cop of the RSI i.s not shown in the sketch but was 
assumed to be the same thickness as in the gap and its effect is accounted for in 
the conduction and storage terms associated with nodes 22 and 43. 

Radiation among the model components was accounted for as indicated by radlosity 
nodes 83 through 124 and 128 through 132. At each time step radiation view factors 
were computed using the crossed-strings method. 

The 0.216 cm layer under the RSI tile consists of a 0.152 cm thick strain iso- 
lation pad with a 0.025 cm RTV bond on top and a 0.038 cm RTV bond on the bottom. 
Effective thermal properties were computed from the actual properties for this com- 
posite layer. Two conductivities were used, one defining resistances in parallel 
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for lateral conductioxi and one for vertical conduction In which the resistances were 
treated as In series. Under this layer Is 0.152 cn rf aluminum, extending to the 
edge of the gap and under this a 0.218 cm sheet of aluminum which functioned as a 
support plate for the entire panel. At the bottom Is 2.54 cm of TG 15000 insulatlo.T. 

Temperature nodes 43 through 60 are located to coincide with the plug thermo- 
couple stack. Temperatures of these nodes were forced to follow the measured tem- 
peratures when the latter were avallal^le, i.e., a temperature boundary was defined 
by the data. When no plug data were available the temperatures of nodes 43 through 
60 were computed. For gaps aligned parallel to the flow, temperature and heating 
conditions shculd be symmetrical on both sides of the gap and were so assumed In the 
analysis . 

Thc! foregoing discussion of the butt Joint model Is generally applicable to 
the overlap Joint model. Figure 68 , except that the latter was expanded to Include 
the filler bar and the waterproof coating on its exposed surface at the bottom of 
the gap. The filler material is FI 600 which is a lightweight fibrous material 
which was compressed at time of Installation and hence is assumed to fill th 
closure alleviating the need to account for air gaps or contact reslstam je 

to the low thermal expansion coefficient for silica RSI, 0.54 x 10 ^ changes 

in gap width during a test were insignificant. Consequently, gap width was held 
constant during the analysis for each test. 

Thermal Properties - Thermal properties were used in the analysis as provided 
by Rockwell International except that those for the RTV bond and strain Isolation 
pad were combined to provide effective properties for a composite layer as dis- 
cussed above In the Thermal Mjdel section. Pressure remained relatively constant 
during test periods, hence pressure dependent properties were defined at a mean 
pressure obtained from measurements during each run along the centerline of the 
calibration plate in the opposite wall of the channel nozzle. Interpolation with 
pressure In the property tables was logarithmic. 

A nominal emlttance of 0.8 was used for the coating fired on the tiles by 
Ames. The channel wall emissivlty was assumed to be 0.77 and that of Che aluminum 
exposed at Che bottom of the butt Joint gap was assumed to be 0.35. 

Inverse Solution Times versus Test Times - Test times for the six runs analyzed 
was from 269 to 2040 seconds, depending on tile thickness, gap size and gap configure 
tlon. The Inverse solution was applied from the start of temperature rise to a time 
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at which teaperatures near the tile top surface approached steady state conditions, 
or to the end of the test time if surface steady state conditions were not achieved. 
Test tlDCS and inverse solution times are tabulated in Figure 69 . 

4.2.3 Bute Joint Heating Patterns ~ Of the four sets of butt Joint data analyzed 
(See Figure 65 ) only that from runs 542 and 544 Included nearly complete plug tem- 
perature distributions. For these luns no data were available for the second thermo- 
couple from the top of the plug (T/C 7) but the remaining data appeared to provide 
adequate definition of the temperature distribution in the plug, hence were used to 
provide a temperature boundary for the inverse solution. Figures 70 and 71 show 
heating rate distributions normalized to the surface value. Figure 70 'employs a 
linear scale and was used to fair curves through the computed points; Figure 71 
presents the same Information In log-log form. 

The same Instrumented 5.08 cm thick tiles were used In these two tunnel runs, 
with the gap width (0.127 and 0.254 cm) being the only parameter changei. The be- 
havior of the solution for Run 542 as the heating rate approaches zero was found to 
be typical of those for the 0.127 cm gap data and considerable effort was devoted to 
understanding it. A curve was faired through the data assuming that the sharp drop 
to a negative heating r'te is spurious. However, it was noted early in the study 
that the 0.127 cm gap data showed temperatures unexpectedly and significantly lower 
than plug temperatures at the same depth. Figure 72 shows the distribution of gap- 
plug temperature differences for runs 542 and 544. For the 0.254 cm gap, the plug 
and gap temperatures are within 14 C deg at the surface, the distribution varies 
smoothly and the gap temperatures are higher than the plug temperatures over most 
of the range. However, the 0.127 cm gap data shows plug temperature significantly 
higher than gap temperatures to a dfpth of 3.43 cm. This difference reaches a sharp 
maxlmim of 99°C at the same location as the spurious — appearing negative heating 
rate in Figure 70 and Indicates an Internal conduction from the plug toward the 
gap. 

The cause of the anomaly has not been determined. The fact that surface tem- 
peratures are higher at tb'' plug than at the gap tor both runs 542 and 544 (••• Figure 
72) suggests a variation in heating conditions acrocs the channel nozzle. Plug 
temperature histories are identical for both runs except for the initial steep rise 
period of about 40 seconds, indicating good run repeatability. The same instrumented 
tiles are used for both runs, so that differences in thermal properties tentatively 
may be ruled out. Gap width is the only known variable and Figure 72 suggests a pos- 
sible variation or nonuniformity in gap width which might affect both the cavity 
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0.127 

600 

330 

533 

5.08 cm OVERLAP 

0.127 

2040 

590 

538 

2.54 cm BUTT 

0.127 

556 

240 

539 

2.54 cm BUTT 

0.254 

269 

220 

542 

5.08 cm BUTT 

0.127 

1724 

1540 

544 

5.08 cm BUTT 

0.254 

562 

460 
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flow and radiation characteristics. A given variation might be a significant frac- 
tion of the smaller gap while having no obvious effect on the larger gap data. 

Results of the analysis on the two 2.54 cm thick butt joint tests, runs 538 
and 539, are shown in Figures 73 and 74 . For these runs, the data were inade- 
quate at the requested location (no. 164) on an upstre 2 un tile and location 169 on 
the center downstream tile was substituted. Since data were available from only 
one RSI plug thermocouple, the temperature at the plug nodes were computed as a 
part of the solution. The surface temperature on the top of the tile was assumed 
uniform although this was found not to be the general case. For the 0.127 cm gap 
(run 538) data from the one functioning thermocouple indicate plug temperatures 
appreciably higher than gap temperatures at the same depth; the same thermocouple 
in the 0.254 cm gap run (run 539) indicates plug temperatures equal to or lower 
than corresponding gap temperatures. The heating distribution in the 2.54 cm tile 
and for the 0.127 cm gap exhibits the same form as that for the 5.08 cm tile, and 
the distribution of the 0.254 cm gap has a small irregularity at the same depth. 

Gap/plug temperature differences are shown in Figure 75 for these runs. 

These distributions have the same form as those in Figure 72 for the 5.08 cm 
tiles. The initial irregularity (within 0.254 cm of the surface) is probably due 
to the ajsumption of constant surface temperatures between the gap and plug. The 
minima at a depth of 0.305 cm again correspond to the erratic points in the heat- 
ing distribution curves in Figure 73 . These distributions both become asymptotic 
to the same slightly negative value, ''' “0.02. Since both of these curves 

were obtained from the same instrumented tiles the effect is presumably associated 
with the ti’.e or instrumentation, i.e., an error associated with a given thermo- 
couple ading or location and/or a deviation in coating thickness. 

».2.4 Overlap Joint Heating Patterns - The heating distribution for run 533 
is shown in Figure 76 and 77 and the gap/plug temperature difference distri- 
bution in Figure 78 . The test used a 5.08 cm tile with overlap in-line joint 
and 0.127 cm gap width. These figures exhibit the same phenomena observed in the 
other 0.12"’ cm gap results, with the plug temperature 130°C higher than the gap 
temperature 0.330 cm below the surface. The fact that the heating distribution is 
as^iiipotic to the same value of -0.02 as seen in Figure 73 is fortuitous since 
the temperature distributions are not analogous. In run 533 plug temperatures are 
higher than gap temperatures to a depth of 1.473 cm and the negative asympote in- 
volves conduction from plug to gap over that depth. 
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The second overlap Joint analyzed was for run 528 at the location where the 
tiles form an aft facing T-slot. As noted previously for this run, location 271 
rather than 15V (Figure 66 ) was used; since by symmetry, conditions should be 
Identical at both locations and locations 156 and 272 provided additional data 
desired In the analysis. Proper analysis of the T-slot Involved formulating a 
three-dimensional model with three-dimensional view factors and Increasing the 
number of nodes. The scope of the study did not permit the formulating of such a 
model, so hypothetical conditions were analyzed using the two-dimensional thermal 
model to establish bounds for the heating at location 271. It ' 's assumed that 
the gap wall at location 271 would see an opposite wall with temperature/ view fac- 
tors representing upper and lower limits. Figure 79 Is an Isometric sketch of 
the T-slot configuration, illustrating the locations where temperature histories 
were available. It was expected that reattaching flow would result in the highest 
temperatures at location 156 and that temperatures at location 271, which sees 156, 
would be Intermediate between chose at 156 and 272. Figure 80 shows actual tem- 
perature distributions at these locations. As expected, those at location 156 are 
the highest. Near the surface, location 271 temperatures are Intermediate but 
crossover at a depth of 0.203 cm and Chen drop Increasingly below those at 272. 

This trend probably results from the radiation relief provided by the Increased 
view factor to the channel wall and, perhaps to some extent, to the bottom of Che 
In-line gap. 

Since temperature bounds were defined by measurements at locatxvins 156 and 
272 only to 0.254 cm below the surface and by measurements at locations 156 and 
271 over the rest of the gap, three cases were computed using the two-dimensional 
model and the following assumptions: 

a. Wall temperatures were assumed identical across the gap, as measured at 
location 271 and the nominal surface emissivity of 0.8 was used. 

b. VJall temperatures opposite location 271 were assumed equal to those at 272 
and the nominal surface emlsslvlCy used. 

c. Wall temperatures opposite location 271 were assumed equal to those at 156 
and an effective slot emltCance of unity was assumed at 271. 

The heating dlstrlb* tions obtained with these three assumptions are shown in 
Figure 81 . Using assumption (a) the gap heating distribution is reasonable al- 
though asymptotic to a slightly more negative value than was found in any of the 
other runs analyzed. The heating distributions obtained using assumption', (b) and 
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(c) show heating dropping to and remaining at significant negative values to the 
bottom of the gap. 

The gap-plug temperature differences are shown In Figure 82 with the plug 
temperatures being higher at all depths. Taken together and In comparison with 
the other 0.127 cm gap runs, Figures 80 through 82 Indicate that the T-slot 
configuration requires better radiation modeling and perhaps better temperature 
definition of radiating surfaces than are obtainable from the present data. 

4.2.5 Comparison of Heating Distributions - A comparison of butt and overlap 
joint heating distributions Is made In Figure 83 . Gap heating rate distributions 
for runs 542 and 533 are shown together to Illustrate the effect of gap depth. Run 
542 has a gap depth of 5.08 cm and run 533 (5.08 cm tile with 2.54 cm filler bar) 
has effectively a 2.54 cm deep gap. The curves coincide closely over most of the 
range with differences only near the bottom of the gap. 

A comparison was also made between data from this test (runs 542 and 544) with 
previous gap heating tests In the JSC 10 MW channel nozzle. Figure 84 Is adapted 
from Figure 31of Reference 1 and shows gap heating distributions on a 5.08 cm 
butt joint computed from a previous test reported in Section 4.1.2 and Reference 
1 . Data from runs 542 and 544 of the present test program have been added and 
are seen to be In good agreement. Cold gap widths In the present test correspond 
to two of those from the earlier cest. However at the test condition, the gap 
widths for the present test vary only slightly from the cold value because of the 
much lower RSI thermal expansion coefficient. 

4.2.6 Sensitivity Studies - The effects of 1; the high sensitivity of the RSI 
thermal conductivity to pressure and 2) coating thickness variations on the computed 
heating distributions were briefly examined. 

On noting the unexpectedly high plug temperatures, case 544 was rerun to 100 
seconds assuming atmospheric pressure in the test section. This rerun was made 
assuming the internal pressure in the RSI might not have adjusted to the test sec- 
tion pressure resulting in higher RSI thermal conductivity (Ref. 2, Fig. 10> As 
expected, the gap heating was not significantly affected but the RSI temperatures 
between the gap and plug (nodes 23 to 39) were appreciably higher than obtained with 
the recorded test pressures. Figure 85 shows results for this rerun case. The 
temperature at node 22 was obtained by interpolating linearly between the gap and 
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plug surface thermocouple data, and nodes 23 through 39 (Stack A) were computed by 
the program. For the test pressure (0.005 atmosphere) run, temperatures at these 
Intermediate nodes were lower than either the gap or plug temperatures. Indicating 
that the thermal properties are not consistent with the recorded data. While for 
the rerun, the resulting temperature distribution was much more realistic. 

Run 542 was rerun to 250 seconds assuming a coating thickness of 0.038 cm 
rather than 0.025 cm, resulting In normalized heating rates up to 5.72 higher In 
the upper section of the gap than reported In Figures 70 and 71. 
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4.3 Analysis of Supplemental LaRC Mach 10 CFHT Tests - Analyses were performed 
of additional gap heating data obtained in the LaRC Mach 10 CFH'^. Analyses of the 
original CFHT tests are reported in Section 4.1.3 herein. The purpose of the supple- 
mental tests was to substantiate the calibration data measured previously, to eval- 
uate the effect of zero gap width, and to determine the effects on heating In the 
tile gaps at higher Reynolds number. The tests employed the same model used In the 
previous tests, namely, a wall mounted Instrumented thin skin tile surrounded by an 
array of uninstrumented RSI tiles. The test conditions, model description, and 
data assimilated are discussed In Section 3.3 and Volume II of this report. The 
test matrix Is repeated In Figure 86 for convenience and gives run number, gap 
width, tile orientation and unit Reynolds number. 

Included in the data analysis are the following: 

a) Gap heating distributions for the supplemental tests including effects 
of gap width and in-line versus staggered tile orientation. 

b) Comparison of supplemental and original heating patterns. 

c) Effect of Reynolds number on heating rates 

c) Investigation of anomalies in heating distributions; specifically, the 
decrease in surface heating near the edges of the thin skin tile and the greater 
heating rates on top of the tile than that measured on the flat calibration plate. 

4.3.1 Supplemental CFHT Heating Analysis - The heating distributions across 
the top of the tile and within the gap are presented herein for the staggered 
tile pattern and the in-line tile orientations. Heating patterns were defined for 
both tile orientations for zero gap (no gap), 0.13 cm and 0.23 cm gap widths. The 
zero gap width was achieved by placing dental plaster in the gaps to an approximate 
depth of 0.64 cm above the top surface. The surface of the dental plaster was 
allowed to harden and then sanded until a smooth transition was achieved between 
the tiles. As for the previous CFHT tests, the in-line arrangement is achieved by 
rotating the test article 90 degrees from the staggered tile orientation (see Figure 
45 ); and for the in-line arrangement, the tiles are in-line in the axial direction 
only. In the lateral direction (normal to the flow) the tiles are staggered. The 
unit Reynolds number per meter during these tests was 3.3 x 10^. All heating data 
have been normalized by the flat plate calibration data. 

Heating patterns for the staggered tile orientations are shown in Figures 87 , 

88 and 89 . The data in Figure 87 were measured at the tile centerline (Y ■ 0) , 
while Figure 88 and 89 data were recorded at the tile edges of Y - +7.2 cc and 
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CFHT TEST 92 - SUPPLEMENTAL RUNS - FEBRUARY 1974 
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Y ■ -7.2 cm respectively. For each Y location, heating distributions are presented 
for che three gap widths Investigated. 

Similarly, re alts are presented for the In-line tile orientation in Figures 
90 (Y - 0), 91 (Y - +7.2 cm), and 92 (Y - -7.2 cm). 

Major conclusions drawn from these supplemental data are that the zero gap 
width condition shows a heating pattern very similar to that produced when a 
physlca' gap existed, heating ratios (H.^/Hpp) at the center of the tile exceeded 
the calibration plate heating, and surface heating decreases near the edges of the 
thin skin tile. The causes for these anomalies were Investigated as reported In 
Sections 4.3.4 and 4.3.5. Also, the In-line gap configuration at all three latera.' 

(Y) locations show a stronger sensitivity to gap width than does the staggered 
configuration. For both configurations, the heating pattern'! at the edges of the 
tile (Y « + 7.2 cm) exhibit a lower level of beating which Is also less sensitive 
to gap width than the centerline case. 

A comparison of the data at the tile edges (Y ■ + 7.2 cm) Indicate the heating 
patterns deviate somewhat rather than being symmetric about the centerline. For 
both tile orientations (Figures 88 and 89 for staggered and Figures 91 and 92 
for in-line), the heating at the leading edge is higher at Y • -7.3 cm than at 

Y > +7.2 cm. This deviation may be due to a slight mismatch In tile height causing 
a small step. 

4.3.2 Comparison to Original CFHT Testi - The heating distributions obtained 
ill these supplemental CFHT tests were compared to the original analyses reported 
in Section 4.1.3 of this report and also In Reference 1 . The comparison was made 

for the staggered tile configuration and tor centerline (Y > 0; heating distributions. 
For both series of tests, tile thickness was 6.35 cm and the thin skin tile had an 
edge radius of 0.3775. Likewise, the test conditions were a Mach number of 10 and 
a unit Reynolds number per meter of 3.3 x 10 . All data were corrected for thin 

skin tile specific heat variations (C ) as described in Section 4.3.5. 

P 

The comparison between r.he supplemental ana original distributions is shown in 
Figure 93 for gap widths of 0.127 and 0.229 cm. Both gap and top of the tile dis- 
tributions are shown. Again, the tile heating has been normalized to the flat plate 
calibration heating. As seen, agreement is poor between the original and supple- 
mental heating distributions, particularly on the top surface of the tile. The 
original data show a higher heating near the leading edge of the tile, while the 
supplemental test heating peaks near the center and aft edge of the tile. In 
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addition the original data are less sensitive to gap width variations. It Is also 

noted that the original data, even after correction for specific heat variations, 
sLill exhibit higher heating on the top surface than the flat plate C'..e., H^,/ 

hpp > 1.0). The reason for the poor comparison between the supplemental and original 
heating distributions is unknown; possible causes could be tile mismatch, different 
instrumentation set points, and different test techniques. 

The gap heating distributions were again compared using a different normalization 
parameter than flat plate heating. Since the heating in the gap and the heating on 
the top of the tile at the gap are interrelated, the feasibility of developing 
correlations in terms of suggested. For example, by employing this 

parameter good agreement was obtained for gap heating (both the downstream and up- 
stream face of the gap) between the supplemental and original data. These comparisons 
are made in Figures 94 and 95 , again for the staggered tile orientation and center 
line distributions. In Figure 94 , the Z distance into the gap is used as a re- 
ference, while in Figure 95 the distance down the gap is measured from the top of 
the tile and includes the edge radius. Using the latter reference improves the 
comparison, especially near the top of the gap. 

4.3.3 Reynolds Number Effect on Heating Rates - One purpose of conducting 
additional tests was to determine the effect of higher Reynolds number on heating 
in the tile gaps. The original CFHT tests were all conducted at a Reyiiolds number 
per metf'r of 3.3 x 10^, while supplemental tests were conducted at both 3.3 x 10^ 
and 7.4 X 10 . Test section Mach number was 10 for all runs. 

The CFHT side wall heating was calibrated for both flow conditions used in 
the gap heating experiments. The heating distributions (Figure 96 ) are similar. 
Although the distributions are similar, individual calibration curves were used to 
normalize heating data for each Reynolds number test condition. The heat transfer 
coefficients (h^) reported in the figure have been corrected for specific heat 
as discussed in Section 4.3.5. 

Tha effect of increasing Reynolds number on tile heating is shown in Figure 
97 for the staggered tile configuration and for the centerline (Y = 0) distribu- 
tion. The heating at the higher Reynolds number flow condition has been referenced 
to the heating distributions at a Reynolds number of 3.3 x 10^/m. For this refer- 
ence case, the data were previously reported in Figure 87. Results are presented 
in the figure for gap widths of 0.0, 0.127 and 0.229 cm. As seen, the effect of 
higher Reynolds number is to double heat transfer coefficients on the tile surface. 
Within the gap the heating increases by factors varing from 1 to 3, depending upon 
gap width and location within the gap. 
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4.3.4 Investigation of Decrease In Surface Heating at Tile Edges - The heat- 
ing data taken In the CFHT during both the original and supplemental test runs 
exhibited characteristics which have not been satislactorlly explained. This 
characteristic is a decrease in surface heating near the edge of the thin skin tile 
as seen in Figure 87 through 93 for both staggered and in-line tile configurations. 

Some of the possible reasons for this lower heating level are: 

a) Flow phenomena produced by the presence of the gaps 

b) Transient conduction effects in the thin skin tile 

c) Thickness variations in the thin skin tile 

Gap Flow Phenomena - Ten gap heating runs (See Figure 86 ) were made during 
the supplemental test program. Two of these runs were made with dental plaster in 
the gaps (zero gap width) to show the effect of gaps on tile surface heating. 

Figure 87 presents the measured heating distribution along the centerline of the 
tile for zero, C.13 cm, and 0.23 cm gap width. It can be seen that the heating rate 
decreases near the front edge of the top surface for the zero gap run as well as for 
the runs with gaps present. This suggests that flow phenomena associated with gaps 
are not the cause of the decrease in measured surface heating rate near the tile 
edge. Since the filled gap runs did not produce a flat plate heating distribution 
on the tile surface, it was decided to take a detailed look at conduction in the 
thin skin tile for a filled gap run. 

Temperature distributions along the centerline of the thin skin tile for a 
typical open gap run and the two filled gap runs are shown in Figure 98 . These 
distributions were measured at the time when heating rates were evaluated. The 
open gap distribution shows a temperature gradient between the surface and the gaps 
which is due to preheating of the tile during tha Insertion procedure. This dis- 
tribution is typical of the distributions e>isting in open gap tests. The filled 
gap was achieved by placing dental plaster in the gaps to an approximate depth of 
0.64 centimeters from the top surface. The surface of the dental plaster was allowed 
to harden and then sanded until a smooth transition between tiles was achieved. It 
is apparent from the measured temperature distributions that the dental plaster was 
much cooler than ambient when the tests were run. The f^rst filled gap run had 
measured gap temperatures as low as 6.1*C, while the top of the tile was at 27.6®C 
at its center. The second gap run (170) also had severe temperature gradients within 
the tile. These initial temperature gradients are the probable cause for failure of 
the filled gap tests to produce a flat plate heating distribution on the tile surface. 
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Transient Conduction Effects In Thin Skin Tile - Run 170 was evaluated to 

determine If correcting the measured heating rates for conduction would result in 

a uniform heating distribution on the tile surface. To evaluate the effect of 

transient conduction on the heating rate, the second derivative of temperature with 

2 ~2 

respect to distance (3 T/3X ) must be determined. The methods used to calculate 
the second derivatives are shown in Figure 99 . One method was to hand fair a 
smooth curve through the measured temperatures and numerically determine the second 
derivative using 3 points from the curve. A range of values was obtained due to 
uncertainties inherent in reeding the hand faired curve. Another approach employed 
a fourth order least squares curve fit using seven data points to determine the 
temperature as a function c£ distance. The second derivative was then evaluated 
at the center data point. The third method used was a second order curve fit using 
three data points. The second derivative was evaluated at the center data point. 

Figure 99 presents a sunnnary of the second derivatives for the surface thermo- 
couples. Also shown is the uncorrectmd and corrected nondimensionallzed heat 
transfer coefficients (H^/H^p) . These coefficients are presented in Figure 100 
as a function of distance from the center of the tile. The conduction corrected 
coefficients are shown as a range of values as computed by the three methods. The 
open symbols display the characteristic surface heating distribution that has been 
observed in all the CFHT tests. In general, the conduction correction results In 
an insignificant change in the heat transfer coefficient. The thermocouple located 
at X -7.62 cm is known to be covered with dental plaster and thus should not 
have experienced convective heating. The conduction correction tends to reduce 
the coefficient toward zero at this location. The heating at locations -7.21 cm 
and +7.40 cm is not significantly affected by correcting for conduction. In 
addition, the level of heating at these locations is approximately the same as is 
measured with open gaps (See Figure 87). 

It appears that neither conduction or gap flow phenomena (l.e., the presence 
of dental plaster) adequately explains the low heating measured at these locations. 
Correcting for conduction made a non-trlvial change in the heat transfer coefficient 
at X ■ -6.85 cm. The remaining thermocouples show very little effect of conduction 
on measured heating. 

The above conclusions are reinforced by the comparison of conduction-corrected 
and uncorrected heating rates for Run 161 shown in Figure 101 . 
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Concurrently with the above analysis, a calculation package was added to the 
data reading program used to Input test results to the Data Bank. The package 
arranges data from adjacent Instrumentation so that the second derivatives In both 
(x) and (y) directions can be computed. The correction on heat transfer coefficient 
Is Implamented as follows: 



h 


cor rected 


T -T 
AW W 




corrected " HMI ' (T -I„) I 2 2 ) 

STORAGE ^ ^ / 

only- 


Thlckness Variations In the Thin Skin Tile - The expected variation in measured 
heat transfer coefficient due to tile skin thickness uncertainties is also shown in 
Figure 100, The thickness of the tile was measured by LaRC personnel and was found 
to be within + 5 percent of the nominal thickness of .0254 cm. The dashed lines 
represent the possible heat transfer coefficient variation around an estimated sui- 
face heating level due to a 5 % uncertainty in tile thickness. By including both 
conduction effects and tile skin thickness uncertainties, it is possible to conclude 
that uniform heating occurred on the surface for all locations but X = -7.21 cm and 
+7.40 cm. However, to do so one must assume worst-on-worst conditions. 

In summary, no definite conclusion has been drawn as to the cause of the 
measured heating drop off near the edge of the thin skin tile. It has been showti 
that the heating distribution is similar for buth the open and filled gap runs. 

This implies that the presence of gaps does not cause the drop off of measured 
heatinj! near the edge of the tile. It also does not appear that conduction alone 
can explain the surface heating distribution. It should be noted that the thermo- 
couples used in these tests were not calibrated as installed, thus, the tempera- 
ture differences between thermocouples are not precisely known. Because of this, 
the conduction corrections that were calculated have an additional uncertainty. 
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It Is suggested that the filled gap runs be neglected due to uncertainties In the 
conduction correction calculations and the measured temperatures. It was also shown 
that tile thickness uncertainties do not alone account for the surface heating 
variation. It is possible by combining conduction effects and thickness uncertain- 
ties, under the most favorable conditions, to conclude that uniform convective heat- 
ing was experienced by the tile surface except for locations X « -7.21 cm and 
+ 7.40 cm. These locations appear to have an additional heat sink which has not 
been explained. 

4.3.5 Investigation of Tile and Calibration Plate Specific Heat Differences - 
Another characteristic exhibited by the heating data obtained in the CFHT tests was 
the greater heating rates on the tile top surface than that measured on the flat 
calibration plate. For both the original and supplemental tests, the measured heat- 
ing rate on the top of the tile (away from the edge) was consistently from 10% to 
20'" higher than that measured on the companion calibration plate. It was determined 
that the specific heat of the thin skin tile (Type 304 stainless steel) used in data 
reduction was approximately 20% higher than that used in reducing the calibration 
plate (Type 321 stainless steel) data. 

The calibration plate was fabricated from 0.0508 cm 321 stainless steel and 
the thin skin tile was made from a sheet of 304 stainless steel (0.0254 cm thick). 
Nominal "alues for the specific heat of both materials were used in the data re- 
duction (Cp^Q^=0.12 and Cp222“0-1019 Cal/gm°C) . To verify the correctness cf these 
values the specific heat of a coupon cut from the two test articles was measured. 

The specific heat was measured at LaRC using the same equipment and technique on 
both specimens. Figure 102 contains the measurements and specific heat ratio. 

It was necessary to extrapolate down to 300°K where the wind tunnel tests were 
conducted. The measured specific heat values used for correcting the data are 
Cp 30^=0. 112 and Cal/gm°C. This reduces the ratio by 9%. The 

computer program has been modified to automatically correct the CFHT data. 

Figure 103 thows graphically the effects of the specific heat correction on 
tile heating. The correction was made on Runs 161 and 165 which are for the stag- 
gered tile configuration and gap widths of 0.13 and 0.23 cm. As seen in the figure, 
use of the Cp correction reduces the heating on top of the tile to less than that 
measured for the flat calibration plate. Subsequently, all test data have been 
corrected. 
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SPECIFIC HEAT MEASURED ON THIN SKIN MATERIALS 
USED IN LaRC CFHT TESTS 
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4.4 Analysis of Edge Radius Tests at NASA JSC - Analyses were performed on 
data obtained from gap heating tests on thin skin tiles with various edge radii. 

The tests were conducted by C. D. Scott of NASA, Johnson Space Center in the JSC 
10 MW Arc Tunnel. Major purposes of the tests were to investigate the effect of 
tile edge radii on gap heating and to compare thin skin tile and RSI tile heating 
distributions obtained In the same tunnel. The test conditions, model description, 
and data assimilated are discussed in Section 3.4 and Volume II of this report. 


The test employed sets of thin skin metallic tiles mounted In a wedge test 
fixture inclined at 15° angle-of-attack. The test article was inserted alternately 
with a calibration panel into the flow field produced by the 20 inch diameter 
conical nozzle. Four edge radii (0.157, 0.3175, 0.635, and 1.27 cm) were paramet- 
rically tested at gap widths of 0.127, 0.254, and 0.381 cm. The joint configuration 
was an in-line butt, and the tile thickness was 4.1275 cm for all tests. 

In the following sections, results from these tests are reported and compared 
with gap heating distributions obtained in the same tunnel using RSI tiles. 

4.4.1 Analysis and Comparison - Twelve tests were conducted on the thin skin 
tiles with four edge radii and three gap widths. Heating distributions for the 
minimum edge radius tiles are compared with data from RSI tiles in Figures 104 thru 
107. The RSI tile heating distributions were computed from test data taken in the 
JSC 10 MW channel nozzle as reported in Reference 1 . Figures 104 and 105 pre- 
sent transverse gap data for gap widths of 0.381 and 0.254 cm respectively. 
Distributions are shown for the RSI tile thicknesses of 3.18 cm and 5.08 cm. The 
thin skin tile data are normalized by both flat plate and surface heating rates. 

The flat plate heating rate was measured on the calibration plate, while the surface 
heating rate was measured on the thin skin tile surface forward of the transverse 
gap. For both gap widths, the heating rates measured on the thin skin tile are 
shown to be higher than those derived irom RSI tile data. 

Figures 106 and 107 compare in-line gap heating data from thin skin and RSI 
tile tests for gap widths of 0.381 cm and 0.254 cm. The thin skin tile data are 
significantly nigner tnan the RSI tile data for che in-line case. 

A comparison of the heating to the upstream and downstream walls of a trans- 
verse gap are shown in Figure 108. These data are from the thin skin tile array 
having a corner radius of 0.318 cm and includes three gap widths (0.127, 0.254, 
and 0,381 cm). The herting on the downstream side of the gap is higher than that 
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for the upstrean side within 0.7!i cm of the surface. Below this point the heating 
appears similar. The effect of width on gap heating is greater for the down- 
stream side of the gap. Increasing gap width results in Increased heatii.g to the 
downstream wall. The effect of gap width on upstream wall heating appears to depend 
upon the distance into the gap. Near the surface and deep in the gap there is 
little affect of gap width. At Z * 0.39 and 0.55 cm the heating increases as the 
gap width varies from 0.127 cm to 0.254 cm. Further increases in gap width do 
not significantly Increase gap heating at these locations. 

4.4.2 Effect of Edge Radius on Heating to Top Surface of the Tile - The four 
tile sets tested by C. D. Scott In the 10 MW Arc Tunnel were Instrumented on the 
top surface near the gap as well as in the gap. Data from the top of the tile near 
the transverse gap were examined in Figure 109 to determine the influence of the 
gap on surface heating. For the upstream side, all the tests except one (Edge 
Radius •• 0.635 cm at a gap width of 0.127 cm) exhibited the same heating Independent 
of the edge radius or gap width. Thiji was as expected because the flow was 
supersonic and the influence of a disturbance should not be felt upstream. 

Most of the tests showed an increase in heating across the gap. Heating on 
the downstream side of the gap showed as much as a 68% increase over the upstream 
data. However, not all the data showed such as increase. The data for the sharpest 
tile at the narrow gap actually show a slight decrease on the downstream gap side. 
This could be data scatter. 

Examination of the downstream data shows no effect of edge radius on heating. 
There is however, an effect of gap width, with the Intermediate gap being highest 
at the tangency point. 
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^•5 Analysis of Single In-Line Gap Tests at Ames 3.5 Foot H.W.T. - Data ob- 
tained from single in-line gap heating tests conducted in the Ames 3.5 Foot Hyper- 
velocity Wind Tunnel (H.W.T.) were analysed to assess the effects of tile orienta- 
tion fo’' several gap width settings. Additionally, the heating patterns for filled 
and unfilled gaps were examined including the effects of correcting for thermal 
conduction. Flow conditions were varied to obtain data in the presen.;e of lamlnrr, 
transitional and turbulent boundary layers. 

Ihe test program used a thin skin model inserted into a carrier plate (See 
Figure 14 ). Joint configuration was a single, in-line gap which was tested at 
four flow orientations, three gap widths, four gap depths, and two gap lengths. 

Some 93 tests were conducted. The test conditions, model description, and data 
assimilated are discussed in Section 3.5 and Volume II of this report. A complete 
listing of the runs is given in Figure 15 . 

The analyses included evaluation of: 

a) Heating patterns for filled and unfilled gaps 

b) Correcting for conduction effects in thin skin tiles 

c) Gap heating distributions including effects of Reynolds number and gap 
width 

d) Different references for normalization of data 

e) Flow orientation effects on in-line gap heating 

4.5.1 Heating Patterns in Filled In-Line Gaps - Initial analysis of the in- 
line gap heating tests consisted of examining the data taken when the gap was filled 
with dental plaster. The test panel configuration and coordinate system is shown 
in Figure 14 . Tests were run at unit Reynolds number per meter of 1.57 x 10 , 

3.32 X 10^ and 6.58 x 10^. Instrumentation was sufficient to define axial and 
lateral heating distributions. 

The heating distribution in the axial direction (X) adjacent to a filled 0.127 cm 
gap is shown In Figure 110 for the three Reynolds numbers investigated. The data were 
measured at a lateral position (Y) of -0.51 cm. The figure shows the influence of 
boundary layer transition on heating at the higher two Reynolds numbers. The increase 
in heat transfer coefficient is very significant at the highest Reynolds number. 

The heating distribution in the lateral direction (Y) is shown in Figures 111, 

112 and 113 for unit Reynolds numbers per meter of 1.57 x 10^, 3.32 x 10^, and 6.58 x 
10^ respectively. In each figure, lateral disti Ibutions are given for five specific 
axial distances. Again, the 0.127 cm gap has been filled with dental plaster. It 
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should be noted that Figures 111 thru 113 are not consistent with the axial dis- 
tribution of Figure 110 in that different Instrumentation was employed to obtain 
the lateral distributions. However, the data presented In Figure 111 for the 
lowest Reynolds number and at Y = +0.51 cm do agree fairly well with heating mea- 
sured (Figure 110) on the other side of the filled gap. For all three Reynolds 
numbers, the lateral heating distributions appear to decrease significantly near 

the gap, suggesting a heat sink produced by the dental plaster filling the gap. 

2 2 

Both finite element analyses and transient conduction effects (3 T/3Y ) were eval- 
uated to explain this discrepancy, but neither could adequately account for the 
low heating measured near the gap. Finally, Figures 112 and 113 for the two higher 
Reynolds numbers show streamwise heating distributions indicative of boundary layer 
transition. 

^.5.2 Conduction Effects In Thin Skin Tile - The previous section has shown 
that the measured surface heat transfer coefficients decrease in the lateral direc- 
tion (Y) near the in-line gap when the gap is filled with dental plaster. Addi- 
tional analysis was performed to evaluate conduction in the test article. The data 
reduction method customarily used in thin skin model testing considers only the 
temperature rise rates and the local heat storage characteristics of the thin skin 
material. Thus the thermocouple/recording system is not generally calibrated against 
known temperatures. It was assumed that this was true for this test and that eval- 
uating conduction effects using second derivatives of measured temperature at the 
time of insertion would be futile. Consequently, ’.-as decided to calculate con- 
servative conduction effects by using temperatures one second after model in- 
sertion into the tunnel. These temperatures were calculated assuming isot.hermal 
conditions at the time of insertion and determining the temperature rise rate from 

the measured heating rates. These temperatures were then used to determine the 

2 2 

second derivative (3 T/3Y ) with respect to distance and then the transient con- 
duction effects. 

Figures 114 , 115 and 116 present the measured and conduction corrected heat 
transfer coefficients for Reynolds number per meter of 1.57 x 10^, 3.32 x 10 , and 
6.58 X 10 , respectively. In general, the conduction corrected coefficients are 
not significantly different from the measured coefficients. In addition, the 
conduction corrected heating distributions also exhibit the drop off near the 
filled gap. It was concluded thai conduction alone would not explain the lateral 
heating distribution near the filled gap edge. 
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4.5.3 Surface and Gap Heating Distributions - Surface and gap heating dis- 
tributions in the direction of flow were examined for the Ames in-line gap tests. 

Both Reynolds number and gap width effects on heating distributions were evaluated. 
The instrumentation for the runs analyzed consisted of one row of thermocouples 
along the surface of Y -0.51 cm and two rows near the top of the gap at Z ■ 0.25 
cm and 0.51 cm. See Figure 14 for the coordinate system definition. 

The measured heat transfer coefficients for a gap width of 0.127 cm and Re^/m 
of 1.624 X 10^, 3.270 x 10^, and 6.359 x ..0^ are given in Figures 117» 118 , and 
119 respectively. The surface heating data (Z«0.0) indicates that the boundary 
layer is either laminar, transitional, or turbulent depending on the distance from 
the leading edge and Reynolds number. Surface heating decreases with distance up 
to 90 cm from the leading edge for the low Re^/m case (Figure 117) and then gradually 
increases with distance, indicative of transitional flow. The heating in the gap 
at Z=0.25 cm shows a similar trend with distance suggesting that gap heating is 
Influenced by surface heating near the gap. The heating at Z=0.51 cm into the gap 
was too low to determine a trend with distance. The Re^/m “ 3.27 x 10^ data *,Fig- 
ure 118) show surface boundary layer trnasition starting at X=75 cm and ending at 
approximately X=125 cm. Again the gap heating is significantly affected by the 
changing surface environment. The heating data for Re^/m = 6.359 x 10^ (Figure 
119 ) shows a short laminar zone with transition starting at X*50 cm and ending at 
approximately X=90 cm. The heating in the gap near the beginning of the gap is in- 
creasing rapidly with distance. This suggests that the flow has not completely 
filled the gap. The heating at Z=0.25 cm is approximately equal to the surface 
heating during the transition zone. The gap heating at both Z«0.25 cm and 0.51 cm 
follows the sn’"face heating trend for the turbulent surface heating zone. Figures 
117 thru 119 illustrate the fact that gap heating is a strong function of the sur- 
face heating at the edge of the gap. The trends observed in the surface heating 
distributions are seen in the gap also. 

The heating distributions fo’" a wider gap (0.254 cm) and Re^/n> “ 3.138 x 10 
are shown in Figure 120 . The surface heating data (Z^Q.C) Indicates transition 
starting at 60 cm from the leading edge and ending at X*120 cm. Transition occurred 
farther forward for the 0.254 cm gap than for the 0.127 cm gap at the same Reynolds 
number as seen by comparing Figures 120 and 118. This indicates that the presence 
of the gap influences the surface heating distribution. The heating in the gap 
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shows a uharp increase and then decrease with distance when the surface boundary 
layer Is l.amlnar. As soon as transition starts the gap heating fellows the sur- 
face heating trend. 

4.5.4 Normalization of Heating Rates - The previouo gap ' atlng data from 
the Ames 3.5 Foot H.W.T. has been ratloed to measured surface heating on a smooth 
flat plate. However, the data presented In Figures 117 thru 120 Indicates that 
(1) the presence of gaps influences the surface heating near the gap and, (2) the 
heating In the gap Is a strong function of the surface heating near the gap. This 
suggests that the gap data would correlate better if ratloed to the measured sur- 
face heating for the run than If ratloed to the smooth flat plate heating rates. 

Figure 121 presents ^av. heating data for a gap width of 0.127 cm and Re^/m 
of 3.27 X 10 (Run 7) ratlcv to both flat plate and measured surface (Z«0.0) 
heating rates. The ratios are similar up to X»80 cm. The data ra,:loed to the 
surface heating forms a ntio which is lower than when it is ratloed to the flat 
plate heating between X»80 cn and 130 cm. This is due to earlier transition on the 
gap ..''(i-h-i than on the smooth flat plate. Figure 122 compares the two gap heating 
ratt. fo* •> gap width of 0.2''4 cm and Re^/m - 3.138 x 10^. Again the ratios 
rt &';roxim o.t I until transition occurs on the surface of ti.e tett article 

neai th . o_ \ ( , i 'tn) , 

"erniii . .UiPi* t heating data to the local surface heating results in a 

mo' ’.t, iiorr. .at^^j as a funnticn of distance where the .axterrial boundary layer is 
trarsicional c-r : i’^hulent. Because of this, it is concluded that normalizing the 
gap heaLi’.ig dac > t. . che measured surface heating near t!ie gap is the preferred form 
for correlation purposes. 

4.5.5 Flow Orientation Effects on In-Line Gap Heating - Heating data obtained 
on the Ames long in-line (or axial) gap were segregated according to laminar, tran- 
sitional or turbulent flow conditions to develop correlations for an in-line gap 
and an in-line gap at small incidence angles (y) up to 15 degrees. As discussed 
above, the presence of the in-line gap promotes boundary layer transition from 
laminar to a transitional state and finally to fully turbulent flow. Therefore, 
rather than using the calibration plate heating distribution to determine transi- 
tion, the heating distribution along the panel measured by the line of thermocouples 
at a distance of O.Sl cm from the gap was used. 

Analyses of the gap heating data for laalrar boundary layer ar.a discussed 
In the following paragraphs. Two sets of data were selected. The data in the 
first set all have a zero flow incidence engle "y" while the second set contains 
the variation In "y" up to 15 degrees. 

m 
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In-Line Gap Heating Laminar Flow, y ■ 0* - Heating measurements along the gap 
length were obtained at several depths Into the gap for three variations in gap 
width, four variations in gap depth and three levels of Reynolds number. Data for 
laminar flow conditions and zero flow incidence are presented in Figures 123 thru 
125 . In Figure 123 , heating distributions into the gap are presented for a gap 
width of 0.127 cm and tile thickness of 2.032 and 3.81 cm. As shown in this figure 
the heating level appears to Increase with Reynolds number and also with axial flow 
length. Increasl'^g gap depth decreases the heating. Some of the data for the 

2.032 cm thick tile was replotted in Figure 12A on different types of graph paper 

to emphasize the non-linearity of the data. 

The effect of gap width on the surface and gap heating distributions is ex- 
amined in Figure 125 for a tile thickness of 2.032 cm. Reynolds number and axial 
location were essentially constant for each gap width investigated. As expected, 

the data indicate that heating increases as gap width is Increased, and that heat- 

ing variation with gap width appears to be non-linear. 

In-Line Gap Heating Laminar Flow, 0 ^ Y ^ 15° - The effect of flow orientation 
angle (y) on the heating distribution in the in-line gap is shown in Figure 126 
and 127 . Figure 126 applies to the low Reynolds number tests, a gap widtn of 
0.127 cm and a tile thickness of 2.032 cm. Curves are presented using both flat 
plate heating and surface heating as the normalization parameter. As seen, the 
heating distributions do show a dependency on ^ae flow angle. However, other 
tunnel conditions and parameters such as distance into the gap, gap width, and unit 
Reynolds number are larger and tend to overshadow this effect. For example, the 
runs at y = 10° show a sharp increase in heating due to gap width especially at 
the high r Reynolds number tested, see Figure 127 . Heating on the downstream side 
of the gap is almost two times that on the upstream side. It is concluded that 
the presence of the in-line gap at a slight incidence angle cause the boundary 
layer to be more turbulent similar to a boundary layer transition phenomenon. 
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A. 6 Analysis of LaRC 8 Foot HTST Tests of Gap Heating Panel - Data from the 
LaRC 8 Foot HTST were analyzed to determine the effects of a large (11) field of 
tiles on heating patterns at critical locations on the tile when subjected to re- 
latively thin (0.61 £ <5* £ 1.43 cm) turbulent boundary layer. A sm oth plate was 
positioned in the test sled in place of the gap test panel to obtain calibration 
heating data to determine the degree of flow non-uniformity and as a reference for 
the gap heating data. The calibration data (Figure 128) showed a minimal deviation 
on the panel especially in center 25 centimeters where the thin skin tile was 
located. The heating across the test panel at various X-stations was essentially 
uniforms with a data spread around a mean. The reference heating curves used in 
the data analysis are snown in Figure 128 . Reference 3 contains oil flow 
patterns and additional heating measurements on the test sled which also indicates 
uniform flow. Flow fences at the sides of the sled were employed during the gap 
panel tests. Heating patterns were examined for both in-line and staggered tile 
patterns at three gap widths. Also the increase in gap heating caused by protruding 
tile were assessed. The panel was instrumented very heavily at the top of the gap 
so good definition of heating on the radiused portion of the tile could be measured. 
Test conditions, model >1escrlption, and data assimilated are discussed in Section 
3.6 and Volume II of this report. The central tile (thin skin stainless steel) was 
held in the same position for all tests so direct comparisons could be made as panel 
configuration was changed. 

4.6.1 Staggered Tiles - One of the most critical gap heating locations is the 
T-slot where an in-line gap terminates in a transverse gap. Indeed, high heating 
was experienced (.see Figure 129) at the stagnation point. 

The radiused portion of the leading edge of the tile experienced a 2.6 fold 
increase in heating. Instrumentation on either side of the radiused zone substan- 
tiate the high heating. Alone the Y “ 3.8 cm cut, heating on the forward portion of 
the tile was more uniform and near the outhoard edge of the tile the heating dropped 
below the smooth plate value. This suggests that the boundary layer flow from the 
upstream tile protected the downstream tile. It should be pointed out that the heating 
on the center of the tile was uniform and nearly equivalent to that measured on the 
calibration plate. This is in contrast to what was measured in the CFHT test facility. 

On the aft portion of the tile the heating pattern reversed. Along the tile 
■ centerline the heating drops off very pronounsly with distance until at the trail - 

I ing edge the heating was halved. At the outer edge of the tile, the heating re- 

I mains nearly constant except for the trailing edge of the tile where high heating 

I 18? 











REPORT MDC E1240 
JSC 096S) 


ntCOOMIVCLt. DOUGLAS ASTDO^AUTiCS COS9f*AIVV . CAST 


/MCOO>V/Wf LL OOL/OLXS ASTIfOISI/HJ . !€J.m AM 


a 


w w w ii 








RSI GAP HEATING ANALYSIS - II 

a 


VOLUME I 


REPOR r MDC.r'. 4 g 
JSC .Obi 


HEATING ON PANEL WITHOUT GAPS (LaRC 8 FOOT HIGH 
TEMPERATURE STRUCTURES TUNNEL) 



TEST RUN 

58 1 

58 3 

58 4 


COMBUSTOR 5 0, 

PRESSURE (N/m^) (DEG) 

7.10x10® 7.7 

6.55x10® 0.0 

17.9x10® 0.08 


NOTE 

0 MEASUREMENTS 
MADE at FIVE 
LATERAL STATIONS 
ACROSS PANEL 

- -DATA SPREAD 
CURVE USED AS 


REFERENCE FOR 
GAP TESTS 





4 


^ RUN 3 


183 


Figure 12b 


MCOOMMULL t>OUGLM9 A9TmOMAUT$C» COM^AJVV • KAmr 



HEATING DISTRIBUTION ON TILE IN 8 FOOT HTST 



-12 -8 -4 0 4 8 12 

J - DISTANCE FROM TILE CENTER (CM) 





i^ouaLAS >isr«?o/vAL/rics co/m#»a^v • cast 


\ 





a 


RSI GAP HEATING ANALYSIS - II 


VOLUME I 


REPORT MDC E)248 
JSC 09651 


(hl/hp? 1.3) was measured. This could be due to an enthalpy dump. Data measured 
by the same thermocouple for an in-line tile configuration did net record high 
heating, hence the energy damp effect must be real for the staggered tile configura- 
tion. 


When the free stream unit Reynolds number was doubled (Figure 130), the heating 
ratio level increased. Heating ratios increased by as mush as 30%. For all data 
analyzed the measured heating was ratloed to that measured on a smooth calibration. 
At the higher Reynolds number condition the edge experienced at heating ratio of 
2.9 and the energy damp zone had a ratio of 1.67. These high heating ratios should 
be of concern when designing a TPS. 

In-line Gap - Gap heating distributions at two stations along the in-line gap 
were examined (Figure 131) to determine the Influence of distance along the gap and 
gap width. As can be noted from the data for the -tarrower gap (W » 0.102 cm), a 
significant heating lnc:.ease was realized at the downstream station. As the gap 
width increased, heating at the upstream station increase and heating at the down- 
stream station decreased. Again these observations are characteristic of a com- 
plicated flow over the tiles and through the gap. 

Effect of Gap Width on Keating at Critical Locations - Next the heating at 
various locations on the tile was examined for a staggered tile pattern (Figure 132) 
and for an in-line tile pattern (Figure 133). These analyses were performed for 
several gap widths. As can be seen in Figure 132 for the staggered tiles, the 
leading edge (T/C 1 and 91) experience significant heating. 

Heating ritios as high as 3.2 are imposed on tiles with wide (0.3 cm) gaps. 
Nominal smooth plate heating experienced at the tile center for all gap widths and 
the heating at the trailing edge (T/C 65) remains a constant (h.j,/hj,p 0.5) inde- 
pendent of gap width. The energy damp at T/C 17 decreased slightly for the wider 

gaps (h_/h decreased from 1.3 to 1.1). 

T FP 

4.6.2 In-Line Tiles - The in-line tile pattern also produced Interesting heat- 
ing patterns. As can be seen from Figure 133 the high heating location was on the 
outboard region of the tile just aft of leading edge rather -han on the leading edge 
at the tile's centerline as was the case for the staggered tiles. The forward cor- 
ner (T/C 6) of the tile experienced heating ratios of approximately ] 85, with only 
a modest increase as the gap was opened. At T/C 66, which was on ' tile's top 
near the leading edge, high heating ws. measured with a stronger dependency on gap 

width (h_/h__ of 1.86 to 2.2). 
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Heating on the leading edge of the tile at the centerline (T/C 91) also In- 
creased with gap width but at a much lower level. For the smaller gap (0.177 cm), 

T/C 91 recorded less than the smooth plate heating, but showed a strong sensitivity 
to gap width with a 30% Increase In heating over the smooth plate for the widest gap 
(W = 0.4 cm). Apparently this region of the gap Is sheltered for narrow gaps but 
becomes subjected to the onslaught of the flow for wide gaps. The center of the tile 
receives, essentially, heating as If there were no gaps. The heating drops off 
significantly rear the trailing edge of the tile as was the case for the staggered 
tile. This could be due to a local thickening of the boundary layer upstream of a 
transverse gap. From the data obtained fron these tests It should be possible to 
formulate a flow field model to describe the heating on the top surface of the tile. 


Figure 134 shows such a flow model. The location and strengths of the crests would 
be dependent on gap width and the ability of the adjoining gaps to alter crest 
strength. Of course Reynolds number, Mach number and other flow field para- 
meters would characterize the nominal boundary layer. 

Steps - Both In-line and staggered tile patterns were tested with the 
tiles being flush and with the central tile protruding 0.254 cm Into the flow. As 
can be readily seen from Figure 135 , heating on the leading edge of the tile and the 
forward facing wall experience very high heating. This data Indicates that it makes 
little difference if the tiles are In-line or staggered; essentially the same high 
heating is experienced when there is a 0.254 cm step. One exception is at the center 
of the leading edge — the staggered tile receives only a 3.2 heating factor while the 
in-line tile receives a 3.45 factor. It should be pointed out that the severe heating 
has "washed" out at the center of the tile where smooth plate heating occurs while 
the trailing edge is relatively cool (h^/hpp 'v 0.5). 

4.6.4 Sp ecific Heat and Conduction Corrections - The preceding analyses were 
accomplished using data as received or. the data tape from the 8 Foot HTST. Examina- 
tion of the temperature histories such as shown in Figure 136, indicate that there 
were temperature gradients across the thin skin tile at the time heat transfer coef- 
ficients were originally computed. The data reading program was modlf.lcd to account 
for thermal conduction In two orthogonal directions. In some Instances the heating 
rates changed by as much as 3.8%. Of course, where there was little or no tempera- 
ture gradient, the heating rate remained unchanged. Figure 137 shows the change In 
a heating dijtr.outlon on the tile when considering thermal conduction. 
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Figure 134 
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Specific heat measurements on the thin skin tile material and the calibration 
plate were performed at the conclusion of the test program. This Information (Figure 
138) was furnished by 1. Weinstein of LaRC so that the data could be corrected. 

The specific heat of the thin skin material was 5.6% higher than that of the calibra- 
tion plate at room conditions (300*’K) . In the original data reduction the same 
specific heat was used for both materials. Hence the heating ratio (h^/h^p) was 5.6% 
low. The Data Bank was corrected for both heat storage and thermal conduction effects. 


4.6.5 Gap Heating Distributions - The staggered and In-line tile patterns 
produce significantly different heating distributions n the downstream wall of a 
gap. Figure 139 shows the sensitivity of gap heating to gap width and Reynolds 
number for the staggered tiles. 

The point analyzed in Figure 139 experiences "T-slot" flow/heating. From thi> 
limited test data, it appears that the heating in the gap is proportional to 

Reynolds number to a power greater than unity. Also the test data indicates chat some 

relief in heating occurs as the gap is opened beyond 0.305 cm. The heating (Figure 
140) in the transverse gap for the in-line tiles is much lower than for the staggered 
tiles. This is because the upstream tile shelters the gap. Gap heating Increases 
with gap width. 

This data shows even for very small gaps that the edge radius experiences 

heating ratios 1.5 or higher. The gap heating for the in-line tiles drops off more 

rapidly than does the measurements for the staggered tiles. The in-line data 
suggests an exponential function of (Z) to correlate the heating. 
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4 . 7 Analysis of RSI Tile (0.635 cm edge radius) Tests In the AMES Arc Tunnel 
Turbulent Duct - Temperature response data were collected on Silica RSI tilf 
subjected to high enthalpy turbulent flow. The details of the test facili and 

the configuration of the test article were described in Section 3.7. The rei>t 
item temperatures were used to determine sensitivity of the heating rates to gap 
width and were input to an inverse heat transfer model to find the local heat flux. 

Due to misreporting of the thermocouple locations the original attempts to 
analyze the data proved difficult. However, radiographs test panel showing 

the exact Instrumentation locations led to more consistent data analyses. Figure 
25 shows the original and corrected thermocouple locations. Figure 141 is a copy 
of the x-ray of the tile. 

Temperature histories of several thermocouples are presented in Figures 142, 14 J 
and 144 for gap widths cf 0.0, 0.127, and 0.180 cm, respectively. These data are 
from tests where T/C //I was located in the upstream position. The figures show 
the upstream tile (T/C #4, 6 and 7) running slightly hotter than the do\mstream tile 
(T/C #9, 11, and 12). Also, there is a different thermal response pattern for T/C 
y/4, and #9. Several, theories were advanced pertaining to this inconsistency, e.g. 
differences in local thermal properties or surface Irregularities. A mold of the 
tile surface in the region of the gap showed that T/C #4 protruded through the RSI 
causing a bump on the surface. Boundary layer disturbances caused by this rough- 
ness element could be the cause of the more rapid heating at this location. 

Figures 145 and 146 shows the effect of gap width on measu >d temperatures at 
200 seconds in the test for the two flow directions. Thermocouples 4, 6 and 7 
indicate conslstantly higher temperatures than thermocouples 9, 11, and 12, 
whether they are in the iq>stream or downstream position. Examination of T/C #4 
and 9 in Figures 145 and 146 indicates that T/C (19 gives consistant trendu in both the 
upstream and downstream position, while T/C (iU gives rather unsystematic results. 
Therefore, it was decided to use T/C to derive gap heating rates and to blame 
surface irregularity for the T/C #4 errantry. The shift in level is due to the 
difference in T/C depth into the slot which was found from the tile x-rays. 

A thermal model of the gap and its edge radius was setup to determine effects 
of edge radius on gap heating using thermal response data measured on the HRSI tiles. 
The details of the thermal model are Illustrated in Figure 147. The model uses hall- 
mark dimensions such as V(2) to define edge radius from which node dimensions are 
computed. 
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This facilitates comparison of heating data. This thermal model was used to 
analyze the test data obtained from the AMES Turbulent Duct Facility and for the 
computation o ' TPS thermal response during entry. 

Figure 148 shows the calculated heating distributions for the downstream face 
of the transverse gap for four gap widths. The data shows a consistent Increase 
In heating with gap depth. The heat flux predictions along with the data were 
used In a sensitivity study to determine that there was jin error In the thermo- 
couple locations. 

Figure 149 Indicates the sensitivity of heating to gap depth and width. Also 
shown on this figure Is the error caused by Inputlng erroneous Instrumentation 
locations Into the data reduction analyses. Figure 150 shows the heating distribution 
In the gap for four gap widths, calculated using the corrected data. Again, heating 
decreases with depth and Increases with width. 


207 


MCOOMMCCI. OOfJOtAS ABTmOMfkUTICS • «A«T 



GAP' ^SURFACE 



BUTT JOINT HEATING DISTRIBUTION HRSI TILES 
AT TRANSVERSE GAP 

• AMES 20 MW 2 X 9 INCH TURBULENT FLOW 
DUCT FACILITY (TEST 3) 

• TILE EDGE RADIUS = 0.635 cm 

• DOWNSTREAM WALL OF TRANSVERSE GAP 

• DISTRIBUTION AT 160 SEC TEST TIME 



1 o 









RSI GAP HEATING ANALYSIS - 11 

I i VOLUME 1 


REPORT MDC EU48 
JSC 09631 


SENSIT'VITY OF T/C LOCATION 
ON HEATING DISTRIBUTION AT TRANSVERSE GAP 
BUTT JOINT - HRSI TILES 

^ AMES 20 MW 2 X 9 INCH TURBULENT FLOW 

# DUCT FACILITY (TEST 3) 

^ TILE EDGE RADIUS = 0.635 cm 

# DOWNSTREAM WALL OF TRANSVERSE GAP 

# DISTRIBUTION AT 160 SEC TEST TIME 

W 


FLOW ^ “Ml-*- 



209 


Figure 149 


MCOONNmt-L. DOUCI.A9 ASTfrO^AtlTiCS COnAr»AMV • KA»T 









1 




MCOONMCLl. DOUGLAS ASTDOfktAUTiCS CO^A^AM^ • MAST 



REPORT MDC E1248 
JSC 096S1 


Q 


RSI GAP HEATING ANALYSIS - II 

VOLUME I 




A. 8 Boundary Layer Flow Conditions - Boundary layer analyses were performed to 
characterize the environments at the various tunnel test conditions. Basic 
boundary layer parameters were computed for use In correlating the heat transfer 
results from the test facilities. These parameters will aid In extrapolating the 
test results to Shuttle flight conditions. The basic boundary layer parameters 
determined for the various test conditions are: 

1) Type of boundary layer (laminar, transitional, turbulent) 

2) Local unit Reynolds number (edge condition) 

3) Local Mach number (edge condition) 

4) Displacement thickness (6*) 

5) Momentum thickness (9) 

6) Laminar sublayer thickness for turbulent boundary layers (6^) 

7) Temperature ratio across the boundary layer 

A summary of the boundary layer analyses that were performed is shown in Figure 
151. The details of most of these analyses are contained in Reference 1 . Data for 
the LaRC 8 foot HTST tests and the Ames 20 MW Turbulent Duct tests are contained 
herein. 

Figure 152 contains the displacement (6*) and momentum (0) thicknesses along 
the HTST panel holder and Figure 153 contains the velocity profiles which were d 
to obtain sublayer thickness. The displacement thickness on the thin skin tile 
ranged from 0.61 cm (0.24 inch) to 1.43 cm (0.56 Inch) and the momentum thickness 
ranged from 0.152 cm (0.060 inch) to 0.214 cm (0.084 Inch). A nominal momentum 
thickness (9 - 0.182 cm) was used In the data bank for all tests. It was determined 
that the displacement thickness and sublayer thickness at the center of the tile 
could be expressed as a linear fmiction of angle of attack (a) and combustor 
pressure (P). 


Displacement thickness; 

6*(inch) - 0. 32-0.01674(0-15) -3. 53636xl0"^ (P-1000) 


Sublayer Thickness; 


£n|6^(10 FT) 


-4, 


- £n 5.8-0.07017(0-15) -4.358x10 (P-1202) 


where o is In degrees 

P is combustor pressure (PSI) 


211 


/MCOO/VWfLL OOUOLM9 MmTmOMMUTiC* • mAMT 





FACILITY BOUNDARY LAYER ANALYSES 







MCOOMNCLL. tJOLfOL/tm J»ST/rOMjHJTtCS COI¥I^AMY • CAST 







VELOCITY PROFILE ON HTST TEST SLED 



DISTANCE FROM SLED - FT X 10 






REPORT MDC El 248 
JSC OViSI 


& 


RSI GAP HEATING ANALYSIS - II 


VOLUME I 


These funcClone were Incorporated Into the Data Bank computations. It should be 
noted that these equations were developed to describe 6* and 9 for a limited set of 
operating conditions In a particular test program and therefore should not be con- 
sidered as general correlations. 

Details of the Ames Turbulent Duct boundary layer are contained In Reference 

4 . Pitot pressure surveys had previously shown the total boundary layer thickness 

to be about 20 mm for a total enthalpy of 3.7 MJ/kg and a static pressure of 
3 2 

3.5 X 10 N/m . Assuming that the Crocco relationship Is applicable the momentum 
and displacement thicknesses are about 2 mm and 5 mm, respectively. 
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4 . 9 Comparison of Gap Heating Data from Arc Tunnel and Wind Tunnel Tests - 
The comparisons performed during Phase ' of the heating distributions obtained from 
arc tunnels and wind tunnels Is still valid. These comparisons are contained herein 
and additional comparisons have been made during Phase II . Comparisons were per- 
formed for both transverse gaps and In-line gaps. The gap heating data which have 
been assimilated come from bith an arc tunnel and wind tunnels. Arc tunnels pro- 
vide high energy (temperature and enthalpy) flow to test articles while wind tunnels 
provide a considerably lower energy flow. The effect of this difference in flow 
energy on gap heating was investigated by comparing data from arc and wind tunnels. 

4.9.1 Transverse Gap Heating Comparisons - Figure 154 summarizes the gap heat- 
ing test environments and transverse gap geometry for which data have been assimilated. 
As can be noted, the wide variety in conditions makes direct comparisons difficult. 

The check marks denote the data which were selected for comparison. Two gap widths 
from the JSC 10 MW arc tunnel tests were compared with the wind tunnel data. These 

gap widths bouad the selected gap widths from the wind tunnels. The arc tunne]. had 

the lowest freestream unit Reynolds number of any of the facilities. Therefore, the 
lowest Reynolds number data available from each v;ind tunnel facility were selected 
for comparison, n laminar boundary layer existed in some of the arc tunnels and 
in some of the wind tunnel tests. 

Gap heating data from each wind tunnel test were compared Individually with 
data from the arc tunnel test in Figures 155 thru 162. A comparison of gap heating 

data from the Ames 3.5 Foot Hypersonic Wind Tunnel and the JSC 10 MW Arc Tunnel is 

shov,n in Figure 155 . Both tests were conducted in a laminar boundary layer environ- 
ment with similar edge Mach numbers. Generally good agreement exists in the level 
of heating down the gap, although the shapes of the heating distributions are dif- 
ferent. Figure 136 presents the comparison of gap heating distributions for the 
arc tunnel and the "freestream" tests in the LaRC Variable Density Tunnel (VDT) . 

Both tests were run with a laminar boundary layer over the test article. The Mach 
number in the V.D.T. was 8.0 while the Mach number in the arc tunnel was 4.2. Also, 
the freestream unit Reynolds number is considerably higher in the wind tunnel than 
in the arc tunnel. The wind tunnel data agree well with the arc tunnel data at a 
depth into the gap of the 0.7 cm and below. The only data taken higher in the gap 
were taken at 0.5 cm and here the dimensionless heating was considerably higher 
than in the arc tunnel. Figure 157 presents the comparison of gap heating obtained 
using two different panel designs (a HCF panel and a thin skin tile set) tested in 
tne JSC 10 MW Arc Tunnel. The HCF panel was tested In the wall of a channel nozzle 
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SUIVIMARY OF GAP HEATING TEST 
ENVIRONMENTS AND GEOMETRY 

(TRANSVERSE BUTT JOINT) 


TESTING FACILITY 

TEST ARTICLE 
POSITION 

GAP DEPTHS 
(CENTIMETERS) 

GAP WIDTHS 
(CENTIMETERS) 

M 

m 

Re./m 

8.L. 

STATE 

JSC 10 MU ARC 
TUNNEL (TP2)* 

CHANNEL 
NOZZLE UALL 

/ 3.18 
5.08 
6.35 

.074 

/.204 

/.333 

.714 

✓4.2 

✓.06X10^ 

LAMINAR 

AMES 3.5 FOOT 
WIND TUNNEL 

FREE STREAM 

1.016 
2.032 
/ 4.064 

✓.254 

✓5.1 

n 

LAMINAR 

LANGLEY VARIABLE 
DENSITY WIND TUNNEL 

FREE STREAM 

/ 2.54 

.159 

✓.318 



✓ 8 

✓2.3X10® 

21.8X10® 

LAMINAR 

LANGLEY VARIABLE 
DENSITY WIND TUNNEL 

TUNNEL WALL 

✓ 2.54 

.159 

✓.318 

✓ 8 

✓1.16X10® 

41.4X10® 

TURBULENT 

LANGLEY CONTINUOUS FLOW 
HYPERSONIC WIND TUNNEL 

TUNNEL UALL 

/ 6.35 

.13 

✓.23 

.46 

.71 

✓10 

✓3.28X10® 

TURBULENT 

JSC 10 MW ARC TUNNEL 
WEDGE 

FREE STREAM 
(15° ANGLE OF 
ATTACK) 

V4.128 

.127 

/.254 

.381 

✓9.6 

✓1.74X10^ 

LAMINAR 

LANGLEY 8 FOOT HIGH 
TEMPERATURE STRUCTURES 
TUNNEL 

FREE STREAM 
{0°, 7.5°, 15°, 
ANGLE OF ATTACK) 

'/6.35 

0.0 

0.10 

✓0.18 

✓0.30 

0.41 

6.2 

7.0 

✓ 2.0X10® 
4.8X10® 

TURBULENT 

AMES 20 MW 2X9 INCH 
TURBULENT FLOW DUCT 
FACILITY 

FREE STREAM 

✓5.08 

■ 

✓3.5 

✓0.3X10® 

TURBULENT 


* TP2 - TEST POSITION NUMBER 2 
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and the thin skin tiles were mounted In a wedge test fixture which was swung Into 
the plume of a free expanding jet. Flow over both test articles was laminar but the 
Reynolds number and boundary layer thickness was much less for the wedge. This could 
be the reason for the higher heating measured uslnt the thin skin tile. 

Gap heating data were also taken in the wall of the V.D.T. to expose the test 
article to a turbulent boundary layer. These data are compared with the arc tunnel 
data in Figure 158 and are considerably higher than the laminar arc tunnel. Figure 
159 compares the arc tunnel data and data taken in the wall of the Mach 10 Continuous 
Flow Hypersonic Tunnel (CFHT) . The boundary layer in the CFHT tests was turbulent 
and the data for these tests are again higher than the arc tunnel data. 

The 8 Foot HTST is a combustion driven facility and the flow over the test 
panel was turbulent with a relatively thin boundary layer (0.61 £ £ 1.43 cm). 

Figure 160 Is a comparison between HTST data and the channel nozzle data obtained 
at JSC. Again the turbulent heating for the turbulent flow was higher but a cross- 
over occured between 0.5 and 0.75 cm into the gap. The heating distributions from 
the HTST were also compared (Figure 161) with another turbulent test program per- 
formed in the CFHT. The boundary layer in the CFHT was much thicker in the HTST 
tests. The heating in the CFHT was higher also. 

Finally the test results for the large edge radiused tiles tested in the JSC 
10 MW Arc Tunnel using the wedge and the Ames Turbulent Duct Arc Tunnel is compared 
in Figure 162. Both these facilities are Arc Tunnels and hence have high enthalpy 
content. The wedge model employed thin skin tiles whereas the duct model used RSI 
tiles. The wedge model had laminar flow and the duct produced turbulent flow. The 
tile exposed to the laminar flow experience a much higher heating than the one in 
the turbulent boundary layer. Apparently the trend to high gap heating in a turbu- 
lent flow is reversed when dealing with large radiused tiles. Additional test pro- 
grams need to be performed to substantiate this trend. 

4.9.2 In-Line Gap Heating Comparisons - Figure 163 lists the data assimilated 
and used to compare heating measurements obtained in the arc, wind and combustion driven 

facilities. The same criteria for selecting data employed in the transverse gap 
applies. Three different test programs were conducted in the JSC 10 MW Arc Tunnel. 
Figure 164 shows the comparison between results using two RSI materials (HCF a 
mullite RSI and Ames silica RSI) and the thin skin tile tested in the wedge. The 
two results from the two RSI materials agree very well whereas the thin skin data 
is higher. The RSI materials were tested in the same channel nozzle whereas the 
thin skin tile used a wedge which had the thinner boundary layer. Figure 165 is 
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SUMMARY OF GAP HEATING TEST 
ENVIRONMENTS AND GEOMETRY 

(IN-LINE BUTT JOINT) 


r 

TESTING FACILITY 

TEST ARTICLE 
POSITION 

GAP DEPTHS 
(CENTIMETERS) 

GAP WIDTHS 
(CENTIMETERS) 

M 

m 

Re^/m 

B.L. 

STATE 

JSC 10 MW ARC TUNNEL 

CHANNEL 

^3.18 

0.074 

✓4.2 

✓0.06x10® 

LAMINAR 


NOZZLE WALL 

5.08 

/0.204 






6.35 

/0.333 







0.714 




AMES 3.5 FOCT WIND 

FREE STREAM 

1.016 

/0.254 

✓5.1 

✓1.6x10® 

LAMINAR 

TUNNEL 


2.032 



2.6x10® 




/4.064 



4.4x10® 


LANGLEY CONTINUOUS FLOW 

TUNNEL WALL 

/6.35 

0.13 

✓10 

✓3.28x10® 

TURBULENT 

HYPERSONIC WIND TUNNEL 



✓0.23 







0.46 







0.71 
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JSC 10 MW ARC TUNNEL 

FREE STREAM 

/4.128 

0.127 

✓9.6 


laminar 

WEDGE 

(15" ANGLE OF 


✓0.254 


■■ill 



ATTACK) 


0.381 


■■■ 


LANGLEY 8 FOOT HIGH 

FREE STREAM 

✓6.35 


6.2 

✓2.0x10® 

TURBULENT 

TEMPERATURE STRUCTURES 

(0", 7.5", 15" 



/ 1 



TUNNEL 

ANGLE OF 



T 




ATTACK) 



7.0 

4.6x10® 





0.41 




JSC 10 MW ARC TUNNEL 

CHANNEL 

2.54 

0.127 

y3.7- 

;4.5- 4 

LAMINAR 


NOZZLE WALL 

/5.08 

✓0.254 

4.5 

8.2x10 


AMES 3.5 FOOT WIND 

FREE STREAM 


P.O 

/5 1 

/l. 64x10? 

LAMINAR 
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0 ’27 
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/0.254 


6.56x10° 
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Figure 163 
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a comparison between the HCF test results and the Ames 3.5 Foot H.W.T. test results. 
The wind tunnel yielded slightly but consistently higher heating ratios. The shapes 
of the heating distributions are similar. Figure 166 compares the effect of 
boundary layer state on in-line gap heating using results from the turbulent CFHT 
and the laminar JSC Duct. This data, too, shows a higher heating when the boundary 
layer is turbulent. 

The in-line gap heating obtained from the 8 Foot HTST test shows a ignificant 
higher heating as evident from Figure 167. The HTST heating pattern is also much 
higher than that obtained from the CFHT which had a much thicker turbulent boundary 
layer. It should be remembered that a field of eleven tiles was used in the HTST 
tests whereas the HCF tests employed four tiles. The flow over a tile field could 
also contribute to Increased heating. Further investigations should be conducted 
to determine if the in-line gap does receive this high heating because it will 
definitely affect TPS design. 

The difference in enthalpy between arc tunnels and wind tunnels does not .ippear 
to significantly affect gap heat transfer when normalized by the reference surface 
heating rate. The laminar wind tunnel data agree reasonably well with the laminar 
arc tunnel data. However, the boundary layer state does affect the heating in gaps. 
The data indicate that turbulent boundary layers result in higher dimensionless 
heating in gaps than laminar boundary layers. These conclusions apply to in-line 
gaps as well as transverse gap with a noticeable exception where a relatively thin 
boundary layer maintains high heating deep within the gap . Data also indicates that 
laminar flow produces higher heating than turbulent flow when the tiles have large 
edge radius. 
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5.0 DATA CORRELATION 


• The gap heating data were correlated in terms of gap dimensions, location of 

i 

[ gap, location In the gap aud boundary layer parameters. A general data management 

I system was set up so a common approach could be applied to data from each available 

I source. Test results were ordered, and combined with physical dimensions, instrument- 

j ation coordinates and boundary layer parameters to form a data bank. Over 24,000 

; individual heat transfer measurements were assimilated into a data bank, which is 

functional on McDonnell Douglas computers. This data management system permitted 
quick access to data sets with similar attributes and direct input to a MRA (Multiple 
Regression Analysis) computer program. Correlations were obtained for transverse 
gaps, in-line gaps, flow angularity, edge radius, and steps in the presence of both 
' laminar and turbulent boundary layers. 

5.1 Correlation Method and Data Bank - Test information supplied by each test 

facility was received on magnetic data tapes, computer tabs, work sheets, and facility 

test reports. Format and methods of transmitting data were suggested to each facil- 

ity and where possible these recommendations were incorporated to be consistent with 

. what was most convenient for that facility. Data from each facility were processed 

into the data bank. 

I A procedure was set up for ordering data and combining of data with flow field 

I parameters. In the following paragraphs the mechanics of the correlation procedure 

are describee with the aid of a data handling flow chart. Also included in the pro- 
cedure are dlscrimindcors used to select data for correlation using the Multiple 
Regression .Analysis program. 

5.1.1 Data Correlation Procedure - Each piece of gap heating data incorporated 
in the data bank was assigned 24 attributes which provided traceable information 
about its origin, instrumentation location in the joint, heat transfer parameters 
and boundary layer parameters. The assigned 24 attributes are listed in Figure 168. 
Traceable information about test program origin, run number and instrumentation 
designation constitutes the first attribute word. The system is formulated so that 

1 information from other tests can be added to the data bank. Information about the 

boundary layer flow over the RSI joints was combined with instrumentation location 
I and the gap heating data to complete the data bank. The 24 attributes are also used 

to select a particular set of data for correlation. The 24 attribute Information is 
stores on magnetic data tape. 

The functioning of the data management system is Illustrated in Figure 169 . For 
example, results from the CFHT tests stored on a data tape are combined with test 
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24 ATTRIBUTES WORDS ASSIGNED 
TO EACH GAP HEATING DATA POINTS 


word! TYPE 


1 


REAL 


2 

INTEGER 

3 

INTEGER 

4 

REAL 

5 

REAL 

6 

REAL 

7 

REAL 

8 

REAL 

9 

REAL 

10 

REAL 

11 

REAL 

12 

REAL 

13 

REAL 

14 

INTEGEf 

15 

REAL 

16 

REAL 

17 

REAL 

18 

REAL 

19 

REAL 

20 

REAL 

21 

REAL 

22 

REAL 

23 

REAL 

24 

INTEGEf 


XX 

r 


OYYY 


zzz 


•T/C OR CHANNEL NUHBER 
RUN NUMBER 


■TEST NUMBER 


TEST NUMBERS THAT HAVE BEEN ASSIGNED 
XX = 1 . JSC WEDGE TESTS CONDUCTED JR MDAC-E 
XX - 2, JSC CHANNEL NOZZLE TESTS CONOIKITED FOR HDAC-E 
XX « 3. OTHER JSC CHANNEL NOZZLE TESTS 

XX - 4. JSC, HEDGE EDGE RADIUS. THIN SKIN, C.SCOH 
XX ■ 10, LaRC, CFHT 


GAP CONFIGURATION 

1 » BUn, 2 • CONTOURED. 3 = OVERLAP, 4 - INCLINED 
INSTRUMENTATION LOCATION 

1 • UPSTREAM SIDE OF GAP, 2 - DOWNSTREAM SIDE OF GAP, 
4 » STAGNATION, 5 • TILE TOP. UPSTREAM SIDE OF GAP, 6 
COORDINATES OF AN INSTRUMENTATION POINT, z = 0 
AT TOP SURFACE OF THE TILE, x = DISTANCE DOWNSTREAM 
FROM CENTER OF TILE, y (RIGHTHAND RULE), (cm) 

DISTANCE FROM LEADING EDGE OF EACH TILE (cm) 


XX • n, LaRC, MACH 8 V.D.T. 

XX • 12, AMES TURBULENT DUCT TEST 

XX • 13, LaRC 8 FOOT HTST 

XX • 14. SO MI WEDGE TESTS 

XX • 15. AMES 3.5 FOOT HHT 

XX - 16, AMES 3.5 FOOT HWT, IN-LINE GAP 


3 - IN-LINE 

» TILE TOP DOWN SIDE OF GAP 


VN. 


FLOW ORIENTATION (RADIANS) 

GAP WIDTH (cm) 

STEP HEIGHT (cm) 

GAP FLOW LENGTH (cm) 

TILE THICKNESS (cm) 

TILE PATTERN, 0 • STAGGERED, 1 • IN-LINE 
LOCAL MACH NUMBER 
REYNOLDS NWBER/METER 
MOMENTUM THICKNESS (cm) 

DISPLACEMENT THICKNESS (cm) 

E UB-LAYER THICKNESS (cm) EDGE RADIUS FOR CARL SCOH ONLY 
EAT TRANSFER COEFFICIENT (h), (K6/M^SEC) 

J TANTON NUMBER 

wal/% TEMPERATURE RATIO ACROSS BOUNDARY LAYER 
OUNDARY LAYER STATE. 1 > LAMINAR, 2 - TRANSITIONAL. 3 • TURBULENT 
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matrix information, T/C coordinates and boundary layer flow field parameters in the 
"RSI C" program to generate a 24 word attribute tape. Companion tapes from the 
other tests are prepared in a similar manner. These tapes are tnen processed by 
another program which selects data according to a list of discriminators specified 
for a particular type of gap to be analyzed. Information from tae JSC tests con- 
tained on data cards can be read by either the "SELECT" program or by the Multiple 
Regression Analysis (MRA) program. The MRA program processes the selected data 
and determines the best fit for candidate correlation equations. 

A list of 35 discriminators was prepared for selecting data fro:- the 24 attribute 
tapes for correlation. Figure 170 shows an input form for the "SELECT" program. 

Upper ond lower limit discriminators were used to facilitate the selection process. 
Tabulations as well as tapes were generated containing the selected data. The 
tabulated information was valuable for identifing trends, anomalies and for checking 
data. 

Because of the large amount of heating data available for correlation, an auto- 
mated multiple regression technique was used to obtain consistent-nonbiased cor- 
relation equations. The step-wise Multiple Regression Analysis (MRA) computer program 
(Reference 1) provided information as to the adequacy of the candidate correlation 
function and the equation coefficients. A modified version of the MRA computer pro- 
gram which accepts information stored in the data bank was used for this study. The 
principal modifications included accepting information stored in the data and eux- 
iliary statistical analyses. The MRA computes a series of multiple linear re- 
gression equations in a stepwise manner. At each step, one parameter is added to 
the equation. The variable added is the one which makes the greatest reduction in 
the variance about the mean. Equivalently, it is the parameter which, if it were 
added, has the highest "F" ratio. Figure 171 lists the form of the correlating 
equation(s) and the statistical parameters use by the MRA to obtain the most 
appropriate correlation. 

Local heating at the gap ratioed to undisturbed flat plate heating or ratloed 
to the value on the tile top near the gap edge was designated as the dependent vari- 
able for all correlations. Correlations were obtained in terms of natural logarithms 
of the dependent variable because gap heating experiences a decrease of several orders 
of magnitude with distance into the gap. For the JSC 10 MW data rhe ratio was formed 
using measured convective heating rates. For all other tests the ratio utilized 
measured convective heat transfer coefficients. The heating rate ratio and heat 
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VARIABLES AVAILABLE FOR SELECTING DATA 


HEATING 

ANALYSIS 

rooiLiM RSI CAP HE/.TINC SCREENING 

COATUTEN SCIfNCtS 

ftcoursYon 0M*Y 

SELECT INPUT FORM o»u i o, 


t 

DESCRIPTION (DO NOT KEY PUNCH) OIMEf'SIOH 


I'.ISTI 



CAbE - 

CASE NUMBER 


»jFn >T. - 

NUMBER OF FILES TO BE SKIPPED Of TAPE 10 KOk PLACEMENT OF CURRENT CAS' 


XLJW - 

X-LOWER LIMIT (cm) 


XlllGH - 

X-UPPER LIMIT (cm) 


Y!.6W - 

Y-LOWER LIMIT(crI 


YliIGH ■ 

Y-UPPER LIMIT (cm) 


ZLtfW - 

Z-LOWER LIMIT (cm) 


ZHICM • 

Z -UPPER LIMIT (cm) 


XBLiJW - 

X LOWER LIMIT (cm) 


XBKICH « 

X-UPPFR LIMIT (cm) 


YBL0W - 

Y-LOWER LIMIT (cm) 


YHHIGH - 

Y-L'PPLK LIMIT (cm) 


AUilAL - 

FLOW ORIENTATION - LOWER LIMIT (RADIANS) 

j ALPILMI - 

Fl.OW ORIENTaTION - UPPER LIMIT (RADIANS) 


CAI'UX ‘ 

GAP WIDTH - LOWER LIMIT (cm) 


CAJ’U'H - 

CAP WIDTH - UPPER LIMIT (cm) 

-H 

STIPL - 

STEP HEIGHT - LOWER LIMIT (cm) 


STEPH - 

STEP HEIGHT - UPPER LIMIT (cm) 


GATFL * 

GAP FLOW LENGTH LOWER LIMIT (cm) 


C.ATFH - 

CAP FLOW LENGTH • UPPER LIMIT (cm) 


THKL - 

TILE THICKNESS - LOWER LIMIT (cm) 


THKil • 

TILE THICKNESS - UPPER LIMIT (cm) 


IPATNdl - 

TILE PATTERN: 0-STACCERED, I-IN-LINE I PATH (9 ) 


ibUIi - 

BOUNDARY LAYER STATE: I-LAMINAR, 2-TKANSITIoNAL. >TURBUL£NT, 191.(9) 


ICAPC(l) - 

CAP CONFIGURATION: 1-BUTT. 2-CONTOURED. 1"OVERLAP. A*INCLINED ICAPC(9) 


ICL0C(1) 

CAP LOCATION: 1-UPSTREAM SIDE OF CAP. 2-DOWNSTREAM SIDE OF CAP. 



3-PARALLEL. 4-STACNATION . 5-TILE TOP. UPSTPJAM SIDE OP 



CAP. 6-TILE TOP. DOWNSTREAM SIDE OF CAP ICL0C(9) 


NT01.(1) • 

TOTAL NUMBER OF DESIRED DATA POINTS ON EACH TAPE TO BE READ: 



0 READS EVTIU PILE NT0L(9) 


IW » 

INTERMEDIATE DIAGNOSTIC PRINT-OUT OF FIRST IW DATA POINTS 


HRL0W • 

LOWEST ACCEPTABLE HEATIMC RATS RATIO 


HRNICH > 

HIGHEST ACCEPTABLE HEATMC RATI RATIO 


KAS( • 

MMO m FOR M INSTRUKNTRTION STICK TNRT IS TO K DIKRMO 


WHC • 

VALUE OF h/h. AT TOP OF GRP FOR EACH STACK 


UELTX • 

icfiia SIS im) mm stack 


MM.II • 

EMC RAOIUS - LOHQI UNIT (ca) 


RADNI • 

OK NRtlW • wvn UNIT (cn) 


SEND 



OSIGWAL PAGE IS 
OF POOR QUALITY 
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MULTIPLE REGRESSION ANALYSIS 


N - SAMPLE SIZE 
Y - DEPENDENT VARIABLE 
- INDEPENDENT VARIABLE 
df » DEGREES OF FREEDOM 


0 STANDARD DEVIATION OF ALL "Y" VALUES: o-^ 

0 CORRELATION EQUATION: Y - X, ♦ Xg Cj Xj + . . . 

0 RESIDUAL • 

0 STANDARD ERROR OF ESTIMATE: S “\rS^^<^mes *^CAL^^ 

0 CORRELATION COEFFICIENT: R-^ 1 - "R^l DENOTES! 


Ri^l DENOTES GOOD FIT” 


0 TOTAL VARIANCE ABOUT MEAN: s/ ■ sj^ ♦ ^ . • • • - 

sZ 

0 F TEST: X, 

F. -2 

^unexplained 

0 TERMS ARE INCLUDED INTO CORRELATION EQUATION STARTING UITH LARGEST "F" 
0 COEFFICIENTS ‘C^”. DETERMINED BY LEAST SQUARES 


Figure 171 
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transfer coefficient ratio become identical for high enthalpy flow produced by the 
JSC 10 MW facility. The independent variables considered in the MRA included boundary 
layer parameters, gap dimensions, locations in the gap and ratios formed from these 
quantities. Loi.al flow properties (e.g. Mach number and Reynolds number) were con- 
sidered as correlation parameters because they affect embedded shock strength, flow 
expansion angle, boundary layer growth, boundary layer structure, etc. Other indepen- 
dent parameters considered cavity geometry and properties of the mass "captured" by 
the cavity relative to the structure and energy level of the recirculating flow in 
the cavity. No individual test program had sufficient variation in free stream con- 
ditions to completely evaluate all candidate parameters. Additional experimental 
data are needed to completely determine the impact of these parameters on gap heating. 

In addition to the final function (or equations) , Intermediate regression equa- 
tions are obtained after each step in the MRA, giving an indication of which vari- 
ables are most important. Also, some parameters in th'* candidate correlation func- 
ion were rejected because they had no significant effect on the dependent variable 
(beating ratio) . 

Statistic.! information is produced regarding goodness of fit, multiple correla- 
tion coefficient (R) and significance of Interaction among Independent variables. 

Of particular importance Is -he standard error of estimate (S) for each step which 
represents the MRS error of prediction (or confidence band around the regression 
line) . In following selections values of "S" and "R" are used to evaluate candidate 
function adequacy. 

5.1.2 h/h ^ Computation - This subsection contains a review of heating ratios 
available for each set of gap heating data contained within the data bank. This 
review was conducted approximately midway through the program to establish which data 
sets contained Che gap heating ratioed to the edge value on the tile top. As a re- 
sult of this review the entire Data Bank was upgraded to include h/h^. At the same 
time the effects of thermal conduction were factored into the calculation of the 
local heat transfere ' -'oefficlent (see Section 5.1.5). 

1) JSC 10 MW Arc. Tunnel Channel Nozzle Tests of Mullite RSI tiles 

These data were referenced to the calculated surface heating rate from the lone 
surface thermocouple on the test article which was normally located 0.127 cm from 
the edge of one of the tiles. 

2) JSC 10 MW Arc Tunnel Wedge Tests of Four Edge Radii Tiles 

These data were measured on thin skin tiles and were initially referenced to 
measured smooth flat plate heating at 15“ angle-of-attack. It can be shown that the 
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surface heating on the thin skin tile is significantly different from that measured 
on the smooth flat plate. The surface heating noar the gap is a strong function of 
surface location and gap width. Therefore computations were added to this data 

set . 

3) LaRC Mach 8 Variable Density Tunnel, Wall and Freestream Tests 

These data were referenced to the measured surface heating at the most forward 
location of the test article. They were assimilated during Phase I. 

Figure 55 shows data referenced to this surface heating as a function of Re^/m 
with the test article In the freestream position. It can be seen that at higher 
Reynolds numbers the flow becomes transltio.ial and turbulent. It is felt that the 
gap heating at these Reynolds numbers would correlate better if referenced to the 
local surface heating near the gap. Figure 56 shows the heating data referenced 
to the measured heating at the forward surface thermocouple (point B) when the test 
article is in the tunnel wall position. The boundary layer is turbulent over the 
entire test article. The surface heating is nearly constant over the entire surface. 
Therefore, the reference heating location is not as critical as with the freestream 
tests. Due to schedule constraints these data were not modified. 

4) Ames 3.5 Foot HWT 

Two sets of data were received from this test facility. The first set employed 

a calibration plate insert and four inserts with different gap orientation (see 

Figurel72). These data were assimilated during Phase I. A set of data for a single 

in-line gap was assimilated during Phase II. From analyses contained in previous 

sections, it appears that b/h^ would be a desirable correlation form for these data. 

The data were originally referenced to the calibration plate data. However, the 

presence of the gaps has been shown to cause earlier transition on the gap inserts 

which affects both surface and gap heating. Both sets of data were put in the form 

of h/h , 
e 

5) LaRC Mach 10 CFHT Tunnel Wall Tests 

These data were referenced to calibration plate heating distributions. The 
boundary layer is turbulent over the entire length of the test article. This re- 
sults in nearly uniform heating along the direction of flow for the calibration 
plate. The heating in the spanwise direction (normal to the flow) increases signi- 
ficantly away from the centerline of the test article. 

The data from these tests indicati that the gap heating is related to the 
heating on the top surface near the gaps. The same trends that apply to the gap 
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heating also apply to the top surface near the edges of the gap. This data appears 

to correlate better when referenced to edge value. Therefore, h/h was added to 

e 

the Data Bank. 

6) LaRC 8 Foot HTST of Gap Model 

These data were referenced to both h/h ... . , ^ and h/h . . The 

calibration plate edge 

test article was tested at 9°, 7.5°, and 15° angles to the flow and two Reynolds 
number factors. Calibration data are available for only three of the six possible 
ar.gle-of-attack/Reynolds number combinations. 

7) Ames 20 MW Turbulent Duct of Silica Tllee 

These data were taken on silica RSI tiles. The data are in the form of thermo- 
couple temperature responses. The heating rates were obtained using the Inverse 
solution technique and the data was already in the form of 

5.1.3 Boundary Layer Transition - During the ~.ourse of analyzing each set of 
data assimilated into this study, the state of the boundary layer over the test 
panel was determined. Each piece of data was marked as to its boundary layer state, 
mainly, based on results from tests on smooth calibration plates. The Ames in- 
line gap tests were handled differently because a continuous row of thermocouples 
were installed along the top of the panel near the gap and could be used to judge 
•^oundary layer transition. 

5.1.4 Specific Heat Correction of Thin Skin Data - Where necessary, corrections 

were applied to the data for differences in specific heat between the calibration 
plates and the test article. For example the heating data obtained on che top of 
the thin skin tile material used in the LaRC CFHT test were consistently higher than 
those measured on a companion calibration plate. The calibration plate was fabri- 
cated from a sheet of 304 stainless steel. Nominal values for the specific heat of 
both materials were used in the data reduction (Cp^Q^“0.12 and Cp2^j^“0.1019 Cal/gm°C). 
To verify the correctness of these values the specific heat of a coupon cut from the 
two test articles was measured. The specific heat was measured at LaRC using the 
same equipment and technique on both specimens. Figure 102 contains the measurements 
and specific heat ratio. It was necessary to extrapolate down to 300°K where the 
wind tunnel tests were conducted. The measured specific heat values used for correct- 
ing the data are Cp^Q^«0.112 and Cal/gm°C. This reducer, the heating 

ratio by 9Z. The computer program was modified to automatically correct the CFHT 
data. Figure 103 shows graphically the effects of the specific heat correction on 
tile heating. 
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A similar procedure was applied to the specific heat measurements on the thin 
skin tile and calibration plate material received from 1. Weinstein at LaRC. The 
data are presented in Figure 138. The specific heat of the thin skin material is 
3.6X higher than that of the calibration plate at room conditions (300®K). In th > 
original data reduction the same specific heat was used for both materials. Hence 
the heating ratio (h/h^p) was 5.6% lev. Results in the Data Bank were corrected 
for specific heat. 

5.1.5 Conduction Correction of Thin Skin Data - An orthogonal conduction 
correction for the thin skin tile data was inserted into uiie data processing p»'o- 
grams . 

This package arranges adjacent instrunentatlon sc that the second derivatives 
in both (x) and (y) directions can be computed. The correction on heat transfer 
coefficient is Implemented as follows: 


h 


corrected 


T -T 
AW ^W 




h , ■> h 

corrected HEAT 


STORAGE 

only- 


kz f ^ 3^T 

^^AW‘V ^ 3x^ 3y^ 


The procedure is setup so that the order of the temperature curve fit and the number 
of data points Jn the least square calculation can be specified. For example, the 
entire set of heating rates obtained on the thin skin tile tested in the HTST was 
corrected for thermal conduction in two orthogonal directions. In some Instances 
the heating rates changed by as much as 3.8%. Of course where there was little or 
no temperature differences indicated by the thermocouples, the heating rate remained 
unchanged. Figure lOI shows the change in a heating distribution on the tile when 
considering thermal conduction. 
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5.2 Correlations Formulated During Phase I - Correlations were obtained for 
both transverse and in-line gaps for the 10 MW Arc Tunnel at JSC. Transverse gap 
heating was best described by the following equation, wherein distance into the gap, 
gap width and the interactlor. between distance and depth were the most significant 
variables in the equation. 

In (r^ ) = -. 3716-5. 0249Z+2.5604Z +1.0733W-.03797ZT -.03654(Z/W) 

‘^SURFACE 

-6.0719(Z/T)^ Equation (9) 

Equation 25 was selected as the most descriptive of in-line gap heating. In this 
correlation, distance into the gap, the ratio of distance to gap width and width 
were the significant variables 


RSI GAF 


ns 


An (■ 


^GAP 

‘^SURFACE 


= -.3319 -4.3979Z+1.5630Z^-.2295(Z/W)+1.0148W 


Equation (25) 

Sets of data from the CFHT and the Mach 8 tests were submitted to the MRA program 
to determine the trends with both boundary layer and gap geometry for turbulent 
flow. The resultant correlation was 


h 

^REF 


0.01384 


2(-1.46402W“^^^) 



-2.549 


-.5362 


Re , ,,.7806 

b ' ” 

10 meter 


Equation (18) 

The independent variables have the following ranges 
0.36 ±1 5 . 74 cm 

0.127£ W < 0,711 cm 
0.29 < T,,/T <0,44 

~ ® ~ -1 6 
1.158X106 £ Re (unit Reynolds number, m ) ^19.37x10 

The CFHT data were also analyzed to determine the effects of flow angularity (y) . 
The following correlation resulted for the transverse gap rotated 90 . 

«,n z — = -0. 1054-3. 71478Z+,60432Z^+1.92179W-.00692Y^+0.04316YW+.01768YG 

hpp 

+.00235 y2z 


Equation (12A) 
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5.3 Correlations for Edge Radius - Sets of four thin skin tiles In a vedge 
test fixture were Inserted Into the flow field generated by the JSC 10 MW Arc 
Tunnel (C. D. Scott, Test Engineer). Four edge radii (0.1575, 0.3175, 0.635 and 
1.270 cm) were tested at three gap widths (0.127, 0.254 and 0.381 r- ). Originally 
the smallest radius tile was specified In the data we received to have a sharp 
edge (E = 0.0), so an entire set of correlations were developed using the Informa- 
tion. Near the end of this program. Information was received that the edge actually 
had a curvature. Unfortunately, scheduling would only allow reworking one correla- 
tion. Multiple Regression Analyses were perfonned to correlate the heating In the 
In-line gap both upstream and downstream of the transverse gap, and also on both 
faces of the transve.rse gap (see Figure 13). Analyses were performed first for 
only the vertical walls and then for the vertical wall plus the edge zone on the 
tile. Approximately 45 data points per gap location were available for the verti- 
cal walls and 63 data points when the edge zone was Included. 

5.3.1 Vertical Wall Correlations - Candidate parameters generated using gap 
width (W) , edge radius (E) and distance into the gap (Z) were used in the correla- 
tion of local heating ratioed to the heating on a plate without gaps. Both heating 
ratio and logarithm of heating ratio were tested. The logarithmic transformation 
produced the best fit. Three sets of candidate parameters (Figure Al ) were used in 
analysis. The first set which Involved (Z) was employed most extensively. The 
other sets, which contained the surface distance (S) were used to include the data 
on the radius. The third set Included Interaction between width and radius. 
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Figure 173 contains the best function (using set 1 parameters) for the vertical 
wall data at the four gap locations. Also the data for the all in-line gap were 
grouped together for an analyses. Similar analyses were performed for the transverse 
gap and a grand correlation for all vertical wall data. The regression correlation 
coefficient, standard error of estimat.*., minimum heating, maximum heat for both 
measured data and the derived functions are listed in the figure. As can be noted 
from the figure, the number of parameters was reduced from 39 to a maximum of 8 by 
the multiple regression analysis procedure The parameters are listed in each 
function in the order of significant contributions to the correlation. For example, 
Z/W was the most significant parameter in Equation 4-1. Each function exhibits a 
drastic decrease in the heating with distance (Z) into the gap. The effect of 
edge radius is coupled either with gap width or distance into the gap. 

Using the logarithmic transformation of heating ratio produced the best fit. 

A regression correlation coefficient (R) of 0.9707 with a standard error of estimate 
(S) of 0.3541 was achieved. The correlating equation for the in-line gap is: 

2 2 

In^ • 0.20113 - 0.66569 | + 0.01285 + 3.96410 - 2.94821 | 

+ 0.15665 + 0.20328 | - 0.83354 (WZ^) 

As can be seen from Figure 174 , the residuals are distributed uniformly when 
presented as a function of measured data and also as a function of calculated heat- 
ing. The fitting of the data was performed using the logarithmic value and hence, 
Figure 174 should be used to judge the degree of bias in the fit. Figure 175 pre- 
sents the residuals for the above function transformed into heating ratios. The 
correlation equation yields a heating ratio that varies from 0.007 to 0.9,11 as com- 
pared with the measured data which varied from 0.006 to 1.020. The correlation 
function has a standard error of estimate (cartesean) of 0.1190. 

To understand the goodness of fit and the data scatter, figures have been pre- 
pared which compared the correlations with the reported data. Figures 176 and 177 
contain the comparison for the in-line gap, upstream of the transverse gap. Although 
Equation 4-1 involves seven terms, the curves generated by the equation are smooth 
with Z, E, and W. Figure 176 is for the extremes of the edge radii tested. For 

246 


RSI GAF 




MCDOfV/yCLL DOUGLAS ASTftOIVAUTICS . EAST 




& 


RSI GAP HEATING ANALYSIS 

□ 


VOLUME I 






REPORT MDC E1248 
JSC 09651 


I 



s 

I 


o 

!S 

o 

s 

UJ 

UJ 

s 


l/> 

!S 


I 

X 

</> 


o 


oc 

o 

o 


o 

O- 

o: 


J 3 ; 



LO O E O 

. = 

O r- U. 

ir> 00 LU * 
• r> ro 


A O 

UJ UJ 

IS IS 


•- »— Lf> 


Q. 

UJ 

: LT U> 

j ^5 

E < 

>- u »- ^ 

■> z »^l/> “ 


» X Q 

I LU »— 

: o c3 Lhj 
: z i-i o 
< 3 < 
J h- U. 
> uO O- QC 
J — < X 
J O O i/t 

■ M N 


’t:: 5 


.2 23 


5 

; QC *-• M 


►— O 
< ^ 


ss 

O 13 
o »— 


o') O 
to 

LU O 
ad UJ 
o tO 
LU *< 
X CO 


to X lO 

UJ X UJ lO 
^ Ul 
U- U. 3 

oa 

cr 3 X > 

S o o 
^ X z 
X x< 

UJ -U lU 

S i/^ 
O UJ 
X X X 
•CC »- < X 
o < o < 

5== 5^ 

«- Z »- Z 
tO O lO o 


V3 


Figure 173 




t^CO€>9^Nei.L, OOL/^L^iS ASri?OAIAC/r#OS CO^#>>l^V • CAST 

247 






RSI GAP HEATING ANALYSIS - II 

r\ ' \ VOLUME I 


I 


REPORT MDC E1248 
JSC 09651 


RESIDUALS (.lx) FOR IN-LINE GAP 
JSC 10 MW, EDGE RADIUS TESTS 



•5 


-.It 


116 -6. in -).m -1.416 -.461 .116 

-6.146 -I.lll ' -t.m -l.«6t ' -.611' 


.tt 


0. 

CTl 


c 


0. 

crl O' 


R 


.1? 




• )A 


• 91 


• U 


ti 


1 1 1 


t 1 It 

1 t t 


I t 


It •I.IM •l.tit .tit 

•t.99t •! 991 -1.919 ••til 


.IJ1 -.ItJ* 

-t,9U -I.9^? -?.9*1 •IT919 -.Stf 


-•If 


• 9P 


*.ll 


t t 

I t 

t t 


••II 


• II «t 


I i 


1 t I 


• tt 


1 tit 

t . 1 t 


t 1 


.91 


•t.tu ‘Kill -.IIJ. 

•fc.tJl -1.919 -.919 


In 



^FP/ MEASURED 


In { ^ 

%'?} CALCULATED 


248 


Figure 174 


Mcootviymt-L. oougl/xs >«sr«POA/Aiyr#cs coiyir»Aiw • mAmr 






RSI GAP HEATING ANALYSIS 
r VOLUME I 


REPORT MDC El 248 
JSC 09651 


RESIDUALS (CARTESIAN) FOR IN-LINE GAP 
JSC 10 MW, EDGE RADIUS TESTS 


P'O'II 


.tii .111 .kit .III .111 t.KI.. 

.Ill .111 .111 .114 .111 .. 


.Ml .III .141 .141 .III .III 

•111 “Till iiii :iii 


im I I I 
ti I 
I }i 1 t 


.III .til .4lt .414 .III l.ltl. 

.lit .111 .111 .114 .111 . 


**FP / MEASURED 



.Ml .III .III .111 .III .III 

.111 .III’" .111 ;4ii " Till' 


FP/ CALCULATED 


Figure 175 


MCOONNmLL DOUGLAS ASTRONAUTICS CONtRAMV • SAST 






1 



RSI GAP HEATING ANALYSIS 

VOLUME I 


Report mdc eu48 

JSC G96S1 


CORRELATION OF IN-LINE GAP HEATING 
(JSC 10 MW TEST) 


A 



L 


1 — 1 

; — 


FLOW 







H 






ao« Minus 


fZ 


VIEM A-A 


VERTICAL WALL. ONLY 


iiiil 

DATA 

EQ. 4-1 

0.381 

0.254 

0.127 

A 

□ 

O 

— 





Z. DISTANCE FROM HEATED SURFACE (cm) 
250 


MCOON^mLL. OOC/CLAS ASTI*K»M/tUTICS COIMIf*AMV . CAST 



























1 



Q 


RSI GAP HEATING ANALYSIS 

Q 


VOLUME I 


REPORT MDC EI248 
JSC 09651 


the wide gap (0.381 cm) the function underpredicts the data near the top of the 
gap. However, while the data scatter Is considerable for the lower heating values, 
the curves pass through the data mean. Figure 177 Is for the Intermediate edge 
radii. The curves pass through the data means except for the large edge radius 
(E ■ 0.635 cm) at the narrow gap width (W 0.127 cm) where the function under- 
predlcts the data. 

Correlations were also computed for the in-line gap (upstream of the transverse) 
using surface distance (S) rather than the Z-dimension. These analyses were per- 
formed to test the hypothesis that heating is a smooth function of (S) . The 
correlation equation (4-2) for the vertical wall region of the gap is contained 
in Figure 173. Using (S) as a correlating parameter decreases the number of terms 
from eight to six at a slight reduction in the goodness of fit. 

5.3.2 Correlation for Edge and Wall - The other correlations prepared considered 
the additional data points for the curved (or edge) portion of the gaps. The best 
resulting correlation (Equation 4-9, Figure 178) has ten terms. In Figure 179 this 
correlation is compared with the measured in-line gap data. The S correlation 
appears able to accomodate heating on the tile edge as well as the vertical portion 
of the gap with no significant loss in accuracy. Figure 180 compares the best 
correlation (Equation 4-17, Figure 178) for the transverse gap with the data. Figure 
178 presents the correlation of set 3 parameters with the edge and wall data. 

5.4 Correlations for Long In-line Gap - Heating data obtained on a long in- 
line or axiai. gap was separated according to boundary layer state and correlations 
were developed for an in-line gap in laminar flow at incidence angles from zero to 
fifteen degrees . 

Two sets of laminar data were selected from the data bank for correlation; 
the first set consisted of 273 measurements at zero incidence "y", the second of 

635 points at y between zero and fifteen degrees. 

5 . 4.1 Correlations for Long In-line Gap (y°0°) - heating measurements were 
obtained for three values of gap width and Reynolds number and four gap depths. 

.\ttempts to linoarize the data proved fruitless so correlations were attempted using 
independant parameters up to the third order. Figure A1 part 2 contains the variables 
investigated. Separate analyses were developed for gap heating referenced to the 
flat calibration plate "hpp" and the edge of the gap "*^edgf" heating 

rate at y=-0.51 cm on the gap. Side B as illustrated in Figure 14. The resultant 
equations are presented in Figure 181. The independent parameters are sequenced 


252 


MCDONIVEH. DOUGLAS ASTDOMAUTiCS CO^tAAMV^ • KAST 



CORRELATION FUNCTIONS FOR GAPS WITH EDGE RADIUS 

(EDGE AND WALL) 


1 


I 






RSI GAP HEATING ANALYSIS - II 

I \ VOLUME I 


REPORT MDC EI248 
JSC 09651 




bb 


E 

:e 5 


o 

•2 


TS 

W ■' 

1 

§S: 


Ssi 

> • o ►- 
>°5 = 

' « h- O 


> r>4 fic e£ 
• O => 
O o *- 


- Ps. o 

' P>4 ^ Z 

* — W*I O 

• • irt 

> o uj a 
S UJ 
I » O 

A X Sic flo 


^ 

III 


UJ W ^ M 


£ 3 £ ^ ^ 

Ck oc < & 

U4 ^ Ui UJ 

S kJ^ »- 
O L«J 

tt 3 O: *- r 

4 < ac 2 

c5 « o < 5 


> VI O 

• p- 


253 


AtCOOWWCLi 


>o. 


ic/rfcs • mM»T 


Figure 178 
Part 1 



CORRELATION FUNCTIONS FOR GAPS WITH EDGE RADIUS 

(EDGE AND WALL) 



= ? S 

«n -i* 

« CO ^ 


REPORT mDC EU4B 
JSC 09651 



**65 


^ ^ 

— X O 

o w ^ ^ 

4 U> O u> OC 
oc ^ — u» <t 

< jK < 
w *- 

K/^ Cl oc ; 
Cl — «t ^ o 

Ui O O 


i i 

g? “ r: 

— wo r: 

►- O W-. i uJ 
< »— W 

_ 3 o ^ 

OCX wx 

<x «r -s < to O 

cxa cx ;rr» 

C 3 O O O ^ 

-U- S-" 5 5 C3 

•3 7 ^ MM 

d *r kO »- S 

— ;; o x o u 4 w X 

k/^c 5 DL 3 oe»^' CO 

i/- 4 1^ «f oc 9 O 

u*4 O n ft LJ 4 K 

cx k^4 ir ,z 7 0 < M 

L? k/> <t < ^ » 

uJ < ►- 7 »- X Q 

CX( 3 i lOOwnO 


Vj o ^ •-• 


Figure 178 
Part 2 


AsrivoMiBvrfcs coiMf^ivv • kmut 













EDG 










































REPORT MDC E1248 
JSC 09651 



if- 



RSI GAP HEATING ANALYSIS - II 

VOLUME I 

according to their individual doir.inance in the correlation. Thus, "!-nZ" was 

the most significant parameter to both the flat plate and edge referenced correlations. 

It is interesting to note that while the second order term was the next most dominant 

3 

parameter for the (h/h ) equation, (inZ) was more important in the other equation. 

r i 

5.4.2 Correlations for Long In-Line Gaps (0 ^Y£l3 )- Analysis and correlations 

were performed to determine the effect of flow orientation angle (y) on the heating 

distribution in the in-line gaps. Heating measurements were obtained for both walls 

(Sides A and B) at 0, 5, 10 and 15 degrees. When data from Side A, alone, were 

considered, there were 516 measurements. Correlations were developed for both data 

sets using (h/h ) aud (h/h„_„„) as the dependent variable. The resulting functions 
ri: hUCjb 

are listed in Figure 182. 

When the gap heating was ratioed to flat plate measurements, the effect of (y) 
was masked in the standard deviation (S) . The correlation coefficient (R) increased 
slightly (from 0.9202 to 0.9399) by excluding Side "B" data. This effect can be 
attributed to a difference in heating on the opposite walls of the gap. When the 
heating was ratioed to the edge value, the effect of (y) appears in the correlation 
(Equation 16-5 and 16-6). Distance into the (Z) , gap width (W) , momentum thickness 
(0) , and unit Teynolds number (Re/10 ) contribute more to the correlation equation 
than does *-he incidence, angle (y) . However, data from one of the low Reynolds 
number tests did show a dependence on (y) . Data for other tunnel conditions, 
widths, depths and (y) must overshadow this effect. For example, the runs at 
y = 10° show a sharp increase in heating due to gap width expecially at the 
higher Reynolds number tested. Keating on the downstream side (Side A) of the 
gap is almost twice that on the upstretm side. The data indicate thac the pre- 
secne of the in-line gap at slight incidence angles causes the boundary layer to 
become turbulent. 

5,5 Correlations for Transverse Tap -(Laminar and Transitional Flow ) - Heating 
rates measured cn the downstream wall of a transverse gap were examined for trends 
with distance into the gap (Z), gap width (W) and unit Reynolds number. Corre- 
lations were then developed. 

The largest group of data examined were obtax.'ed from the Ames 3.5 Foot H.W.T. 

A total of 668 data points were obtained using transverse gap models with gap depths 
of 1.02 cm, 2.03 cm and 4.06 cm. Heating rates measured on a smooth calibration 
plate were used as a reference and also to determine the onset of boundary layer 
transition from laminar flow. Of the 668 measurements, 316 were for laminar flow. 
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Results obtained from two test programs conducted in the NASA JSC 10 MW Arc 
Tunnel facility were also included in the analysis. Seventy eight measurements were 
used from the HCF type HRSl tiles tested in the Laminar Duct, and eleven measure- 
ments were used from a thin skin tile model tested in a wedge test fixture. These 
heating data were also referenced to heating on a smooth plate. Three HCF tile 
thicknesses (3.18 cm, 5.08 cm and 6.35 cm) were employed in the Laminar Duct tests. 

The thin skin tile tests utilized a tile 4.13 cm thick. 

In Figure 183, heating measurements from the AMES 3.5 foot H.W.T. are plotted 
two ways in an effort to obtain the simplest correlation. When plotted on the Log- 
Log scale the data shows a sharp break which is difficult to describe mathematically. 
Of note is the inconsistant trend witn unit Reynolds number. This inconsistency 
V7hether real or due to experimental technique makes the correlation development more 
difficult and contributes to the standard error of estimate. 

In Figure 184, the heating at two depths into the gap is presented in the two 
formats. Data from three gap widths (0.127, 0,254 and 0.508 cm) are consolidated on 
this figure. The trends are flatter on the semi-log scale. It should be pointed 
out that the curves in the figure are uncorrelated estimates and the presence of 
several (abnormal) points was neglected. Also, there is an apparent non-monatonic 
trend with gap width. These observations will manifest themselves in a larger 
standard error of estimate. 

Plots of the NASA-JSC 10 MW Arc Tunnel tests results are contained in other 
sections of this report. 

Correlations for the Laminar /Transitional data set were developed. First the 

entire set (.668) of data from Ames was correlated using Parameter Sets 4 and 6 (Figure 

Al, part 2). The resulting functions (Equations 15-1, and 15-2) are contained in 

Figure 185 with a correlation coefficient (R) of less than 0.89. Equation 15-1 was 

developed from the set of candidate parameters (Set 4) involving the logarithm of 

boundary layer parameters and physical dimensions of the gap. The wall to boundary 

layer temperature ratio (T /T ) contribtjted most to the correlation. Next Parameter 

w e 

bet 6 which con::ains both logarithmic and non-logarithmic terms was employed in a 
second correlation (Equation 15-2). Unit Reynolds number, boundary layer displacement 
thickness as well as physical dimensions contributed to the correlation. The AMES 
tests were conducted at four unit Reynolds numbers (10 ^xRe/m = 1.4, 2.2, 3.3 and 
4.5) and had a boundary layer with a displacement thickness of 0.21 to 0.48 cm and 
a momentum thickness of 0.127 to 0.052 cm. Additional experiments are needed to 
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verify the dependency of gap heating on boundary layer parameters. An improved fit 
was obtained using a parameter set (number 7) assembled from the trends observed in 
Figures 183 and 184, 

Set 7 

1) cubic variation with (Z) 

2) quadratic variation with (W) 

3) linear variation with (Re) 

4) linear variation with (T) 

in h/hj^p = (aQ+a^Z+a,Z^+a^Z^) (bQ+b^^W+b^W^) (C°+C^Re) (D^^+D^T) 

The parameters resulting by expanding this expression were used in the Multiple 
Regression Analysis. Only a small improvement was obtained for the 668 data points 
(see Equation 15-3). 

When Parameter Set 7 was applied to just the laminar data (316 measurements) 

from Ames, the regression coefficient increased to 0.9364 and the standard error of 

estimate decreased to 0.1439 (cartesian). Equation 15-4. 

The other laminar data was then added to the AMES data and 405 measurements 

were used to develop Equation 15-5. Distance into the gap (Z) contributed most to 

2 

the goodness of fit and (W ) the .least. The influence of (Z) is very evident and 
also enters into the most complicated terns involving (W) , (Re) and (T) . Figures 
186 and 187 are the residual plots for Equation 15-5 and show a uniform spread in- 
dicating an unbiased fit. Figure 186 is for the £n(h/h^gp) which was used to develop 
the correlation and Figure 187 presents the same data in terms of 

Figure 188 is a comparison of Equation 15-5 and a sample of measured data. Only a 
small amount of data from the AMES 3.5 ft. H.W.T. and Jf-C Arc Tunnel tests is con- 
tained on the figure. The preponderance of data is from the A.MES Tunnel and hence 
the correlation passes through the center of the data. Equation 15-3 also describes 
the data from the JSC Laminar Duct. This is not the case for the wedge tes: whicl. 
employed a thin skin metallic tile. Equation 15-5 under predicts the data from the 
thin skin tile. Figure 189 shows the effect of each gap variable on heating to the 
downstream wall of a transverse gap. 

The correlation study for the transverse gap (laminar flow) was explored 
further by adding a pseudo heating vain** of unity to each stack of instrumentation 
and rhpn. gencraLing another correlation. For example, the 405 measurements were 
obtained from separate 114 stacks of instrumentation. A stack consists of those 
thermocouples located on the side of a tile at the same X and Y coordinates. The 
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same set of candidate correlation parameters used to generate Equation 15-5 was 
used. The derived correlation function (Equation 15-6) ic contained In Figure 190. 
The correlation coefficient (F.) Improved slightly (from 0.9373 to 0.9562) and the 
standard error of estimate (S ) decreased (from 0.1451 to 0.1335). Also, a 
correlation (Equation 15-7) was developed using the gap heating ratioed to the value 
on the top of the tile near the edge of the gap. The correlation coefficient 

7 

decreases slightly and the (S ) Increased slightly. 

Several correlations were developea for heating on the downstream wall of a 
transverse gap submerged in a laminar boundary layer. Either Equation 15-5 or 
15-6 is suitable for predicting heating in transverse gap. 
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5.6 Correlationo for Transvra* Gap - Turbulent Flow - Heating measured on 
the downstream wall of the transverse gaps tested lit the LaRC « 8 Variable 
Density Tunnel (116 measurements), the LaRC Confnuous Flow Hypersonic Flow Tunnel 
'47 measurements) and In the LaRC 8 Foot High Temperature Structures Tunnel (42 
measurements) were used In the Multiple Regresjlon Analysis. Figure 191 contains 
the resulting correlation function. 
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5.7 Correlations fot Gaps with Steps - Figure 192 shows the effect of step 
height on heating in a transverse gap. The data indicates that the heating rate 
increases with both step height and gap width. Separate correlations were developed 
for gaps with steps for laminar, transitional and turbulent flow. These equations 
are applicable to the downstream wall of gaps with both rearward and forward facing 
steps and are shown in Figure 193 . The laminar flow correlation was derived using 
704 data points mostly from the Ames 3.5 ft tunnel. The transitional correlation 
was based on 461 data poin“s from Ames and 89 from LaRC Mach 8 tunnel. The turbulent 
function was based on data from three tests at LaRC. Figure 194 comperes the test 
and calculated values for the LaRC CFHT data for two gap widths. Figure 195 shows 
ta and calculated curves for the Ames laminar flow data. Figure 196 shows the 
regression residuals for the calculated heating ratios. 

The correlation coefficient is low for all three of these fits, especially for 
the laminar cases. Examination of the data seems to indicate that there is some 
interaction between gap width, step height and boundai 7 layer thickness not accounted 
for in the correlation. 
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EFFECT OF 


STEP HEIGHT ON HEATING 


IN TRANSVERSE GAP 
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0.4 
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A W = 0.46cm / 
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□ W = 0.13CM I 

<> AMES 3.5 FT HWT W = 0.127cm 

JSC CHANNEL NOZZLE (Z=0.3cm) 

▼ W = 0.72cm 
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• W = 0.2cm 
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6.0 GAP HEATING CALCULATION PROCEDURE 

A general calculation procedure has been developed using available gap heating 
correlations. This procedure Is suitable for sizing TPS, determining system per- 
formance and stro-.tural temperatures. The calculation procedure is formulated as 
a set of subroutines which can be used with other thermal model components which 
describe the thermal characteristics of the gap between tiles, the tiles themselves 
and the structure. At the present, a set of 24 correlation expressions and control 
logic have been set up in a group of subroutines so that the pat’. agt is self con- 
tained with well defined interfaces and readily identified input and output. The 
package is designed to be compatable with general heat transfer computer programs 
such as the MDAC-E HEATKAN and SINDA. The input list consists of the location 
within the gap where the convective heating is to be computed, gap geometry des- 
criptors and boundary layer descriptors. Figure 197 describes the interface of 
the main subroutine with the calling program. In addition to the FORTRAN calling 
list, a labeled common is used to store the argument list ARGL(25) and other para- 
meters passed between the subroutines within the calculation procedure. The argu- 
ment list is sized so that the subroutine can be expanded as more correlations are 
added. The key parameters with.n the argument list ARGL(25) have preset default 
values to insure proper functioning of the subroutine in case a parameter is not 
supplied by the calling program. 

The main subroutine "GAPH" is written with an option to input any parameter in 
the ARGL list when the subroutine is first accessed by using NAMELIST input system. 
Not all parameters described on the ARGL list are needed at the present time, so 
the formulation of the desired correlation function should be examined to determine 
the needed parameters. The features of the subroutine package are highlighted In 
Figure 198. The parameter "IC(jR" can be used to specify a correlation equation cr 
by setting IC(5R*=0 the logic within GAPH can L“ used to select the correlation equa- 
tion. 

In Figure 199 the correspondence between correlation sequence number, the 
equation designation, the gap configuration and flow conditions are tabulated. The 
designated corn lation equations are found earlier in this report. Several correla- 
tions were developed during Phase 1 ( 01 ) and are included in the gap heating pro- 
cedure . 

In general, the correlations were developed from data measured by instrumenta- 
tion covering only a finite zone of the gap. Hence, some correlrtions may yield 
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GAP HEATING SUBROUTINE "GAPH" 


0 FORMULATION IS SETUP FO "GAPH" CAN ACCEPT MANY CORRELATIONS 
0 SUBROUTINE WORKS WITH HEATRAN AND SINDA 

SUBROUTINE GAPH (0, NTAB, ZZ, Y, QQ) 

ARGL IS IN CABLED COMMON /COMGAP/ 

J = -1 SETS UP SUBROUTINE NAME (HEATRAN, ONLY) 

= 0 READS SUBROUTINE INPUT CARDS (HEATRAN, ONLY) 

= 1 COMPUTES GAP HEATING 

NTAB = NUMBER OF POINTS IN ZTAB (25 MAX) 

ZZ = TABLE OF GAP Z-COORDINATES (OR S-COORDINATES) WHERE HEATING RATES ARE 
TO BE COMPUTED (25 MAX) 

S = SURFACE DISTANCE FROM UPPER TANGENCY POINT INTO GAP (cm) 

Y = GAP Y-COORDINATES WHERE HEATING IS TO BE COMPUTED (cm) 

QQ = TABLE OF COMPUTED HEATING RATES CORRESPONDING TO ZZ TABLE (25 MAX) 

ARGL, ARGUMENT LIST (DIMENSIONED TO 25) 


ARGL 

(1) 


IBL 

ss 

BOUNDARY LAYER STA'‘E: 1 * LAMINAR, 2 = TRANSITIONAL, 

3 = TURBULENT 

ARGL 

(2) 


IGL0C 


GAP LOCATION: 1 = UPSTREAM SIDE OF GAP, 2 = DOWNSTREAM 

SIDE OF GAP 

3 = IN-LINE GAP, 4 = STAGNATION POINT, 

5 = TILE TOP, 6 = LONG IN-LINE GAP 

ARGL 

(3) 


IC0R 

• 

CORRELATION NUMBER TO BE USED: IF ZERO IBL AND IGL0C 

DETERMINE CORRELATION TO BE USED 

ARGL 

(4) 


JOINT 

CONFIGURATION: 1 = BUTT, 2 = COUNTUURED, 3 = OVERLAP, 

4 = INCLINED 

ARGL 

(5) 


E 

= 

EDGE RADIUS (cm) 

ARGL 

(6) 

= 

GAMMA 

= 

FLOW ORIENTATION (RADIANS) 

ARGL 

(7) 

= 

GAPW 

= 

GAP WIDTH (cm) 

ARGL 

(8) 

= 

GAPD 

= 

GAP DEPTH (cm) 

ARGL 

(9) 

= 

STEP 

=5 

STEP HEIGHT (cm) 

ARGL 

(10) 

= 

GAPFL 

= 

G;? flow length (cm) 

ARGL 

(11) 

= 

H0HL0W 


LOWER LIMIT ON HEATING RATIO 

ARGL 

(12) 

= 

H0HHI 

= 

UPPER LIMIT ON HEATING RATIO 

ARGL 

(13) 

= 

AMACH 

= 

LOCAL MACH NUMBER 

ARGL 

(14) 

= 

REPM 


REYNOLDS NUMBER/METER 

ARGL 

(15) 

= 

THETA 

=: 

MOMENTUM THICKNESS (cm) 

ARGL 

(16) 

= 

DTHK 


DISPLACEMENT THICKNESS (cm) 

ARGL 

(17) 

= 

SLYR 

= 

SUBLAYER THICKNESS (cm) 

ARGL 

(18) 


HFP 


HEATING RATE OR HEAT TRANSFER COEFFICIENT FOR LOCAL 
CONDITIONS ON A SMOOTH VEHICLE 

ARGL 

(19) 

= 

TW0TE 

= 

T^jALL/^e* temperature RATIO ACROSS BOUNDARY LAYER 

ARGL 

(20) 

= 

IH0H 

= 

HEATING RATIO: 1 = H/Hpp, 2 •■= H/H^ 

ARGL 

(21) 

To 

(25) 

= 

FOR FUTURE EXPANSION 
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ORGANIZATION OF GAP HEATING 
CALCULATION PROCEDURE 

1. SUBROUTINE FORMAT (CALL GAPH) 

1.1 USER SELECTS CORRELATION OR SUBROUTINE SELECTS CORRELATION 

1.2 BRANCH IS MADE TO CORRELATION EQUATION SUBROUTINE 

1.2.1 CONSTANTS FOR EQUATION 

1.2.2 ON SUCCESSIVE PASSES EQUATION EVALUATED 

1.3 IF DISTANCE INTO GAP IS LESS THAN "ZMIN", CONTROLLED EXTRAPOLATION IS 
USED 

1.4 IF DISTANCE INTO GAP IS GREATER THAN "ZMAX" LINEAR EXTRAPOLATION IS USED 

2. HEATING DISTRIBUTION CAN BE INPUT THROUGH NAMELIST 

3. CORRELATION EQUATION SUBROUTINES ARE NAMED USING EQUATION NUMBER 

EXAMPLE: SUBROUTINE EQ4D1 4 -^EQUATION 4-14 
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unrealistic values when extrapolated, especially very near the top of the gap or deep 
within the gap. For those functions with these tendancles, an extrapolation pro- 
cedure for the very top zone (A) of the gap and zone (B) within the gap has been 
implemented. The procedure is depicted In Figure 200. An "SMIN (or ZMIN)" dis- 
tance defining zone A is contained in the subroutine. A second order equation bridges 
the "A" zone from S (or Z) ■ 0 to SMIN. At the S-0 the second order equation passes 
through h/hj^gp - 1. For the deep zone "B", a correlation cut-off distance SMAX (or 
ZMAX) is used to anchor a linear extrapolation. The constants for the extrapolation 
are evaluated in th<^ equation subroutine and transferred back to the main subroutine 
where the actual extrapolation calculation procedure is implemented. 

The calculation procedure is listed in Appendix C. The main subroutine "GAPH" 
is structured with an input section (NAMELIST), printing cf input parameters, logic 
for selecting correlation function if none is specified and calls to correlation 
equation subroutines according to sequence number. Each equation subroutine is 
given a name derived from the "Equation Designation", for example; Equatio., 4-9 has 
a subroutine EQ4D9 which contains the formulation of constants which are evaluated 
during their first access. In turn each equation subroutine calls a function sub- 
routine which actually computes the heating ratio. All equation subroutines are 
setup with essentially the same calculation flow for easy comprehension by the reader. 
The subroutine "GAPH" computes heating ratios or actual heating rates if designated 
for the entire set of S (or Z) values input prior to printing the gap heating dis- 
tribution. If one of the gap parameters or flow field parameters -.hanges during 
the course of an analysis of a mission, "GAPH" should be re-called. 
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Figure 200 
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7.0 INFLUENCE OF GAP HEATING ON TPS SIZING 

The effect uf gap heating on thermal protection s) stem (TPS) requirements Is 
a major conclusion to be drawn from this study. The thermal protection system for 
the Shuttle consists of a strain Isolation pad (SIP) bonded to the aluminum surface 
with the RSI tiles bonded to the SIP. Aluminum structure temperature penalties are 
expected because of the possibility of having large gap widths and tiles with small 
corner edge radius. The nominal TPS configuration has tiles with edge radii of 
0.152 cm and a gap width of 0.254 cm between tiles. The greatest tile thickness 
Increase over a one dimensional model (no gap) was 42% and corresponded to tiles 
having a corner edge radius of 0.152 cm and a gap width between tiles of 0.508 cm. 
However, If the tile corner edge radius Is Increased and the gap width reduced, the 
tile thickness increase can be reduced to 12%. Heat leakage in the gap Is a complex 
combination of convectioa, radiation, conduction within the RSI tile and coating, 
and duration of heat soak. 

A set of TPS sizing calculations was made fur body point 1040 which Is located 
on the lower surface of the Shuttle fuselage where surface temperi».;ures reach 
nominally 1278®C. The effect of gap width and edge radius on TPS requirements 
were Investigated. The thermal model used in this analysis is shown on Figure 201 
and Is basically the same as described in Section 4.7 except for substitution of a 
radiation heat sink for the channel wall. The RSI used was LI-900, having a density 

3 

of 144 kgs/m and was covered with waterproof coating. Analyses were perfonned for 
the current Shuttle baseline entry trajectoiv (14414). Figure 202 Is the reference 
heating rate and Figure 203 is the local pres.sure. The reference heating was con- 
verted > ) a local heating rate using a multiplying factor of 0.3646 and was imposed 
on the top of the tile with r.ie gap heating correlation applied to the gap walls. 

“'he gap heating correlation (Equation 4-17) for a transverse gap submerged in a 
laminar boundary layer flow was used. This correlation was developed from data mea- 
sured on thin skin tiles with various edge radii tested in a wedge in the JSC 1C MW 
Arc Tunnel. The method whereby the correlation was obtained is described in Sec- 
tion 5.2. Figure 204 shows the heating distributions. As car. b? seen on this figure 
the curves do not extend deep into the gap. Figure 205 shows ;.he .nethod of extra- 
polating the correlation cu.-^v'e toward the bondline. In the gap near the upper sur- 
face, a second order curve was fit between an assigned tangency point on Equation 
4-17 and passing through ln(q/qj,p)-l .0 at the tangency point of the edge ladlus and 
tile flat surface (S - 0.0). The technique of extending the correlation curves is 
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THERMAL MODEL OF AN RSI TILE JOINT WITH EDGE RADIUS 



THERMAL MODEL NOTES: 

1. CONVECTIVE HEATING ON SURFACE AND GA.P WALL 

2. COATING THICKNESS - 0.0381 cm 

3. COATING EMISSIVITY = 0.80 

4. ALUMINUM STRUCTURE = 0.203 cm 

5. SIP + RTV ADHESIVE = 0.445 cm 
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REFERENCE ENTRY HEATING RATE 

TRAJECTORY 14414 



LOCAL HEATING = (a^/^REF^ 

FOR BODY POINT 1040 = * 0.3646 
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Figure 202 
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LOCAL PRESSURE ON LOWER SURFACE 
OF SHUTTLE DURING ENTRY 



TRAJECTORY 14414 
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LOCAL PRESSURE = (1/3) TOTAL PRESSURE 
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TECHNIQUE FOR EXTRAPOLATING 
CORRELATION FUNCTIONS 


1.0 


O" 

CT 


c 



ZONE A 


INSTRUMENTATION 

ZONE 


ZONE A FIT WITH 2ND ORDER 
FUNCTION, TANGENT AT A AND 
PASSES THROUGH 1. 

ZONE B FIT 
SLOPE OF 1 CYCLE 
PER 0.5 cm FROM Z. 


TO BOTTOM OF GAP 


"MAX 


ZONE K 


ZMIN ZMAX 

Z OR S (DISTANCE INTO THE* GAP) 
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also describee', in Section 6.0. The local pressure was used to obtain the proper 
thermal conductivity of LI-900 which Is a function of pressure and temperature. 

Complete trajectory analyses were performed by computing transient temperatures 
using the 2-Dlmenslonal thermal model until the aluminum structure (0.203 cm) reaches 
maximum temperature. Solutions were obtained for gap widths ranging from 0.0 to 
0.508 cm and Inner edge radius (E^) ranging from 0.1524 cm to 0.635 cm. The sensi- 
tivity of aluminum temperature to tile thickness is shown in Figure 206 for the 
nominal gap width of 0.254 cm and > 0.152 cm. By comparison, a one dimensional 
thermal analysis Indicates 7.315 cm of RSI would suffice with no gap present. In- 
formation from the previous figure Is also presented In Figure 207 as the ratio of 
TPS thickness with a gap to thickness without. Approximately 33% more RSI Is needed 
considering a nominal gap. 

Gap width between the RSI tiles has a significant effect on TPS requirements. 
Figure 208 shows that Increasing the gap width from 0.254 cm to 0.508 cm results In 
a tile thickness Increase of seven percent for a tile with an edge radius of 0.152 cm. 
Figure 209 shows the same Information relative to the 1-D thermal model requirements . 
Likewise, the tile edge radius has a significant effect on TPS requirements as can 
been seen on Figure 210. Increasing the tile edge radius (not Including coating 
thickness) from 0.152 cm to 0.305 cm reduces the required RSI thickness by five per- 
cent for tiles with a gap width of 0.254 cm. Figure 211 shows this Information nor- 
malized to the 1-D thermal model requirements . 

Typical temperature-time history plots (Figures 212, 213, 214) show the response 
of the tile, the coating on the tile top and gap wall, and In the center of the tile. 
At the time of peak heating, the In-depth temperatures of the coating In the gap and 
temperatures in the RSI adjacent to the gap coating (not plotted) are hotter than the 
in-depth temperature of the RSI toward the center of the tile. This shows the effect 
of gap heating on temperature. 

Figure 215 summarizes the amount and percent change of RSI (LI-900) required to 
limit the aluminum structure temperature to 177°C for several combinations of gap 
widths and edge radii. The percent change over a 1-D model (no gap) ranges between 
12 and 42%. 
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ALUMINUM STRUCTURE MAXIMUM TEMPERATURE 
VS. RSI TILE THICKNESS 
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GAP WIDTH EFFECT ON TPS REQUIREMENTS 
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Figure 209 
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SURFACE, EDGE RADIUS AND 
ALUMINUM STRUCTURE TEMPERATURES 
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ANALYSIS 

EDGE 

RADIUS 

GAP 

WIDTH 

X- 

RSI 

X2-D/X1.D 

NUMBER 

cm, (inch) 

cm. 

(inch) 

cm. 

(inch) 

1 

1-0 MODEL 

NO GAP 


7.315 

(2.880) 

- 

2 

0.152 

(0.060) 

0.127 

(0.050) 

9.119 

(3.590) 

1.25 

3 

0.305 

(0.120) 

0.127 

(0.050) 

8. 738 

(3.440) 

1.19 

4 

0.635 

(0.250) 

0.127 

(0.050) 

8.204 

(3.230) 

1.12 

5 

0.152 

(0.060) 

0.254 

(0.100) 

9.703 

(2.820) 

1.33 

6 

0.305 

(0.120) 

0.254 

(0.100) 

9.271 

(3.650) 

1.27 

7 

0.635 

(0.250) 

0.254 

(0.100) 

8.687 

(3.420) 

1.19 

8 

0.152 

(0.060) 

0.508 

(0.200) 

10.414 

(4.100) 

1.42 

9 

0.305 

(0.120) 

0.508 

(0.200) 

9.957 

(3.920) 

1.36 

10 

0.635 

(0.250) 

0.508 

(0.200) 

9.474 

(3.730) 

1 .30 


THE ABOVE RESULTS WERE OBTAINED USING EQUATION 4-17. 


300 


Figure 215 


/MCOO/V<VCC.I. OOUCt-AS ASmOMAUTICS COI¥ir*ANV • AAST 









RSI GAP HEATING ANALYSIS - II 

m VOLUME I 


REPORT MOC EI248 
JSC 096 5 1 


8.0 CONCLUSIONS 

Convective heating on the surface and within the gaps between representative 
Space Shuttle RSI tile configurations has been extensively Investigated. Data 
from 10 test programs In 6 NASA facilities provided the experimental base for 
this Investigation, which was conducted In two phases. Major conclusions from this 
Investigation are summarized below. 

A number of major conclusions from Phase I were reinforced during Phase II 
and none was Invalidated. Some of these are: 

o Gap heating Is a three dimensional phenomenon. Both transverse gaps and 
In-line gaps (as well as gaps of Intermediate sweep angles) experienced 
significant lengthwise and wall-to-wall variations In heating In addition 
to depthwlse variation. 

o Interactions exist between gap heating and heating on the top surfaces of 
tiles. 

o Gap heating Increases with width; however, for a given width and distance 
from the tile surface, heating decreases as the gap depth Is Increased. 

o Tiles having forward facing steps experience higher gap heating than tiles 
with rear-facing steps or no steps. 

o Gap cross section geometry Is significant for wide gaps (greater than 0.25 cm) 
and less Important for narrow gaps. 

o With a turbulent external boundary layer, gap heating at flow Incidence 
angles between 0 and 7t/2 Is higher than at either 0 or 7t/2. 

Additional significant conclusions which were drawn from the Phase II studies 
Include: 


1 


i 



o In-line gaps promote boundary layer transition. Examination of data from 
the Ames 3.5 Foot HWT Indicates transition is Initiated within the gap and 
propogates outward. 

o Gap heating appears to Increase more rapidly than tile surface heating as 
Reynolds number is Increased. (A heating increase phenomenon similar to an 
"energy dump" was observed on a tile trailing edge tested in a turbulent 
boundary layer at the LaRC 8 Foot HTST.) 

o Increasing edge radius Is observed to Increase gap convective heating slightly. 
However, lower tile temperatures result, sli.ce the heating Increase Is more 
than offset by the greater radiant heat rejection afforded by the Increased 
radius . 
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Data correlation activity during the study yielded 21 equations describing gap 
heating distributions as functions of boundary layer state and gap geometry. Several 
of these were originally developed during Phase I and improved by utilizing the larger 
data base made available during Phase II. Among the more important correlations are 
the following: 

o Transverse gap correlations which have been improved through the incorpora- 
tion of additional data 

o Correlations for edge radius 

o Correlations for long in-line gaps including sweep angle effects up to 
15 degrees 

o Correlations for gaps with steps 

A new gap heating subroutine called GAPH was prepared during Phase II. The 
subroutine is designed to be compatible with thermal analysis programs such as SIWDA 
and HEATRAN. The subroutine incorporates che correlation equations developed during 
the study together with internal logic to select the appropriate correlation equa- 
tion. The desired correlation can also be selected directly by the program user. 

Parametric evaluation of gap width and radius effect shows increase in TPS 
thickness caused by presence of gaps varies between 12% and 42%. Heat rejection by 
radiation from gap regions is enhanced by increasing tile edge radius. For instance, 
increasing radius from 0.152 cm to 0.305 cm reduces RSI thickness requirement by 5%. 

Finally, though not a study objective, it is worth noting that during the data 
analysis a persistent tendency existed in all the thin skin tile models for the 
apparent heat transfer rate on the upper surface to decrease near the tile edge. 

Tests with filled gaps seem to indicate this does not result from a fluid dynamic 
phenomenon produced by the gaps. Attempts to explain the mystery on the basis of 
lateral conduction in the thin skin indicate the magnitude of skin conduction is 
not sufficiently large to account for the observed effect. The phenomenon remains 
unresolved . 
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S.O RECOMMENDATIONS 

Gap heating data from ten test programs In six test facilities provided the 
' asls for this study. The analysis of these data revealed a number of surprises 
for which satisfactory explanations have not been found. Not all of the unexplained 
phenomena are of sufficient Importance to the Space Shuttle Program or possess suf- 
ficient Intrinsic merit to deserve further attention. In addition to the unexplained 
T)henomena, the present body of data contains omissions which deserve attention. The 
ecommendatlons which follow are Intended to address significant phenomena which are 
lOt understood and to address significant voids In the existing data. 

As stated In the Conclusions (Section 8) , gap heating Is observed to be a three 
dimensional phenomenon. Significant lengthwise variations In heating rates were ob- 
1 erved for all gap orientations In the absence of nominal pressure gradients on tile 
surfaces. In addition to the observed three-dimensional Ity of heating within gaps, 
results from several of the tests show evidence of strong interactions between heat- 
ing within the gap and adjacent tile surfaces. (Similar results have been observed 
using teflon models In the Ames facility. Reference 5.) The three dimensional nature 
of gap heating and its interaction with tile surface heating merits further explora- 
tion. Such exploration should employ pressure measurements as well as oil flow and 
thermal map’ x,ig techniques. One or more of the existing t«#'t panels could well be 
iiiodlfied to perform such an Investigation. Also of merit is the development of an 
analytical computational method for computing flow fields with the gaps between 
tiles which include viscous effects, real gas effects, pressure gradients, flow field 
sunks and sources, and heat transfer. 

This currenL program Investigated In-line and staggered butt type gaps at flow 
orientations . anglng from 0 to tt/2 radians In the presence of a turbulent external 
boundary lyer (LaRC Mach 10 ClUT Wall Test) . Peak heating over large portions of 
the Shuttle will occur In the presence of a laminar boundary layer at flow incllna- 
ti ns of approximately tt/ 4. The tests previously conducted in the LaRC CFHT should 
ue repeated with a laminar external boundary layer, using. If possible, the same 
test apparatus. 

The following activities are recommended to be pursued, either individually 
or (preferably) in conjunction with the above: 

o Expand correlation activity to Include T slots and the most forward 

corner on tiles oriented at flow Inclination angles in the vicinity of r/4, 
o Investigate the tendency of heat transfer rates measured using thin skin 
models to decrease near the tile edge. 
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Investigate differences between results obtained fron thin skin end RSI 
models by building duplicate models and testing In the same facility. (The 
need for one of the duplicate models could be satisfied by using results 
from the current program.) 

Perform gap heating tests using companion RSI tile models each with a 
different edge radius to further define the benefits (or disadvantages) 
of radlused tiles especially when exposed to the high enthalpy flow expected 
during flight. 
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APPENDIX A 

LIST OF CORRELATION PARAMETERS 


Figure A-I presents the lists of the Independent correlation parameters sub- 
mitted to the MRA program. These parameters are divided into seven separate sets 
so that a particular set could be used in the analysis of a specific gap heating 
effect. 



McoowwmaL oouaemm Asrarowmi/rfcs coMnewr - ecmmr 



a 


RSI GAP HEATING ANALYSIS > II 

J^! ,, .1 VOLUME I 


SET 1 
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PARAMETER SETS INVESTIGATED FOR 
HEATING CORRELATIONS 


Z, 

z.^ 

Z,^ Jz^ 

tnZ, 



E, 

E.^ 

E.^ yr. 

llnEt 



w. 


w,^ /v. 

i,nW, 

„1.5 


ZE, 

zw. 

Z(E^) 

z(vh. 

Z/E, E/Z, (W/Z), 

(W/Z) 


EW, 

E(Z^) 

E(W^), 

Z/W, E/V. (W/E), 

(W/E) 



W(Z^), 

W(E^), 

(Z/E),^ (E/Z)^ 






(Z/W)^, (E/W)^ 



SET 2 


SET 3 


s. S, 

^ s.^ 

i-nS, 

(JlnS)^, Wins, W(£nS)^ 

E, E. 

2 

UnE, 

E£nS, E(£nS)^ 

w. 


£nW, 

E^(inS)^ 



WE, 

WE£nS, WE(£nS)^ 



W(E^), 

W(E^)£nS,W(E^)(£nS)^ 

s/w, 

E/W, 

W(S^). 

(S/W)^, (E/W)^, E/S, (E/S)^ 

S 1 S 1 

^ s,^ 

inS, 

(£nS)^, Wins, W(lnS)^ 

E, E, 

^ 2E, 

inE, 

EtnS, E(lnS)^ 

W, 2U 


InW, 

E^(lnS)^ 



WE, 

WElnS, WE(lnS)^ 



W(E^), 

W(£^)lnS,W(E^)(lnS)^ 

S/W, 

E/W, 

W(S^), 

^S/W)^, (E/W)^, E/S, (E/S)^ 

(W+E), 

(W+E)^, 

£n(w+E), to (W+E) ^ 

(W+2E) 

. (W+2E)^, 

ln(W+2E), ln(W+2E)^ 


NOTE: SURFACE DISTANCES ARE ASSIGNED A POSITIVE SIQ( FOR PURPOSES OF THE 

CORRELATION 

Figure A1 
Part 1 
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SET 4 


InZ, (InZ)"', (InZ)^ 
UnL , ii.nRe 

JlnO, iln 8,*, Hn T /T 
’ W' E 

. '§* Z L i* Z 

«,n ^ ^ ’ r" L W 


0 i* iP 

L W 


A.n y ^,n y 

«,n W, ilnUp, (UnZ) (UnW), (£nZ) (ilnW)^ 
KlnZp («,nW)^, (£nZ)^ (£nW) 


SET 5 


Same as SET 4 with the addition of y, y . (COSy) 


where: Z - Distance into the gap (cm) 

L ■ Gap length (cm) 

Re “ Local unit Reynolds number (lO^m ) 

G - Momentum thickness (cm) 

6* “ Displacement thickness (cm) 

W “ Gap width (cm) 

T - Gap depth (tile thickness), (cm) 

y ^ Flow incidence angle (radians) 

— “ Temperature rates across boundary layer (T = wall temperature) 
E ^ 


Same es SET 4 with the addition of Z, Z , W. L, Re, 0, 6 *, T^^/T^ 


1) cubic variation with (Z) 

2) quadratic variation with (W) 

3) linear variation with (Re) 

4) linear variation with (T) 

iln h/hj^p- (a^+a^Z+a 2 Z^+a 3 Z^)(b^+b^W+b 2 W^)(C^+CjRe)(D^+D^T) 


A3 
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Figure A1 
Part 2 
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APPENDIX B 


ooco 


aocist 


0000:3 

mv:i 

ii 

CdUll 


062 a 


QCw:2% 

006030 


C800^3 

CQC06? 

:06iio 


00611^ 


000110 

QOQllT 


TABULATION OF CORRELATION EQUATION 15-5 [TRANSVERSE GAP, DOWNSTREAM WALL, 

laminar flow] 

2il5 t'. IMm, i.m./ 




DATA T/ Ij, 2,, i., 
DATA AE/ i.Ol. 0.1. 


1 .. 1 ./ 


it.a 

j*>h 

h?th8c»i 


MOH 

eO(2>«lUHtMNS 6*r 

me THfCKNCss 


00 

w«; 

MP] 


.oil EQ 


REYAOIOS AUHOEN 

30. to: jj^^jp.2 


H*lTEt6,«CM 
MRITEIt,9:2t T(I), Tin 
MR]TE(i,4i3) REIJI. *114*11 
HP|TEI6;»C?) 

TE(6 .4011 
TEI0.4C‘ 

T|lI,9C 
Tf I6,4i 


HP 

HP 

MR 

HR 


.4001 

,40 

,4iT> 


«MILL> ,Uh 1,LH) , INItLI ,Up1,LMI 


2- distance 

00 76C K>1,K2 

GAP AICTN 
00 2S0 L'l.LN 

JJ* J-1 

00 TkC l«*l,2 
L»L*S* 


8b*Vo^«*««i8U»,HETpao 

EQUATION is*. 

tgi 0* ULI*-0.A5S7T-6 .*CW}pE . 

! -C.11073*I2«*I*»JI •Nlil*»tl 

: :n?j 

4 •(.]MS)*2UI**C(JJ» 


•*u> 

:(4Ji 


01)0173 

I0C17S 

lOCtTt 




GO TO 740 

104 continue 


0002:3 

000233 

000243 
00. Hi 
S00270 


100273 





710 continue 

713 CONTJnuI 
•00 CCNTInuC 


igtILL) ,Ll*t,L«l,2UI, !Q*aL**l,LL*l,.MI 
IQR(LU,LL*1,LM>,2IKI, IQ*IU*.I.LL*i,N.Hl 


STCP 1 


4.3 PORPATIlHlI 
401 roRPATijol* 
roRPATl|l|< 


lUN * 2 Aia) 

... . ... DrPTH IT, CN 

4.3 EO»PAtlil<«*UNTT PET NO, I*? N- 


404 roRPATI* IlOP S(P,J(* 

4fS FORPAtlOl* ICN) »Sr7.3, 

401 FORPlTl2l|Iir7.S,2'll) 

4,7 FORPATIIHI 

4{{ FOfPltllM*,SCI,* 2IM $)•,.(* 
401 Fo*pATi2«f«pr.l,ix.irr.4.i«ii 


H 

2«U 


If 




•11 


ENO 




peecedwo paob 

B1 
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tOUki lON T«MS C4P 


CiP CEPtH 
UNIT iCV NO 


It, CH 

l»f N-l, 



C*P OEPTM 
UNIT ICT NO. 




‘‘.m 


no* s> 

IbNt 


.86$ 


M 

.'.17 


.2«6 


M M 

.m .808 



ENt 

.88 


.86$ 


.117 


.781. 


.181 


.888 


:\n 

’M 

;2?* 

'.88C 

.830 

.780 

.808 


1. 

1.20 

1.8 

2.28 



CCUl.ION 18-8.(.8M.T*8NS 88P 


r,«f DtPtM 

UNIT l(r NO, 


7 IDP 8 t 
ICH 

. 8 ;i 

.ots 

. C83 

. >98 

8 J«C 

. SCO 

. <« o : 

• s.c 

8^00 

;.'!!? 

• TCC 

uccc 

1. ’CC 

i. 

1.^:. 

li ^do 


8 

8 r«9 

• Hi A 

8 : 9^ 

8 ' lb 

8l t79 

.11 A 

8-bb« 
.0 SAf 

8 0 ? J : 

. K9l 
8 i.'iq 
.lyAC 
. ]C?C 
.)wll 
. JJM ^ 

* r. « A 

.nOb 


r.AP oCpth 
UMT A£t HC. 


(To Cn |i 


M 

.1?^ 


2.o:c 

• blQ 


.9)9 
, 9*5^ 

. b^b 

• kfi 

• 

8 ^bS 
8 lAb 
.IH 

,:hst 

8 

.Cl^b 
8 ..Ub 
8 : Jb*} 
, 3 Jbf 

liCCb 


<»? **-i, 


N 

« 

w 

288 

. 18 1 

.s;0 

ur 

1.2b*i 

1.282 

Off 

I.itz 

l.l^b 

ACf 

• Abb 

.Alb 

881 

.bbb 

.bbt 

bj? 

.b^C 
8 JbS 

.bfl 

. isf 

jrn 

??; 

.288 


.Ibi 

.lib 

O'^f 3 



1AM, 

• Cbbb 

C 3b 1 

.03bb 

,9szn 

3213 

.0288 

.0200 

:i3fl 

• Clb6 

8 : 111 


• dllA 

.ou« 

0 :t>3 

• OOAl 

• OObS 

• c les 

. Q0«1 


:V.t\ 

o:iJ 

• q:zz 

.0w3A 

\iV, 

:oi^l 



!■ 

2.099 

SC-SI* 

1.0..0 

H 

W 

M 

?Sb 

.301 

.SOI 



:i\n 

:S!I! 

• cosc 

• fibt 

• o:)b 

• 6^37 

• iobt 


GiP OC^TH 
UNIT ICY NO 




!sa. 


SI 


• o:i 

I CIS 

• bl»3 

• IOC 

• ISO 
I3J0 

• bOO 

• S34 

m 

i8i: 

It (04 


1.880 .ecc7 

urn lira 

88. OC.TH 
UNIT 8C1 NO 




'Ml \:\ll l.m l:(U 

.828 .788 .I«8 .l»t 

.*28 .82° .88$ .888 

.287 .]|r .838 .898 

:m :lh :J« :ifi 

'.8«l8 


*408 e 

o* *008 ooiilnj 
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t 

p 

1 




COUA" 

ION IS 

-5, LAN. 

GAP 

UNIT 

ccpth 

AFT NO, 

Ilf s: 

1 ■ 

It 10-AI« 

S.QOO 

.010 

1 

7COP SI 

N 

M 

H 

W 

w 

1 

ICil 

9 JUS 

»izz 

.2SA 

.381 

• S3A 

1 

. c :i 

,rAQ 

9 9 j9 


1.2hS 

1.2A2 

9 

. c u 

9F%e 


l.OSA 

1.172 

i.lA9 

1 

9 

9 bF 

.hT<9 

.8"jr 

.89A 

• 89A 

1 

9 ICO 
9 15C 

9*»(,A 

.ATI 

• 5AJ 

• 6Af» 

• 6AA 

9?*»^ 

9 ui 

.A22 

.A70 

.A71 

i 

.?3C 

.:i6 

.ZA9 

. 109 

.3AS 

• 1A7 



9 ISA 

:f?o" 

• 2SS 

• 2‘8 


9 ^3l 

911A 

.lU 

.191 

• 19A 


9 ^C§ 

9jbt' 

• GTb9 

,09bZ 

.1093 

• 112F 


9O JAA 

.OAAT 

.OSAA 

• Of*A9 

.0(01 



9^:21 


• 3?A1 
9 0?1 3 

• 0399 

• 0A28 


9 rc: 

9 */ 1 >•( 

9 :if 6 

.OT'Sa 

.0280 


9 "Ct 

92CSC 

.31CA 

9 3 1 JN 

1 c. 8 

.Cl9t 
• 013* 


.o:c 

.CuSA 

9 jqts 

• 0392 

• Olis 


1. ccc 

• cLliC 

. jOat 

9 :ef>3 

.0081 

.Cy99 


i.rcc 

9 J u?C 


.0033 

• 0<^AS 

. OOAl 


KPT. 

:ccJ^ 

• lsIa 

• y 2 2 

*U329 

. 0 i A i 


1. bJC 

9 :jc9 

.0313 

• 00 22 

. 08 3A 


Z* 0 


• wOCS 

9 300 9 

• OOli 

:VM 


>.Z3C 

9WCA 

. JCtS 

.3009 

•uolo 


GAP 

OCPTH 

IT. CH 

l« 

3.030 


UNIT 

AEV NO. 

Cr£ m- 

1. lo-si* 

1.000 


Z COP SI 

M 

M 

N 

M 

u 


ICN) 

.ubS 

9 1Z7 

9?9A 

• 381 

• SJO 


9 o': 

9 7A9 

.909 

}-5n 

* 9 jbc 

1.2AS 

1.2A2 


9^‘». 

• AS2 

i.lF»7 

1.16S 

.879 


.S*>6 

.A4.2 

9 

• 881 


GIF OE8TH ITi CN >• 

UNIT llV NO. INf N>1, 10-6I* 

Z(0i S> N N M 

ICNI .06S .IZT .2S0 

;8!i ‘M\ ‘.til \:in i 

.0«0 .96T .»f3 .tuk 

.'00 •‘•gr .»6s .««o 

iip 

.2!9 .iSe .103 .227 

IlOO .117 ..3S .160 

\ni :ms ‘.nil 

,60C .C229 .0266 .0330 

{no .01*1 .,16S .0212 

(Oog .oojo .0105 .Ji:‘ 

MOO . 0059 .Il0c4 .0051 

ilcgg .go*g .00*7 .gg63 

i;200 .0020 .302* .0033 

i)*og .0011 .oci* ..o{o 

lUoO . 0007 . 0009 .30)^ 

lim 


2*5 1.2*2 

17t 1.1 * 


COP OCPTH «t, cm 
UNIT NET NO. INC M- 


2109 SI 
ICNI 


It 10-61* 




• 8SJ .Us 

•^05 .$65 

• «lO0 

:IS? :m 

.151 .152 

.0027 .0**9 

Ml 

.Cld2 .31^3 


eos4 .0:8% 

0:3<i .9036 


•0010 •GOOt 
•0009 .OOO* 


r.IP C£PTH 
U*»IT tf» NO. 


? COP 
ICii 


EQUATION l9-$. LAM. TRANS 6AA 


"-1, U-hl* 


CAP OCPTN (T. CN 
uRxr lev no. ii! m- 


eioi SI 

ICNI 


I 9 


f.AP CCPIM 
U»>lf ait kc. 


CPE N-l# K-fcl» 


K W 

.?*«. .1*1 


liOOtf •OOOA .0003 NiOfl 
e*200 •OOOA .lOOf .0019 

GAP OePTM ITj CM I 

UNIT ACT NO. IRC M-1» 19-A) 


ZIOI SI 
fCMI 


i!ii m 

.130 •ifi 


#303 .C9A 

•400 .wlrO 
•roo .Oiii 


itli i:ii| ‘ 

.*26 . 52 ! 

.297 .3n» 


120 .Oli* 



liiu 

• iP9 

• 909 

• <.C9 

. 0**2 

. 0290 
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APPENDIX C 

GAP HEATING CALCULATION PROCEDURE LISTING 


A complete discussion of the calculation procedures is 
contained in Section 6.0. The following 23 pages contain 
the detail listing of the gap heating calculation procedure. 
It is written in Fortran IV and consists of a labeled common 
(COMGAP) , a main (QRATIO) and 36 subroutines. A list of 
program components is shown at the right. QRATIO is a 
demonstration main program and should be re-written to 
interface with the other parts of a gap thermal model. A 
sample output from the procedure is shown at the bottom of 
this page. 


l-mviiOlNG PACE blank NOT HUtlsa 


/CONGAP/ 

QRATIO 

GAPH 

CSCOTT 

EQ41716 

ORF 

TBLU 

EQ4I39 

SHAXF 

SMINF 

ORF409 

CQ4010 

QRF4O10 

QRF4011 

EQ4011 

EQ4017 

0^4013 

SJfSJia 

IrfS^is 

mn\ 

EQ9 

I3R>'9 

EQ18 


.05 


>01 


EQl( 

QRF: 

EQK 
QRF1602 
EQ13D34 
QR16034 
EQ16056 
QR16056 
EQ17123 
QRl 7123 


SAMPLE OUTPUl 


6fSn8H"tlisWluti8l5'''rfi6 


I 
f 


6 

7 

8 
9 

1C 

11 

12 

13 

14 


OC 




.40000 

.50000 

.60000 

.70003 

.80000 

.90000 

,.00000 


.?0576 

:nm 

.0p31 

.00061 

0.00000 

0.00000 

0.00000 


O.OOOOQ 

0*00000 

0.00000 

0.00000 




.254 CM 


.01331 

•00061 

1*00000 


o.Qooeo 
0.00000 
o.ooor 


0*00 
0.00 
0.0009 
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PROGRAM QRATIO(INPUT=10B, OUTPUT, TAPESsI NPUT.TAPE6=OUTPUT) 

GA'' HEATING FACTORS GENERATED BY CORRELATION EQ. 

MARCH J.975 H.E. CHRISTENSEN 

TEST MAIN 

COMMON /COMGAP/ AR6L(25I, ZHIN,FMIN,DFOZA , ZMAX, FMAX-OFOZB 

DIMENSION Z2<25l. Q0I25I 

DATA iaL,IGLOC, IGAPC.IHOH/I.a. 1,1/ 

DATA ZZ /O.O. 0. 05, d.iZ6,0,10«e.l5,Q.20. 0.30*0. 4,0 .5*0.6, 

1 0*7, U.8, 0.9, 1*0, 1.2, 1.4, 1.6, 1.8, 2. 0.2*5, 3. 0.3. 5. 4.0/ 

DATA NTA8/23/ 

DATA ICCR.H, T, REPM /12, 0*254, 4*0, l.OE+06/ 

EQUIV'LENCE <ARGL ( 3 ) ,I C CR I , ( ARGL ( 7» , W f , ( AR6L(8),T», 

1 (ARGL(14).R£PM), ( ARGL(1),IBL). 

2 ( ARGL (2 ) , TGL OC I , ( ARGL (4 ) , IGAPC ) 

ITRIG=1 

ITRIG = -1, ALL PARAMETERS NECSSARY TO SELECT ANO EVALUATE 
GAP HEATING EOUATIOA ARE CONTAINED IN (ARGL) 

ITRIG = 1. NAMELIST INPUT OF ADDITIONAL (ARGL) VALUES. 


ITRIG = 


IS SET TO (2) BY SUBROUTINES AFTER FIRST PASS* 


DEFAULT— TRANSVERSE GAP (EQ.9) 
CALL GAPKITRIG,NTA9,ZZ,Y,QQI 

STOP 1 
END 


HCF/JSC TESTS (EQ.9) ICOR =12 
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c 

c 

c 

c 

c 


c 

c 


c 

c 

c 


c 

c 


c 

c 


SUBROUaNE GAPHC ITRIG, NTAB, ZZ, Y , QQ » 

G.P HE, TING CAIC^ATION METHOD CHRISTENSEN 75 

COMMON /COMGAP/ ARGLC25), SMIN.FMIN, OFDSA,SMAX, FMAX» OFDSB 

□ IMENSICN ZZ(25)t QQC25I* IRGL(25), ZTAB J25 » »QTA8(25) » 

1 SS(25) 

EQUIWfl-ENCE CARGL.IRGL»t (ZtS) 

EQUIVALENCE * ARGL ( 1 ) .IBU . (ARGL (2) . IGLOO, {ARGLCO, ICOR), 

1 CARGL 14 », I6APC» , CAROL 151 »E)f^ , 

2 (ARGL (7 ) vGAPH ) , ( AR CL (8 ) ,GAP0 ), (ARGL(9), STEP ), 

3 IARGUIOUGAPFlI, <ARGL<iI»,HOHl6mI, URGLCpj ,HOHhI» , 

4 CARGLC13),AMACHI, <ARGL (14) ,REPM) . ** Sgii (If i J ^HET A I , 

5 (ARGL IlGlJoTHK ), (ARGL (17) ,SLYR ) , C ARGL ( 18 ) ,HFP ). 

6 (ARGL(19) iTHOTE), (ARGL(20) ,IHOH ) 

DATA SS/-9999./ 

DATA EjGAMHA^GApJl^GAPbt^STEP^^GAPFL, HOHLOH, HOHHI.AMACM, 

0. 00 dol 0*^2 fOlz^L^sIOfoIoClO, 0*0.0801,3.0 ,5.0,1.QE*06,0.1, 

0.0 01 , 1 . 0 , 1 . 0 / 

AVERAGED GAP HEATING FROM AMES 3.5 FT H.H.T. 

DATA ITQ /12/ 

DATA HaB^/ 0.0, 0.01, 0.015, 0.023, 0.036, 0.060, 0.10, 0.14, 
^OATA^5TAB^/I.6,^i! 2^*1^ G, 0.8, 0.6, 0.4, 0.2, 0.1, 0.05, 0.025, 


1 

2 

3 


1 0 . 0 , 

DATA 


0 . 0 / 

ID /12/ 


NAMELIST/LISTI/ ARGL,IRGL,IBL,IGL0C,IC0R,£, GAMMA.GAPW.GAPD, 

1 STEP, GAPFL, HOHLOH, HOHHI, AMACH, REPM,OtHK, SL >R ,HFP, 

2 tHOT£,IHOH,2TAB,QTAB, ITQ, NTAB, ZZ,Y,SS 


IF (ITRIG.EQ.2) GO TO 10 

IFdTRIG.EQ. -1) READ (5,LIST1) 
HRITE(6,9Q0) 

IF (SS(1) .LT.-9990. ) GO TO 4 
DO 2 I=1,NTA9 
ZZ(I)= SS<I) 

CONTINUE 
00 5 1 = 1, ITQ 
SS(II = ZZ (I) 

MRITE(6,LIST1) 


C3 
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c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 


c 

c 


VOLUME I 

EQUATION DECISION 

10 IFdCOR .GT. d) GO TO 50 

IFtlBL.NE.l .OR. E.EQ. C. 01 GO TO 12 

CORRELATIONS FROM EDGE RAOlUSEO MODELS 
IF (IGL0C.EQ.2 ' ' 

IFIICOR.NE.IO 
IFdGLOC.EQ.2 
IFdCOR. NE. ID 
IFdGLOC.EQ.3 
IFdCOR. NE. ID 
IFdGLOC.EQ.3 
IFdCOR. NE. ID 
IFCIGLOC.EQ.l 
IFdCOR. NE. ID 
rF(IGLOC.EQ.3 
IF (IGOR. NE. ID 
rF(lGL0C.E0.3 
IFdCORcNE.IO 
IFdGLOC.EQ.l 

CORRELATION FOR TRANSVERSE GAP, COVN-STREAM WALL, LAMINAR (15-5) 
12 IFdBL.EQ.l .AND. IGL0C.EQ.2) 

IN-LINE GAP HCF/JSC (tO.25) 



.AND. 

IHOF.EQ.l ) 


ICOR 


3 


60 TO 

50 





.AND. 

IHOH.EQ.2) 


ICOR 


4 


GO TO 

50 






.AND. 

IHOH.EQ.I . AND.GAPFL.LE. 

15.24) 

ICOR 


5 


GO TO 

50 






.AND. 

IHOH.EQ.I .ANO.GAPFL.GT. 

15.24) 

ICOR 

s 

6 


GO TO 

50 






.AND. 

IHOH.EQ.I) 


ICOR 

r 

7 


GO TO 

50 





.AND. 

IHOH.EQ.2 .AND. GAP'LL. LE. 

15.24) 

ICOR 

- 

8 


GO TO 

50 





.AND. 

IHOF.EQ.2 .ANO.GAPFL.GT. 

15.24) 

ICOR 

z 

9 


GO TO 

50 






.AND. 

IHGH.EQ.2 ) 


ICOR 

= 

10 


ICOR = 11 


IFdBL.EQ.l .AND. IGLOC.EQ.3) ICOR = 13 

IFdCOR. NE. ID ) GO TO 50 

TRANSVERSE GAP, DOWN STREAM WALL, TURBULENT (EQ. 18) 

IF(IBL.EQ.3 .AND. IGL0C.E0.2) ICOR = 14 

IFdCOR. NE. ID ) GO TO 50 


EFFECT OF GAMMA ON TRANSVERSE GAP, TURBULENT (EQ. 12A) 
IF(IBL.EQ.3 .AND. IGLOC.EC.2 . AND . GAMMA. NE. 0 . 0 ) ICOR 

IFdCOR. NE. 10 ) GO TO 50 
IF (IGL0C.NE.6) GO TO 14 
TFdCOR.NE.ID ) GO TO 50 


LONG IN-LINE GAP CORRELATION (EOc 

IF(GAHMA.EQ. 0.0 .AND. IHOH.EQ.I) 
IFdCOR. NE. ID ) GO TO 50 
IF(GANMA.EQ, 0. Q .AND. IHOH.EQ.2) 
IFdCOR. NE. ID ) GO TO 50 


16-1 TO 


IFdGLOC.EQ. 1 
IFdCOR. NE. ID ) 
IFdGLOC.EQ. 2 
IFdGLOC.EQ. 1 
IFdGLOC.EQ. 2 
IFdCOR. NE. ID ) 


AND. IHOH.EQ.I) 
GO TO 50 

.AND. IHOH.EQ.I) 
.AND. IHOH.EQ.2) 
.AND. IHOH.EQ.2) 
GO TO 50 


EQ 16-6) 
ICOR 

ICOR 

ICOR 

ICOR 

ICOR 


= 15 


16 

17 

18 

19 

20 
21 


STEPS 

14 IF (STEP.EQ.C .0) 
IFdBL.EQ.l) 
IF(I8L.EQ.2> 
IF(IBL.FQ.3) 


GO TO 50 


ICOR 

ICOR 

ICOR 


22 

23 

24 
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c 

c 


50 


c 

c 

c 


c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 


55 

70 

80 

90 

94 


1010 

1011 

1012 

1020 


1022 


C 

c 

c 

c 

c 


c 

c 

c 


c 

c 


c 


1030 

104C 

105G 


lUATION SELECTED 
^ONTINUE 

Z=ZZ(H 
GO TO 94 

IF( S.LT.SMIN) GO TO 70 
IF( S.GT. SHAX) GO TO 80 
GO TO 94 

HEATING RATIO EVALUATED 

SECOND ORDER CURVE FIT BETWEEN SHIN AND S=0 
XLNQR = (BHIN CHIN*S>*S 

GO TO 90 

LINEAR EXTRAPOLATION INTO GAP BELOW INSTRUMENTATION 
XLNQR = FMAX DFOSB* ( S -SMAXl 

Q = EXP<XLNQR) 

GO TO 20 00 

CONSTANTS COMPUTED DURING ITRIG=1 PASS 

BRANCHING TO EQUATION SUBROUTINES 


GO TC 


{ 1 010 1 1020* 10 30. 1040,10 50 «1060« 10 70. 1080, 1090«1 


100 , 

200 . 


table LOOK-UP 
IF(ITRIG.EQ.2) GO TO 1012 
WRITE(6,901> 

DO 1011 K=1,ITQ 

WRITE 16,950) K.ZTAB(K) ,QTA8(K) 

CONTINUE 
ITRIG = 2 
IMEM = 2 
GO TO 55 

CALL TBLUdOUM, IMEM , Z, Z T AB, Q. QT AB) 

GO TO 2000 

SCOTT AVERAGE OF TRANSVESE GAP LP AND DOWN HALLS 

rF(ITRIG,EQ.2 ) GO TC 1022 
HRITE(6.902) 

IR = 1 

IFfE.GE* 0.3048 ) IR = 2 
GO TO 55 

CALL CSCOTT(IR.S.Q) 

GO TO 2000 


TRANSVERSE GAP, DOWN IEQ.4-17), LAMINAR, (H/HFP) 

IF (ITRIG, EQ. 1) MRITE(6,903) 

CALL EQ41716(ITRI6,S,Q, ICOR) 

GO TO 2000 

transverse gap, down (EQ. 4-16) , laminar, (H/HE) 

IFdTRIG.EQ.l) HRITE(6.904) 

CALL EQ4l716(TRIG,S,Q,lCOR) 

GO TO 2000 

IN-LINE GAP UP-STREAM OF JUNCTION, LAMINAR (H/HFP ), EQ .4-9 
IFdTRIG.EQ.l) MRIT£(6,905) 

CALL EQ409(ITRIG,S,0,ICCR) 

GO TO 2000 

C5 
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IN-LINE GAP OMN-STREAH. LAMINAR CH/HFP), £0.4-10 
1060 IF«ITRIG.E0.1» HRITEC6,§06I 
CALL EQ4010aTRIG«S,Q,ICORI 
GO TO 2000 


REPORT MDC E1248 
JSC 09651 


1070 


C 

C 

C 

c 


2000 

TRANSVERSE GAP, UPSTREAM, LINEAR 
IF CITRIG.EQ.l) HRITEf6, 907l 
CALL EQ4D11CTRIG,S,0) 

GO TO 2000 


JH/HFP), (EQ.4-11) 


iOdO 


IN-LINE GAP IH/hE>, UP, !EQ. 

TFaTRIG.EQ.lI MRITE^6,'.J8> 
CALL EQ4013(ITR1G,S,Q) 

GO TO 2000 


4-131 


n 


1090 


C 

C 

c 


iiOU 


C 

C 

C 


C 

C 

C 


IN-LINE GAP (H/HE), OOMN, ( EQ.l— 141 

IFdTRIG.EQ.ir MRITE(6,909< 

CALL EQ4014CITRIG,S,Q) 

GO TO 2000 

TRANSVERSE GAP(H/HE», UP CEQ.4-1SI 

IF (ITRiG.EGill NRITF(6,910I 
CALL EQ4015(ITRlG,5,a; 

GO TO 2000 

TRANSVERSE GAP, OOMN, H/HFP, LAMINAR, 3TEST PROGRAMS, i£G 15-9) 

IF (ITRIG.EQ.l) WRITE(6,911» 

CALL EQ1505(ITRI3,S,Q) 

GO TO 2000 

TRANSVERSE GAP, DOWN, H/HE, HCF/JSE(£Q. 9) 

1120 IFIITRIG.EQ.I) MRITE(6,912) 

CALL EQ9(1TRIG,S,Q) 

GO TO 2000 


C 

C 

c 


1110 


IN-LINE GAP (H/HE) 


1130 


,HCF/JSC (EQ.25) 
(6,913) 


C 

c 


1140 


C 

C 

C 


IF(ITRIG.EQ.l) MRITE 
CALL EQ25 (ITRIG,S,Q) 

GO TO 2000 

TRANSVERSE GAP TURBULENT (CFAT,Ms8), 
IF(ITRIG.EQ.l) MRITE (6,914) 

CALL FQ18(ITRIG,S,Q) 

GO TO 2000 


(EQ. 18) 


1150 


C 

C 

C 

C 

c 


1150 


c 

c 

c 


1170 


INFLUENCE OF GAMMA, TRANSVERSE GAP (EO. 12A) 

IF (ITRIG.EQ.l) MRITE(6,S15) 

CALL EQ12A ( ITRIG , Z , Y, Q ) 

GO TO 2C 00 


H/HFP , LONG IN-LINE GAP , LAMINAR, GAMMA=0 (EQ.16-1) 

IF (ITRIG. EQ.x) MRITE(6,£16» 

CALL EQ16D1(ITRIG,Z,Q) 

GO TO 2000 

H/HE , LONG IN-LINE GAP, LAMIFAR, GAMMA=0 (EQ.16-2) 

IF (ITRIG.EQ.l) MRITE(6,917) 

CALL EQ1602(ITRIG,Z,0) 

GO TO 2000 
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H/HFP, LONG IN-LINE GAP, LAHINAR, GAMMA UP TO 15, SIDES A/B EQ16-3 

1180 IFCITRIG.EO.l) MRITEC6,918I 
OALL EQ16034(ITRIG,Z,Q) 

GO TO 20 00 

H/HFP, LONG IN-LINE GAP, LAMINAR, GAMMA UP TO 15, SIDE A (EQ.16-4) 

1190 IF(ITRIG.EQ.l) WRITE (6,9191 
CALL EQ16034(1TRIG,Z,QI 

GO TO 200C 

H/HE, LONG IN-LINE GAP, LAMINAR, GAMMA UP TO 15, SIDES A/B EQ.16-5 

1200 IFdTRIG.EQ.ll WRITE (6,920) 

CALL EQ16056(ITRIG,Z,Q) 

GO TO 2C00 

H/HE, LONG IN-LINE GAP, LAMINAR, GAMMA UP TO 15, SIDE A (EQ.16-6) 

121C IFdTRIG.EQ.ll WRITE (6,9211 

CALL EQ16D56( ITRI :,Z,QI 
GO TO 2000 

1220 IFdTRIG.EQ.ll WRITE(6,S22I 
CALL EQ17123(ITRIG,Z,Y, Cl 
GO TO 2000 

123u if ilTRIG.EG.l) WRITE(6,923I 
CALL EQ17123(ITRIG,Z,Y,(]| 

GO TO 2000 

1240 IFdTRIG.EQ.ll WRITE(6,924I 
CALL EQ17123(ITRIG,Z,Y,QI 
GO TO 2000 
1250 CONTINUE 
1260 CONTINUE 
1270 CONTINUE 


203C CONTINUE 

IF (ITRIG.EQ.2) GO TO 2010 


ITRIG = 2 

CMIN = -(FMIN -SMIN^OFDSA I/{SMIN*SMIN) 
0MIN = (FMiN-CMIN*SMIN*SMlNI/SMIN 

GO TO 55 
2010 CONTINUE 
C 
C 

IF (Q.LT. HOHLOW) Q =0.0 
IF(Q.GT. HOHHI) Q=HOHHI 
C 

QGd) = Q 

IF( I.EQ. NTAB) GO TO 3000 
I = !♦! 

Z - ZZ(I) 

GO TO 55 

3000 MRITE(6,952) GAPW 
OO 3010 1=1, NTAB 
QR = 00 ( I) 

QQdl = QO(II*HFP 
WRITE(6,E54) I , ZZ ( 1 1 «QQ ( 1 1 , QR 
3010 CONTINUE 
RETURN 
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FORHATC*! USER SUPPLIED PARAMETERS TO GAP HEATING DISTRIBUTION*) 
FCRMATC* HEATING RATIO INPUT TABLE*/* I Z(CM) QRATIO*/) 
FORMAT (* AVERAGED GAP HEATING IN TRANSVERSE GAP CC. SCOTT DATA)*) 
FORMAT (• CORRELATION OF TRANSVERSE GAP , E-EFFECT <EC. A-17)*) 
FORMAT!* CORRELATION OF TRANSVERSE GAP , E-EFFECT (EQ.4-1G H/HE)*) 
FORMAT!* IN-LINE GAP»H/HFP, UP-JTREAM OF JUNCTION !EQ. *♦-9)*) 
FORMAT!* IN-LINE GAP.H/HFP, OH-STKEAM OF JUNCTION ! EQ . 4-10)*) 
FCRMAT!* TRANSVERSE CAP, UP. H/HFP !EQ. 4-11)*) 

FORMAT!* IN-LINE GAP, UP, H/HE IEQ.4-13)*) 

FORMAT!* IN-LINE GAP, CN H/HE !EQ.4-14)*) 

FORMAT!* TRANSVERSE GAP, UP, H/HE !EQ. 4-15)*) 

FORMAT!* transverse GAP, OWN, H/HFP, LAMINAR !EQ,lS-5)*) 

format;* transverse gap, down, H/HE , HCF/JSC (EQ.9)*) 

FORMAT!* IN-LINE GAP (H/HE), HCF/JSC CEO. 25)*) 

FORMAT!* TRANSVERSE GAP TURBULENT CFHT,M=8 (EQ.18) *) 

FORMAT!* TRANSVERSE GAP TURBULENT CFHT,M=8 (EQ.18) *) 

H/HFP, LONG IN-LINE GAP, L AMIN AR.GAMMA = 0 (EQ.16-D*) 
FORMAT (• H/HE, LONG IN-LlNE GAP, LAMINAR, GAMMA=0 (£0.16-2)*) 
FORMAT(* H/HFP, long IN-LINE GAP, LAMINAR, GAMMA=UP TO 15, SIDES A/B 
( EQ »1& — 3 )*) 

FORMAT (* H/HFP, LONG IN-LINE GAP, LAMINAR, GAMMA UP TO 15, SIDE A 
. (EQ.lo-4)*) 

FORMAT!* H/HE, LONG IN-LINE GAP, LAMINAR, GAMMA UP TO 15, SIDES A/B 
( EQ . 6-5 )* ) 

F0RMAT^*^H/HE, long IN-LINE GAP, LAMINA fi, GAMMA UP TO 15, SIDE A 

FORMAT (• TRANSVERSE STEP, LAMINAR FLOW*) 

FORMAT (♦ TRANSVERSE STEP, TRANSITION *) 

FORMAT! *TRANSVERSE STEP, LTURBULENT *) 

FCRMAT ( I4,2F10. 5) 

FORMAT!* COMPUTED HEATING DISTRIBUTION GAP H ID TH= *F8 . 3 , • CM*/ 

♦ I Z(CK'/ Q QRATIO*) 

FORMAT(I4,3F10.5> 
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SUBROUTINE CSCOTT I IR.SCH, QR> 

GAP HEATING OISTRI0UTIOK DATA PROM CARL SCOTT 
FOR GAP WIDTH = 0.10 INCH 

Q/QREF = FIS, EDGE RAO» (DIMENSIONS ARE IN INCHES) 

DIMENSION SCOTTSI18.2) , SCOTTQ (1 8.2 ) 
giHeNSiON^CSO(2).CSi(2) ,CS2(2),CS3I2) 

DATA ISCOTT / 18 / 

DATA SCOTTS / -100. 0, -0.26, -0.14. -J.IQ, -0.08, -0.04, -0.02, 

1 —0.01, 0.0* 0.06, O.lO, 0.16, 6.20, 0.28, 1.0, 2.0, 

2 10,0, lOfl.u, 

3 -lOC.O. -0.30, -0.2C, -0.15, -0.10, -0.06, -0.02, 0.0, 0.02, 

4 0.06, 0.10. 0.15. 0.20, 0.22. 0.30. 0.36, 0.40, 100.0 / 

DATA SCOTTQ / 0.93, 0.93. 0.94. 0.96, 0.99, 1.02, 1.03, 1.00, 

1 3.96, 0.80, 0.48, 0.36, 0.24, 0.0, 4*0.0, 

2 0.95, 0.95, 0.96, 0.97, 0.98, l.Ou, 1.02, 1.01, 0.98, 0.95, 

3 0.90, 0.80, 0.55, 0.45, 0.27, 0.08, 0.0, 0.0 / 

DATA CSC / 0.965, 0.975 / 

DATA CSl / 0.035, 0.0 25 / 

DATA CS2 / 0.06, 0.30 / 

DATA CS3 / 2.88, 2.575 / 

S s SCM/2,54 

IF(S.LT.SCOT1S.2,IRI) GO TO 156 

call TBLUIIOUM.IMEM, ISCOTT, S, SCOTTS 1 1, 1 RI ,QR, SCCTTQI I, IR) ) 

GO TO 157 

156 QR = CSO(IR)-CS1IIR)*COSIPI*C-S-CS2(1R) )-’CS3CIR)) 

157 CONTINUE 
RETURN 
END 
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SUBROUTINE EQ41 71 6( IT RI G , S ♦QR ♦ ICORR) 

GAP HEATING DISTRIBUTION USING EQUATIONS 4-17 AND 4-16 

ICORR = 3, EQUATION 4-17 
ICORR = 4, EQUATION 4-16 

COMMON / COMGAP / ARGL ( 25 ) .SM IN •FMIN , OF CS A, SM AX ,F MAX t DFOSB 
DIMENSION CI5» 

EQUIVALENCE (AR6L(5l*EI , ( AR6L C7 I • M I * (ARGL(21>tT) 

DIMENSION XEC4) ,XSMAX(4 > »XSMINC4) 

DATA XE / 0.15748* 0.318. 0.635. 1.270 / 

DATA XSMIN / 0.400. 0.540. 0.810. 1.350 / 

DATA XSMAX / 0.90 0. 1.050. 1.850. 3.000 / 

DATA IS / 4 / 

11 FORMAT (6X.12HC0EFFICIENTS.2X. 5F10.5/6X. lOHSMI N, QRMIN . 4X . 2F10 . 5/ 

1 6X,10HSMAX.ORMAX.4X.2F10.5» 

12 F0RMAT(6X.9HCMIN,BMIN.5X.2P10.5» 

F = C1*C2*S*C3*S** 3+C4/S**2*C5*LN<S» = LN(Q/QREF) 


INITIALIZE 

IF (ITRIG.GE.2> GO TO 100 

IF(IC0RR.£Q.4> GO TO 20 


•CONSTANTS FOR EQUATION 4-17 


0(1» = C.59775^3.05060*E 
C(2I = -3.67827-0.09307/M 
C(3> = 0.14962 
C(4) = -0.27898*Z*e 
C(5) = 0.0 
GO TO 40 


CONSTANTS FDR EQUATION 4-16 


20 CONTINUE 

cm = 0.42466«-3. 84088* (M^ei-2.18027*M*E*E 

C(2) = -5.14719 

C(3» = 0 .19450 

C(4) = -0.3C623*E*E 

C(5) = 2.27613*^ 

DETERMINE VALUES AT S EXTREMES OF CORRELATION 
INTERPOLATE TO OBTAIN SMAX AND SMIN = F'.E) 


40 CONTINUE 
IMEM = 2 

CALL TBLUdCUM.IMEM.IS. E.XE. SMAX. XSMAX) 

Call TQLU( IOUH, IMEM, is. E.XE. SMIN. XSMIN) 

QRMAX = QRFIC.SMAX) 

ALNQRMX = ALOG(QRMAX) 

DALNQS = 2.0*(T-SMA'')*ALOGf0.1D) 

OFOSB = DALNQS 

COEFFICIENTS FOR CORRELATION 

AT LOW ENC OF CURVES 
MATCHING OFOS AT SMIN 

OFOSA = C(2)*3.0*C(3)*SMIN*SMIN-2.0*C(4)/ (SMIN**3)*C(5)/SMIN 
FMIN = ALOG(QRMIN) 

GO TO 160 

DETERMINE GAP HEATING RATIO 

100 CONTINUE 

QR = QRF(C.S) 


DALNQS = SLOPE OF 1 CYCLE PER 0.5 CM FROM SMAX TO T 


160 RETURN 
END 
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FUNCTION QRFCC.S) 

CORRELATION EXPRESSION TO EQUATIONS 4-17 AND 4-16 
DIMENSION CC^ ) 

XLNQR = C(l)i-C(2»*S*C(3)*S*S*S*CI4|/<S*S»>C(5)«ALOG(S) 
QRF = EXP(XLNQR) 

RETURN 

ENO 


SUBROUTINE TBLU (lOU M, I . IT , X, XT AB , Y , YTAB ) 

ONE DIMENSIONAL TABLE LOOKUP WITH A MEMORY 

DIMENSION XTAB(l) .YTABCl) 

C FIND BRACKET 

IF(X.LT.XTAB(IT-in GO TO 8 
I=IT 

GO TO 33 

8 IFCXTABm.GT.XI GO TO 20 
10 I=I^1 

IFd.EQ.ITl GO TO 30 
IF(XTA3(I».GT,XI GO TO 30 
GO TO 10 
C 

20 IF (XTAB (I-l) .LE.X » GO TO 30 
1 = 1-1 

IF(I.EQ.2l GO TO 30 
GO TO 20 
30 CONTINUE 

C I DESIGNATES THE XTAB POINT TC THE RIGHT 

AB*<X-XTAB(I) l/(XTABCI-ll-X'; ABl in 
Y=AB*YTAB(I-D* (1.0-AB) »YTAB(li 
RETURN 
END 


cn 
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SUBROUTINE EQ4091 ITRIG, S, Q) 

IN-LINE GAP. UPSTREAM OF INTERSECTION. >- AMINAR^^^^g/gFP)^^^ 

£ = EDGE RADIUS <CM) 

M = GAP UIOTH (CM) 

S = SURFACE DISTANCE INTO GAP <CH» 


COMMON /COMGAP/ ARGUZ 5 ) , SHIN ,F HI N, OFOSA . SMAX ,FMAX. OFOS0 
EQUIVALENCE (ARGLtSl.EI. (ARGLf7)fM) 

OIMENSICN C(5) 

IF (ITRIG.GE.2 ) GO TO 100 

CONSTANTS FOR EQUATION 4-9 

C(l» = 0.39439 ♦0.76083*E/M - 0 . 02602*E*E/ ( M*M ) 

C(2» = -0.71672/W-0. 51625 
C(3) = 0.0132/(M*M» , . ^ 

C(4) = 2.53332*E*E -3.62726*M*E 
C(5» = -1.33050*E ♦ 2.53332*N*E*E 


INTERPOLATE TO OBTAIN SMIN ANO SMAX 


SMIN = SMINFIEI 
SMAX = SMAXF(E) 

QRMIN = QRF409(C,SMIN» 

FHIN = ALOG(QRHIN) 

QRMAX = QRF409(C.SMAX) 

OFDSA = C*??) ‘!’^2.0*C(3) »SMIN ♦ CC4)/SMIN ♦ 2 , 0*C <5 ) * ( ALOG (SHINU / 
1 ^HTN 

OFDSB = C(2) ♦ 2.0*C!3I *SMAX ♦CUI/SMAX ♦ 2 . 0 *C (5 ) * ( ALOG ( SM AX )» / 

1 SMAX 

RETURN 

100 Q = QRF409(C.S) 

RETURN 

END 


FUNCTION SMAXF(E) 

MAX. SURFACE DISTANCE FOR INSTRUMENTATION IN C. SCOTT TESTS 
DIMENSION X£ (4) ,XSMAX(4 ) 

DATA XE / 0,00254, 0.318, 0.E35, 1.270/ 

DATA XSMAX / 0,900 . 1,050, 1.850, 3.000 / 

DATA IS.IMEM /4,2/ 

CALL TBLUdOUM, IME M,I S , E , XE, AB , XSMAX) 

SMAX = AB 

return 

END 


FUNCTION SMiNFfE) 

MIN. SURFACE DISTANCE FOR INSTRUMENTATION IN C. SCOTT TESTS 


DIMENSION XE(4). XSHIN(4) 
dUa XE / 0.00254, 0.318. 0.€3«, 

DATA XSMIN / 0.400 , 0.54C, 0.810, 
DATA IS.IMEM /4,2/ 


1.270/ 

1.350/ 


CALL T9LL(I0UM,IMEM,IS, E.XE.AB, 

SMIN = AB 

RETURN 

END ^^2 


XSMIN) 


' * %*'**'i-*' • r AST 
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FUNCTION QRF409CC,S» 

C CORRELATION FUNCTION 4-9 

DINENSION C(l) 

ALNS = ALOG(S) 

XLNQR = C(ll ♦ (C(2» ♦ C(3I*S»*S ♦<C(4lf C (5» *ALNS) »ALNS 
QRF409 = EXP(XLNQR) 

RETURN 

END 
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SUBROUTINE EQ401(i ( ITRIG ,S,QJ 

IN-LINE GAP, OOMNSTREAM OF INTERSECT ION , L AMIN AR , (Q/QFP) 

H.E. CHRISTENSEN 75 

E=E0GE RADIUS (CM) 

W=GAP WIDTH (CM) 

S= SURFACE OISTA ICE INTO GAP(''M) 

COMMON/COMGAP/ ARGL(25) , SHIN, FMIN ,0F OSA ,SM AX, FM AX , DF0S8 
EQUIVALENCE (ARGL (5) ,E) , (ARGL (7I,W) 

OINENSICN C(5) 

IF(ITRIG,GE.2 ) GO TO 100 


CONSTANTS FOR EQUATION 4-10 


C(l) =-1.21264 ♦(7,82944 - 2.30047 ♦£)•£ 

C(2) =-0.060260/W -4.27622 

C(3) = 0.125Q1 

C(4) =-a.3939Q*E*E 

C(5) - 1,0 9199“M 


IN7ERPLETE TO OBTAIN SHIN AND SMAX 


C 

C 


SMIN = SMINF(E) 

SMAX = SMAXF(E) 

ORMIN = QRF4O10 (C.SPIN) 

FMIN = ALOG(QRMIN) 

QRMAX = QRF401C (C.SMAX) 

FMAX = ALOG(QRMAX) 

OFO^A = U(2I ♦ 3,0*C(3I*SMIN*S 'N -2 . 0 • C ( 4) / ( SMIN** 3 ) ♦ 2.0*C(5)* 
1 ALOG(SMIN) /SMIN 

DF0S8 = C(2)^ 3,B''C(i) *SMAX*SMAX -2 . 0*C (•♦) / (S MAX**3 ) ♦ 2.0*C(5l* 

1 ALOG(SMAX I /SMAX 

RETURN 

100 Q = QRF4O10(C,S) 

RETURN 

END 


FUNCTION QRF401C(C,S) 

XLNQR^i*'nii ♦ic(2) ♦C(3)*S*S.*S ♦ C(4)/(S*S) ♦ C (5 ) * ( ( ALOG ( S ) ) * 2 ) 
QRF4010 = EXP(XLNQR) 

RETURN 

END 
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VOLUME I 


REPORT MDC EI24H 
JSC 096SI 


SUBROUTINE EQ40 IH ITR IG , S » Q» 

TRANSVERSE GAP UP-STREAM HALL, LAMINAR, lO/QFP) 

H.E. CHRISTENSEN 

E = EDGE RADIUS IQM) 

W = GAP WIDTH fCM) 

S = SURFACE DISTANCE INTO GAP (C^i 

COMMON / COMGAP/ ARGL (2 5) , SMT N, FMIN, OFDSA ,SM A X , FM AX , DFDSB 
EQUIVALENCE (ARGLI5),E)« (ARC (71, M) 
iTIMENSICN C(3) 

IF (ITRIG.GE.2) GO TO 100 

CONSTANTS FOR EQUATION 4-11 

C(l) =-0.31375f> * 2.27958*E ♦ 0 , 33317* ALOG ( W) 

C«2» =-4.1428 ♦ 1.3586l*e 

C(3» =-1.42212 

3MAX = SMAXF(El 

QRMAX = QRF4011(C,SMAX) 

FMAX = ALOG(QRMAX) 

SMIN = SMINF(E) 

QRMIN = QRF4011 CC.SMINl 
FMIN = AL0G(QRMIN» 

OFDSA = C(2)/SHIN ♦ 2. C*C (3) *ALOG (SMIN) /SMIN 
OFOSB = C(2)/SMAX ♦ 2 . C*C ( 3) * ALOG (SM AX ) /SM AX 

100 Q = QRF4011(C,S) 

RETURN 

END 


FUNCTION QRF4D11(C, S) 

OIHENSICN C(l) 

AQ = ALOG(S) 

XLNQR = C(1)*(C(2)* C(3!*A0I*AB 
QRF4011 = EXP(XLNQR) 

RETURN 

END 
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SUBROUTINE EQ4013 (ITRIG ,S,QI 

IK-LINE GAP, UPSTREAN OF INTERSECTION, LAMINAR, (Q/OEJ 


COMMON /COMGAP/ A RGL C2 5 I , SHIN ,FMIN, OFOSA, SM AX ,FM A X, OFOSE 
EQUIVALENCE (ARGLI5),E), (ARuL(7t,M) 
dimension CC6) 

IFCITRIG.GE.Z) GO TO 100 


-CONSTANTS FOR EQUATION 13 


C(l) 

C(2» 

C(3l 

C(4) 

C(5J 

C(6) 


=-0,01221 ♦0.27744*E/N 
=-1,37104 -0,4iq40/H 
= 0, 00725/<M*M» 

= C, 04555 
= 0,80703*e 
=-0,47496*E*E 


♦ 0,70519*( <E»MI**2> 


SHIN = SrtlNFJEl 
SMAX = SHAXF(E) 

QRMIN= QRF4013IC,SMIN) 

FMIN = ALOG(QRMIN) 

QRMAX = QRF4013(C,SMAXI 
FHAX = ALOGIQRHAXl 

OFOSA= C(2»^ 2,0»C ( 3MSMIN ♦ 3. 0*C '4 »*SMIN*SMIN -C (5 » / C SMIN*SMI K) 
1 -2.0*C(6>/(SMIN**3) 

OFDS0= C(2) v2,0*Ct3)*SHAX ♦ 3,0*CI4) ♦SMAX*SNAX - C < 5 1 / ( SMAX* SMA Xf 
1 -2, C*CI6l/(SMAX**3l 

C 

RETURN 

C 

100 Q = QRF40l3fC,SI 
RETURN 
END 


C 


FUNCTION QRF4013(C 


,SI 

CORRELATION 


FUNCTION 


4-13 


OIMENSICN C«1J 

XLNQR= Cm*lC(2»+ {C<3 » ♦Cf4»*SI*S! *S 
QRF4D13 = EXP(XLNQR) 

RETURN 

END 


♦CC<5»* C<6?/SI/S 
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GAP HEATING ANALYSIS - 

VOLUME I 


I 


SUBROUTINE EQ4014I ITRIG ,S,0> 


IN-LINE GAP, OOHNSTREAM OF INTERSECTION, LAMT ;jAR, (Q/QE > 

COHMON/COHGAP/ ARGL (25 ) , SHIN, FHIN, OFDSA ,SHAX, FMAX ,OFOSB 
EQUIVALENCE (ARGL(5),E), (ARGL(7),Vn 
01 PENSION C(5) 

IF(ITRIG.GE.2 ) GO TO 100 

CONSTANTS 

C(l)=-0. 16717 ♦ 1.45275*( (M+E»**2» ♦ 0. 42637*ALCG(EI 
C(2)=-2,6609 
C(3)= 0.06483 
C(4»- 0 ,33840^E*E 
C(5)=-1.3285 
C 

SMIN = SMINFCEI 
SMAX = SMAXF(EI 
C 

QRMIN= QRF4014(C,SMIN) 

FMIN = ALOG(QRMIN) 

QRMAX = QRF4D14(C,SMAX) 

FHAX = ALOG (QRMAX) 

C 

DFOSA = C(2) ♦ 3,0*C(3)*SMIN*SMIN -2 ,0 *C (4) / ( SMIN**3 ) ♦ C(5)/SMIN 
OFDSa = C(2) f 3. 0*C (3) •SMAX^'SMAX - 2 . 0 *C (4 ) / (SMAX** 3 ) ♦ C(5)/SMAX 
C 

RETURN 

C 

130 Q = 0RF4D14(C,S) 

RETURN 

ENO 


FUNCTION QRF4D14(C,3) 

C CORRELATION FUNCTION 4-14 

DIMENSION C(l) 

XLNQR = C(l) ♦ (C(2) ♦ L(3)*S*S)*S ♦ C(4)/S ♦ C(5)*ALOG(S) 
QRF4014 = EXP(XLNQR) 

RETURN 

END 
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VOLUME 


REPORT MDC E1248 
JSC 09651 


subroutine: EQ<*015 (ITRIG«S,Q) 

TRANSVERSE GAP <UP-STREAH SIDE) t LAMINAR, CQ/UE) 

COMMON /COMGAP/ ARGL (25 ) , SH IN, F MIN, OFOS A, SMAX ,FMAX, 0F0S8 
EQUIVALENCE (AR6L(5),EI ,(AR6L(7),M) 

DIMENSION C(3) 

IF (ITRIG.GE.2* GO TO 100 

CONSTANTS 

C(l)= 0.35564 ♦ 1,18292*E*E 
C(2l=-2.7736 
C(3I= 1.83519»H*E*E 

C 

SMIN = SMINF(E» 

SMAX = SHAXF(EI 

C 

QRMIN= QRF4015(C,SHINI 
FMIN = ALOG(QRMIN» 

QRMAX= QRF4015(C,SHAX» 

FMAX = ALOG(QRMAX) 

OFDSA = C(2I ♦ C(3I/SMIK 
DFOSB = C(2» C(3»/SMAX 
C 

RETURN 

C 

IGO Q = QRF4D15(C,S» 

RETURN 

END 


FUNCTION QRF4015(C,S) 

dihensicn cm correlation function r-ie 

SiifSoM * C.3...LCG.S. 

RETURN 

END 
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REPORT MDC E1248 
JSC 096S1 


RSI GAP HEATING ANALYSIS - II 

\ VOLUME 


SUBROUTINE EQ15C5 ( ITRIG » Z , 0) 

TRANSVERSE GAP, 30UNSTR fAP HALL (LAMINAR ) (H/HFP) (EQ 5-15» 

COMMON /COMGAP/ ARGL (2 5 ) , ZMIN,FMI N. DFOZ A, ZMAX, FHAX, DFDZ3 
EQUIVALENCE (ARGL(3),E) , (ARGL (7 ) , Ml , (ARGL (6 ) , T) , ( ARGL (14) ,REPH) 
OIHENSICN C(4) 

RE = REFM/(1.0£+36) 

IFdTRIG.GE.Z) GO TO 100 

CONSTANTS FOR EQ.15-5 

C(l)=-D, 40577 ♦ 2.gi733*M -3,03105*M*H 
C(2)=-6, 80369 - 0.36933*RE 
C(3)= 1.5057 ♦ 0.38980»RE ♦ 1.79532*M*H 
C(4)=-0 .11073*H*RE*T 

ZMIN = -1,E*3C 
ZMAX = *.,0 
HRITE(6,990I ZMAX 
QRMAX = QRF1505 (C.ZMAX) 

WRITE (6, 9001 ZMAX,C,QRMA '< 

FMAX = ALOG(QRMAXl 


DFOZB = C(2) ♦ (2.0*C(3) *■ 3. 0«C (4) *ZMAX) ♦ZMAX 

WRITE(6,900) DFDZ9 

RETURN 

100 Q = QRF15D5 (C.Z) 

RETURN 

900 FORHAT(* ZMAX, C, QRMA X*/8E13. 3 ) 

END 


FUNCTION QRF1505(C 


'I 


) 

ORRELATION FUNCTION 15 


5 


OIMENSICN C(l) 

XLNQR = C(l) ♦ (C(2)»(C (3)>C(4) *Z)*Z)*Z 
QRF1505 = EXP(ZLNQR) 

return 

END 


C 

C 


SUBROUTINE £Q9( IT RIG,Z , 0) 

TRANSVERSE GAP, LAMINAR (E09) HCF TESTS AT JSC 

COMHON/COMGAP/ ARGL(25), ZMIN ,FMIN,uIFOZA , ZMAX, FMAX , DFOZB 
EQUIVALENCE (ARGL(7),W), (ARGL(8),T) 

DIMENSION C(3) 

IFdTRIG.GE.Z) GO TO 100 

CONSTANTS 

C(H =-C.3716 «-1.0733*H 

C(2) =-5.0249 -0.3797*T*T -0.03654/H 

C(3) =2.5604 -6.-7l9/(T^T) 

RETURN 

108 Q = QRF9(C,Z) 

RETURN 

END 


CIS 
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RSI GAP HEATING ANALYSIS - 

1^1 \ VOLUME 


FUNCTION QRF9(C,ZI 

CORRELATION EQUATION 9 

01 KENS I CN Cfl) 

XLNQR = C(l) ♦ «C(2) + C13»*Z)*Z 

QRF9 = EXPtXLNQRI 

RETURN 

END 


REPORT MDC E1248 
JSC 09651 


c 

c 

c 


c 

c 

c 


c 

c 


100 


SUBROUTINE EQ25 (I TRIG, Z »Q) 
IN-LINE GAP, HCF TESTED AT JSC 


(H/HFP), LAMINAR 


COMMON/ COMGAP/ A RGL (2 5 ) , ZMIN ,F Ml N , OFD Z A , ZM AX , FM AX ,DFD ZB 
EQUIVALlNCE (ARGL(7),Mi 
DIMENSION C(3) 


CONSTANTS 


Ctn =-C.3319 ♦ 
CC2» = -4.3979 
C(3) = 1.563 

RETURN 

Q = EXP( C(l) ♦ 

RETURN 

END 


1. J148*H 


2295/H 


(C<2» ♦ C(3)*Z)*Z) 


C 

C 

c 


c 

c 

c 

c 

c 


c 

c 


SUBROUTINE E3 1 8 ( ITR TG, Z ,0 ) 
TRANSVERSE GAP (CFHT M=1Q, VDT M=8) 


TURBULENT FLOW (EQ.18) 


100 


COMMON/COMGAP/ AR GL (25 I , ZMIN, FMIN , OFDZA , ZMAX, FM AX ,OFOZB^^^ ^ 
EQUIVALENCE IAR6L ( 7»,W», (ARGL(19), TWOTEJ, IARGL(14) ,REPM> 
DthenSION C(3> 

IF (ITRIG. GE.2 ) GO TO 100 
CONSTANTS FOR EQUATION 18 

C(l) = 4.28021 -0.53620*(ALOG(R£PM/1.0E+06M *0 .7 806* ALOG 5 M) 
C(2» = 1.46402*(W**(-2.0/7.0n 
C(3) =-2.549 

RETURN 

Q = EXP (C(l) fCIZ) *ALOG( Z> ♦ C ( 3) *ALOG (THOTE )) 

RETURN 

END 


C 

C 

C 


C 

c 

c 

c 

c 


IQO 


SUBROUTINE EQ12 A ( TRIG, Z , Y , Q ) 

INFLUENCE OF FLOW ANGLE ON TRANSVERSE GAP HEAT ( TURB ),( H/HFP ) 

COMMON /COMGAP/ ARGL ( 2 5) , ZMI N, FMI N.OFDZA , ZMA X, FM AX , OFD ZB 
EQUIVALENCE (AR6L(7),N), (ARGL(e>,Cl 

IF (ITRIG, GE.2 ) GO TO 100 

CONSTANTS 

C = -0.15C4 ♦ 1.92179*M 

RETURN 

Q = £XP(C -3. 714*7 ♦0.60432*Z*Z -0.00692*Y*V ♦0.04316*Y*M 
1 ♦ 0.iH 8*Y*6 ♦ 0,00235*V*y*Z» 

RETURN 

END C19 
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c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


SUBROUTINE EQ1601 ( ITRIG , Z t Q> 

LONG IN-LINE GAP, LAMINAR, GAMMA-0 (EQ.16-1) 

CCMMON /COMGAP/ ARGL C2 5 ) , ZM IN . F HI N, OFOZ A. ZMAX ,FMAX, DFDZB 
EQUIVALENCE lARGL (7 » ,W) . f ARGL (8 ) ,T) , (ARGLUO K GAPFL) , 

1 (ARGH15J ,THEtA), (ARGL116) ,OTHK) 

DIMENSION C(5) 

IF fTRIG.GE.2) GO TO 10 0 

CONSTANTS 

C(l) = 0.70912 ♦0.13155*ALOG(T) 41.13475*AL0G(M) ♦0.ia345*ALOG(G 7PF 
1 Li 

C(2) =-1.2942 -1.13475 -0 . 07728* ( t ALOG ( VO ) **2 ) 

C(3i =-0.52376 
C(4) = 0.8671 
C(5) =-0.1835 


ZMIN = 0.04 
ZMAX = 3.5 

QRMIN = QRF1601(C,ZMIN) 

FMIN = ALOGCQRMIN) 

QRMAX = QRF1601 (C,ZMAX| 

FMAX = ALOG(QRHAXl 

OFOZA = Ct2)/ZHIN ♦ 2 . 0 *C ( 3 )* ALOG (ZMINi /Z M N 

OFDZB= C(2l/ZMAX ♦ 2 . 0*C ( 3 ) ♦ ALOG ( ZMAX ) /ZMAX 

RETURN 

100 Q = QRF1601(C,ZI 
RETURN 
F.NO 


FUNCTION QRF1601(C,Z) 

C 

C EQ 16-1 

C 

OTMEN5ICN CCl) 

COMMON /COMGAP/ ARGLC25), ZMI N, FMIN, OFOZA , ZMAX, FMAX , OFDZB 
EQUIVALENCE (ARGH15) , THETA) , ( A KGL C 16 ) , DTHK ) 

AS = ALOG(Z) 

XLNQR = cm ♦ (C(2)*C (3)*ABi*A8 ♦ C ( 4 i • ALOG ( THET A) ♦ C(5)*AL0G( 
I DTHK) 

QRF16D1=EXP<XLNQR) 

RETURN 

END 
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c 

c 


SUBROUTINE EQ1602 CITRI G , Z ♦0» 

H/HE, LONG IN-LINE GAP, L AMINAR ,GAMMA=0 (EQ16-2) 

COMMON /COMGAP/ ARC L C25 ) , ZMIN ,F MIN , DFOZ A. ZMAX .FMA X, OFOZ E 
EQUIVALENCE (ARGLIZ).MI, (ARGL(fil,T), lARGLdOl, GAPFL), 
1 IARGLI15) ,THETA), CARGL C16) ,0THK) 

OIMENSION CI5I 

IFCITRIG.EQ.2) GO TO 100 

CONSTANTS 

C(ll = 0.58794 ♦1.00779*ALOGIM)-0.19332*ALOG(GAPFL) 

1 *0. 13326»AL0G(T) 

C(2) =-0.91294 -1.00779 -0 .06665* (C ALOG (H) ) ••2 ) 

C(3) = 0.86146 
C(4> = -0.19332 
C(5> = 0.11529 

ZMIN =0.04 
ZMAX = 3.5 

QRHIN= QRF1602(C,ZMIN) 

FMIN = ALOGCQRMIN) 

QRMAX= QRF1602(C, ZMAX) 

PMAX = ALOGCQRHAX) 

OFOZA = C(2)/ZMIN ♦ 3. 0 *C <5 » ♦ ( ( ALOG ( Z MI N» ) **3 ) 

OFOZB = C(2>/ZNAX ♦ 3 .0 «C (51* ( f ALOGC ZMAX) )*«3 I 

RETURN 


100 Q =QRF16D2(C,Z) 
RETURN 
END 


FUNCTION QRF1602(C,Z) 


EQ 16-2 

DIMENSION- C(2) 

COMMON/COMGAP/ AR6Lf25) , ZMIN.FMIN, OFOZA. ZMAX, FMAX .OFOZB 
EQUIVALENCE ( ARGL (15 ) , THET A) ,(ARGL(16) ,OTHK) 

AB = ALOG(Z) 

XLNQR = C(l) ♦ C(2)*AB ♦ C ( 3) * A LOG ( THET A ) ♦ C (4 )* ALOG (OTHK) 
1 ♦ C(5)*AB*AB*AB 

QRF1602 = EXP(XLNQR) 

RETURN 

ENO 
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^ RS! GAP HEATING ANALYSIS - II 

I VOLUME I 


SUBROUTINE EQ16D34 ( ITRI G , 2 » Q) 


H/HFP, IN-LINE GAP, LAN INAR , 6ANHA UP TO 15 

IGOR = le EQ. 16-3, SIDE A AND 0 

ICOR = 19 — EQ. 16-4, SIDE A 

COMMON /COMGAP/ ARGL(25), ZMIN, FMIN, DFDZA, ZMAX, FMAX, DFOZ0 

EQUIVALENCE CARGL(3», ICORI, (ARSL(7). HI, (ARGL(IO), GAPFL), 

1 (ARGL (15I,THETA), ( A RGL (16 1 , OTHK) , (ARGL(8) ,TI 

DIMENSION C15I 
IFdTRIG.GE.ZI GO TO IOC 
AB - ALOG(H) 

IF(ICOR.£Q.19l GO TO 10 

CONSTANTS FOR EQ . 16-3 

C(l) = 1.88198 ♦2.86909*M ♦0. 13 326*ALOG (GAPFL I ♦ 0 . 451 39*AB* AB 
1 ♦.10532*ALOG(T) ♦ C.60717* (LOG(THETA) 

C(2) =-2.86909*1.33939*AB 

C(3) =-0.44674 

C(4I = 0.00376 

C(5) = DTHK/(GAPFL*HI 

GO TO 20 

CONSTANTS FOR EQ. 16-4 

10 C(l) = 2. 65299 ♦ 2.98795*AB «-0 . 4789*AB* AB ♦0 . 69f-61*AL0G (THET A I 
1 *0. 10187*ALOG(T) 

C(2) =-2.65299 ♦1.49684*AB 

C(3) =-0.50150 

C(4) = 0.00215 

C(5I = 0THK/(GAPFL*HI 

20 ZMIN = 0.04 
ZMAX = 3.8 

QRMIN= QR16D34(C,ZMIN) 

FMIN = AL0G(QRMINI 
QRMAX= QR16034(C, ZKAX) 

FMAX = ALOG(QRMAX) 

OFOZA= C(2)/ZMIN ♦ 2. 0*C ( 3) ♦ALOG ( ZMIN I /ZMI N *■ 

1 3.0*C(4)*( ( ALOG( C(5)*ZMAX) )••2)/ (C(5) *ZMIN, 

RETURN 

100 Q = QRie034(C,ZI 
RETURN 
END 


FUNCTION QR16D34 (C,Z) 

EQ.16-3 AND 63.16-4 

OIMENSICN C(5I 
A8 = AlOG(ZI 

ALNQ = Cd) ♦(C(2)> C(3)*ABI*AB ♦ C ( 4 I * d ALOG ( C (5 ) • Z ) I • »3 ) 
QR16034 = EXP(ALNQ) 

RETURN 

END 
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VOLUME I 


SUBROUTINE E Q16 056 ( ITR I G, Z ,Q) 

H/ME, IN-LINE GAP, LAMIKAR, GAMHA UP TO 15 


ICOR - 20 EQ. 16-5 

ICOR = 21 EQ. 16-6 


SIDES A AND B 
SIDE A 


COMMON /COMGAP/ ARGL<25I, ZHIN, FHIN, DFOZA, ZMAX, FMAX, OFQZB 


EQUIVALENCE (ARGL( 3), ICOR 
1 ,(ARGL( 6), gamma 

OIMENSICN C(3) 

IF (ITRIG.GE.2IGO TO IQO 

AB = ALOG(H) 

8C = ALOG(TMETA) 

IF(ICOR.EQ.21) GO TO 10 

— — - CONSTANTS FOR EQ . 16-5 


), (ARGL(7),M ), (ARGL (15' , THETA) 

), (ARGLI14) ,REPM) 


C 

c 


c 

c 


C(l) = 4.16869 ♦1.20204*AB ♦1.39857*BC -0 . 39913*ALOG ( REPM/l , 0E+ 06) 
1 ♦ (2.64364- 8. 30287*6AMMA)*6AHHA 

C(2) =-1 .20204 ♦0.54123*AB ♦O . 0 00 26* AB* AB 
C(3) = 0 .9374 

GO TO 2C 

CONSTANTS FOR EQ . 16-6 


10 CONTINUE 

C(l) = 1.48434 
1 *(2.4384 

C(2) =-1 .2492 
C(3) = 0 .10130 


♦1.2492*AB ♦ 0.74418*00 
- 7.964*6AMMA) *GAMMA 
♦0. 5469O*A0 


20 ZMIN = C .04 
ZMAX = 3 .8 

QRMIN= 0R16056(C,Z) 

FMIN = A lOG(QRMIN) 

QRMAX= CR16056(C,Z) 

FMAX = ALOG(QRMAX) 

DFOZA = C(2)/ZMIN ♦ 3. 0 *C (3)* ( ( ALOG (ZMI N) ) **2 ) /ZHIN 
OFDZB = C(2)/ZMAX t 3. 0 *C ( 3) • ( ( ALOG ( ZMAX ) ) **2 ) /ZM AX 

RETURN 

100 Q = QR16D56(C,Z) 

RETURN 

END 


C 

C 


FUNCTION QR16D56(C,Z) 


16-5 AND 16-6 


CORRELATION EQUATIONS 
A0 = ALOG(Z) 

ALNQR = C(1)*(C(2) *C(3) *Ae*AB)*AB 
QRlo056 = EXP(ALNQR) 

RETURN 

END 
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RSI GAP HEATING ANALYSIS - 

I VOLUME I 


REPORT MDC E124B 
JSC' 0V6SI 


SUBROUTINE EQ17 12 3 < ITRI G. Z , Y, Q I 
TRANSVERSE GAP WITH FORWARD FACING STEP 


21 

laminar 

EQ. 

17-1 

22 

TRANSITION 

EQ. 

l?-2 

17-3 

23 

TURBULENT 

EQ. 


(DOWN STREAM WALU 


COMMON /COMGAP/ ARGL(25), ZMIN, FMIN, DFDZA, 7MAX, FMAX, DFDZB 

EQUIVALENCE (ARGL( 7),W ),(ARGL(9), D)» 

1 (ARGL( 1),IBLI 

0 = STEP HEIGHT (CM » 

W = GAP WIDTH (CMI 

0IMENSI0^ C(3) 

IF (ITRIG .GE.2I GO TO 100 

IFdBL.NE.ll GO TO 10 

LAMINAR FLOW CONSTANTS 

C(l» =-0.19252 ♦(S. 67183-11. 37065 *Oji*D 
C(2) =-t*. 54537 ♦23.5240 7*W*0 -0 . 00 192*Y*Y*Z 
C(3> = 1. 26103 -6.57773*W*0 
GO TO 30 

10 IF(I8L.N£. 2) GO TO 20 

TRANSITIONAL FLOW CONSTANTS 

C(l) =-0.65828 ♦ 8.93015*0 - 2.63955*0*0 ♦ 3.18C26*W*W 
C(2) =-4.00086 

C(3) = 1.16511 - 0.Ci083*Y*Y*W *0.01461*Y 
GO TO 30 

TURBULENT "LOW CONSTANTS 

20 C(l) =-0.81624 ♦7.58098*0 ♦0.17100*Y -0.28504*Y*W ♦1.437l2*W*W 
C(2) =-2.5465 

C(3» = 0.35211 - 1,11972*0*0 

30 CONTINUE 

ZMIn =-1.0E*30 
7MAX = 3.8 

QMAX = QR17123(C, ZMAX) 

FMAX = ALOGIQMAXJ 

RETURN 

100 Q = QR17123(C.Z ) 

RETURN 

ENO 


FUNCTION QR17123(CiZ) 

STEP HEATING CORKELATIN EQ.17-1 

DIMENSION C(3» 

QR17123 = EXP(C(1) ♦(C(2)tC 3)*Z»*Z) 
RETURN 


EQ.17-1 , EQ.17-2, EQ.17-3 
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RSI GAP HEATING ANALYSIS - II 

O 


VOLUME I 


APPENDIX D 

The units used In this volume are SI Units, except Section 2 where the raw 
data presented was In English Units. Conversion factors required for used 
herein are given In the following table: 


Physical quantity 

U.S. 

Customary Unit 

Conversion 

factor 

(*) 

SI Unit 

Convective Keating Rate 

2 

Btu/ft -sec 

1.134 X lO'^ 

2 

Watts/m 

Convective Heat Transfer Coef- 
ficient Based on Taw/T^ “ 0.895 

2 

Ibm/ft -sec 

4.88 

Kg/m^-S 

Enthalpy 

Btu/lbm 

2.324 X 10^ 

J/Kg 

Heat Transfer Coefficient 

Btu/ft^-sec-°F 

2.042 X 10^ 

Watts/m^-°C 

Lf'ngth 

In. 

2.54 X 10“^ 

m 

Length 

ft. 

0.305 

m 

Pressure 

Ibf/ft^ 

47.88 

N/m^ 

Pressure 

Ibf/ft^ 

47.88 X 10^ 

KN/m^ 

Pressure 

Ibf/ft^ 

47.88 X 10^ 

MN/m^ 

Reynolds Number 

Re/ ft 

3.28 

Re/m 

(RHO V)1 Freestream Density - 
Velocity Product 

2 

Ibm/ ft -sec 

4.88 

Kg/m^-S 

Temperature 

°F 

5/9(°F-32) 

°C 

Temperature 

°F 

5/9(°F+460) 

°K 


*Multiply value In U.S. Customary Unit by conversion factor to obtain equivalent 
value in SI Unit. 
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